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1 3.  ABSTRACT  (M sjttmum  200  words) 

A  three-year  Balanced  Technology  Initiative  (BTI)  program  to  develop  tests  for  determining  the  failure  criteria  of 
ceramics  under  ballistic  impact  is  described.  The  phenomenology  of  the  penetration  of  a  confined  ceramic  block  by  a  long  rod 
pcnctrator  was  established  by  examining  tlie  fracture  damage  in  targets  alter  ballistic  testing.  The  material  properties  that 
govern  the  penetration  resistance  of  ceramics  were  deduced.  Two  plate  impact  tests — the  high-strain-ratc  press urc/s hear  test 
and  the  symmetric  pressurc/shcar  lest— were  then  developed  to  measure  the  granulation  and  granular  flow  properties  of 
candidate  annor  ceramics  at  strain  rates  and  pressures  representative  of  those  at  the  tip  of  an  advancing  pcnctrator.  In  the 
development  of  the  tests,  a  laser-Dopplcr  vclocimeter  was  applied  for  the  first  time  to  measure  the  transverse  motion  of  the 
target  plate.  The  tests  were  applied  to  aluminum  nitride,  boron  carbide,  and  Coors  AD-85  and  AD-995  aluminas;  and  the 
resultant  behaviors  were  compared.  Two  one -dimensional  computational  models,  one  based  on  a  continuum  Mohr-Coulomb 
description  of  the  flow  of  granulated  ceramic  and  the  other  (developed  in  a  companion  BTI  project)  based  on  the  frictional 
sliding  of  blocks,  were  applied  to  interpret  the  symmetric  pressurc/shcar  tests,  which  revealed  the  near-penetrator  failure 
behavior.  The  observations  and  data  were  used  as  input  for  two-dimensional  penetration  simulations  using  the  block-sliding 
model. 
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summary 


The  objectives  of  this  three-year  Balanced  Technology  Initiative  program  were  to  identify 
the  material  properties  that  govern  the  resistance  of  ceramics  to  penetration,  develop  tests  for 
measuring  these  properties,  and  apply  the  tests  to  evaluate  several  ceramics  that  are  candidates 
for  armor  application.  This  information  is  needed  to  design  armor  structures  that  can  protect 
against  advanced  penetrators. 

To  identify  die  material  properties  governing  penetration  resistance,  we  needed  to 
understand  how  ceramics  resist  penetration  and  how  they  fail.  We  sought  this  informauon  by 
examining  ceramic  armor  targets  that  had  been  penetrated  or  partially  penetrated  and  reached  the 
following  conclusions: 

•  Penetration  of  a  confined  ceramic  block  by  an  impacting  long  rod  projectile 
occurs  by  the  comminudon  of  a  zone  of  the  ceramic  at  the  leading  edge  of  the 
projecule  and  the  flow  of  the  fine  fragments  out  of  the  projectile  path.  (We 
named  this  zone  the  Mescall  zone  in  recognition  of  John  Mescall,  who 
postulated  its  existence  based  on  his  hydrocode  calculations  in  the  1980s.) 

•  The  material  properties  governing  the  penetration  resistance  of  a  confined 
ceramic  block  include  the  dynamic  compressive  strength,  the  comminution  and 
granular  flow  resistance  under  high  rate  and  high  pressure  conditions,  and  the 
abrasiveness  of  the  ceramic  fragments. 


To  develop  tests  for  measuring  these  properties,  we  sought  ways  to  produce  conditions  of 
stress,  strain,  and  strain  rate  representative  of  those  produced  during  penetration.  By  performing 
hydrocode  calculations  of  simulated  penetration  scenarios,  we  found  that  stresses  of  several  tens 
of  kilobars  and  strain  rates  of  1 05/s  were  typical.  We  then  analyzed  several  types  of  tests  that 
produce  appropriate  conditions  to  determine  their  suitability  for  providing  the  granulation  and 
flow  properties  of  candidate  armor  ceramics.  We  chose  three  tests  as  most  promising:  (1)  high 
strain-rate  pressure/shear  plate  impact,  (2)  symmetric  pressure/shear  plate  impact,  and  (3) 
explosively  loaded  contained  rod  impact. 

We  performed  the  high  strain-rate  pressure/shear  (HSRP/S)  plate  impact  test  by  t  For 


impacting  a  thin  specimen  plate  against  a  hard  anvil  plate  inclined  at  an  angle.  This  thin 
specimen  plate  was  backed  up  by  a  thick,  hard  flyer  plate,  which  created  a  state  of  high  pressure 


and  high  shear  strain  rate  in  the  specimen.  The  flyer  and  anvil  plates  were  designed  to  remain 
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elastic  while  the  specimen  deformed.  During  the  early  stages  of  the  test,  before  release  waves 
reach  the  center  of  the  plate  from  the  periphery,  stress  components  in  die  specimen  were  inferred 
from  the  projectile  velocity  and  particle  velocities  at  the  rear  surface  of  the  anvil.  Thus, 
measuring  particle  velocity  history  at  the  rear  surface  allowed  us  to  construct  a  shear  stress-shear 
strain  curve.  The  HSRP/S  test  was  adapted  to  test  aluminum  nitride  by  substituting  boron 
carbide  for  the  usual  steel  flyers  and  anvils.  The  transverse  particle  velocity  at  the  rear  surface  of 
the  anvil  was  monitored  by  using  a  transverse  displacement  interferometer,  which  uses  the 
interference  of  light  beams  of  different  Doppler-shifted  frequencies  to  indicate  the  transverse 
velocity. 

The  HSRP/S  test  could  not  be  used  for  specimens  of  the  very  strongest  ceramics,  such  as 
SiC,  B4C,  and  T1B2,  because  even  stronger  materials  such  as  diamond  or  cubic  boron  nitride  are 
required  for  the  flyers  and  anvil  plates.  The  second  test  method,  the  symmetrir  pressure/shear 
(SP/S)  test,  circumvented  this  requirement  by  impacting  a  ceramic  anvil  in  a  pressure/shear 
configuration  with  a  flyer  composed  entirely  of  the  same  ceramic.  In  this  test  the  impact  velocity 
and  angle  and  the  histories  of  normal  and  transverse  particle  velocity  at  the  target  rear  surface 
were  measured.  These  histories  were  compared  with  those  from  a  numerical  simulation  using  a 
postulated  material  model,  which  was  iteratively  modified  until  the  experiment  and  the 
simulation  agree.  If  the  match  is  good,  the  model  is  at  least  valid  for  the  test  conditions.  The 
SP/S  test  is  not  so  easy  to  analyze  as  the  HSRP/S  test  and  provides  the  shear  stress-strain  curve 
only  indirectly,  but  it  nevertheless  allows  very  strong  ceramics  to  be  evaluated  under  conditions 
of  high  confining  pressure  and  strain  rate  representative  of  armor  penetration. 

The  SP/S  test  was  applied  to  aluminum  nitride,  boron  carbide,  and  Coors  AD-85  and  AD- 
995  aluminas.  In  our  adaptation  of  the  test,  the  transverse  rear-surface  particle  velocity  of  the 
target  plates  was  measured  by  using  a  laser-Doppler  velocimeter  (LDV)  rather  than  a  transverse 
displacement  interferometer.  To  our  knowledge,  this  is  the  first  application  of  LDV  in  shock 
physics  experiments.  The  LDV  has  the  advantage  of  not  requiring  a  diffraction  grating  on  the 
rear  surface  of  the  target;  a  diffraction  grating  is  difficult  and  expensive  to  etch  onto  ceramic,  but 
it  is  required  (and  must  survive  shock  loading)  when  using  the  traverse  displacement 
interferometer. 

We  tried  several  rod  impact  configurations  in  an  attempt  10  develop  a  test  simpler  than 
the  pressure/shear  tests.  Ceramic  rods  were  surrounded  by  metal  or  plastic  containment  and 
impacted  axially  with  explosively  driven  flyer  plates.  The  confinement  successfully  elevated  the 
mean  stress  in  the  deforming  ceramic.  Unfortunately,  we  were  unsuccessful  in  recovering  the 
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rods  or  containment  to  measure  posttest  deformation.  Although  some  contained-rod  impact 
configuration  could  be  made  to  work,  we  did  not  pursue  this  problem  further. 

Both  types  of  pressure/shear  experiments  revealed  a  wide  range  of  mechanical  response 
among  the  ceramics  tested.  The  transverse  velocity  for  AIN  nearly  reached  the  elastic  level, 
wb  eas  the  velocity  for  B4C  was  only  35%  of  the  elastic  level.  For  AIN  and  B4C,  velocities 
decayed  with  time;  for  AD-85  and  AD-995  AI2CV3,  velocities  increased  with  time.  The  former 
behavior  is  characteristic  of  faster  decay  waves  overtaking  the  initial  shear  wave  from  behind, 
whereas  the  latter  behavior  is  characteristic  of  slower  following  waves.  The  different  ceramics 
seem  to  have  different  mechanisms  for  failure  and  postfailure  flow,  and  these  mechanisms 
operate  at  different  threshold  loads. 

Two  one-dimensional  computational  models  were  developed  to  infer  strength  behavior 
from  the  symmetric  pressure/shear  particle  velocity  records.  One  model  is  based  on  a  continuum 
Mohr-Coulomb  description  of  the  flow  of  granulated  ceramic,  and  the  other  is  based  on  frictional 
sliding  of  blocks.  With  the  exception  of  AIN  tested  in  SP/S,  the  shear  strengths  inferred  from  the 
pressure/shear  tests  were  lower  than  the  strengths  reported  by  others.  We  think  the  reason  is  that 
we  have  essentially  measured  the  response  of  comminuted  ceramic,  which  may  be  appropriate 
for  evaluating  materials  for  armor  applications.  The  longitudinal  wave,  which  is  strong  enough 
to  cause  damage  even  in  compression,  runs  ahead  of  the  shear  wave.  The  shear  wave,  which  we 
measured,  is  always  running  through  damaged  material. 

The  Mohr-Coulomb  simulations  reproduced  the  transverse  velocity  histories  at  least 
qualitatively  and,  in  the  cases  of  AIN  tested  in  HSRP/S  and  AD-85  and  AD-995  AI2O3  tested  in 
SP/S,  quantitatively.  The  simulations  were  able  to  predict  the  SP/S  response  of  AD-995  from  the 
SP/S  response  of  AD-85,  but  were  unable  to  predict  the  SP/S  response  of  AIN  from  the  HSRP/S 
results.  We  concluded  that  the  results  of  the  computations  are  valid  at  least  ^perhaps  only)  for 
the  conditions  of  a  given  test.  Further  work  is  needed  to  refine  these  models. 

We  established  at  the  outset  that  any  tests  for  measuring  the  ballistic  behavior  of  armor 
ceramics  must  measure  their  behavior  under  penetration  conditions.  Both  the  HSRP/S  and  SP/S 
tests  have  met  this  requirement,  because  mean  stresses  and  strain  rates  characteristic  of  the 
Mescall  zone  were  achieved  in  both.  Thus  we  conclude  that  pressure/shear  tests  in  these  two 
configurations  with  rear  surface  laser  diagnostics  and  with  appropriate  constitutive  analysis  are 
suitable  for  measuring  the  high  pressure,  high  strain  rate  comminution  and  granular  flow 
properties  of  armor  ceramics. 
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Future  work  should  focus  on  resolving  certain  measurement  issues  in  the  tests  and 
refining  the  simulation  models.  Further  work  is  also  needed  to  design  a  simple,  inexpensive  test 
for  examining  ceramics  under  penetration  conditions.  Although  the  pressure/shear  tests 
described  here  were  successful,  they  require  a  sophisticated  gas  gun  facility  with  accompanying 
interferometers  and  high-speed  recording  equipment.  We  recommend  that  the  search  for  a 
simpler  test  begin  with  investigation  of  the  cylindrical  cavity  expansion  test.  A  simpler  static  test 
may  be  possible,  if  a  large  enough  mean  stress  can  be  applied,  but  a  dynamic  test  may  be 
necessary  to  produce  such  a  large  mean  stress. 
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Chapter  1 
INTRODUCTION 


Armors  that  incorporate  ceramics  are  known  to  provide  higher  resistance  to  penetrators 
than  metallic  armors.  However,  designing  armor  packages  having  optimal  performance  is 
difficult  because  the  reasons  for  the  enhanced  ballistic  protection  are  only  partly  understood.  To 
attain  the  information  necessary  to  design  superior  armor,  new  tests  are  needed  that  measure  the 
ceramic  properties  that  govern  penetration  resistance.  Simple  ballistic  tests  can  only  rank  the 
various  ceramics  and  armor  configurations;  they  do  not  provide  the  information  necessary  to 
guide  the  development  of  ceramic  microstmctures  or  armor  geometries. 

OBJECTIVES 

The  objectives  of  this  project  were  to  identify  the  material  properties  that  govern  the 
resistance  of  ceramics  to  penetration,  develop  tests  for  measuring  these  properties,  and  apply  the 
tests  to  evaluate  several  candidate  ceramics. 

APPROACH 

As  the  first  step,  we  sought  to  understand  the  phenomenology  of  failure  under  ballistic 
penetration  conditions.  After  reviewing  the  classified  and  unclassified  literature  to  assess  the 
current  knowledge,  we  obtained  penetrated  and  partially  penetrated  ceramic  armor  targets  (in 
conjunction  with  a  related,  parallel  project1)  and  sectioned  and  examined  them  to  determine  the 
mechanisms  of  penetration.  We  then  performed  hydrocode  calculations  of  ceramic  armor 
penetration  to  determine  the  conditions  of  stress,  strain,  and  strain  rate  produced  during 
penetration.  These  calculations  were  necessarily  qualitative  because  accurate  material  models,  a 
goal  to  be  attained  at  the  end  of  the  project,  were  unavailable. 

Once  the  mechanisms  and  conditions  of  failure  and  postfailure  flow  were  determined,  we 
began  to  consider  test  techniques  that  could  invoke  the  salient  material  failure  mechanisms  under 
representative  load  conditions  in  well-controlled  ways.  We  further  considered  diagnostic 
techniques  that  could  measure  the  governing  stress,  strain,  and  strain  rate  variables.  Given  the 
high  deformation  rates  involved  in  penetration,  our  attention  naturally  focused  on  impact  and 
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explosive  loading  tests.  We  analyzed  several  candidate  tests,  including  normal  and  slanted  plate 
impact,  rod  impact,  and  cylindrical  cavity  expansion. 

We  chose  three  candidate  tests  as  the  most  promising:  (1)  high  strain-rate  pressure/shear 
plate  impact,  (2)  symmetric  pressure/shear  plate  impact,  and  (3)  explosively  loaded  contained  rod 
impact.  We  implemented  each  test  and,  from  the  results,  decided  whether  to  pursue  each 
technique  further.  We  had  the  most  success  with  the  pressure/shear  tests  and  so  applied  most  of 
our  developmental  effort  to  them. 

In  Chapter  2,  we  present  a  qualitative  and  quantitative  picture  of  the  conditions  ahead  of 
an  advancing  penetrator.  In  Chapter  3,  we  describe  the  high  strain-rate  and  symmetric 
pressure/shear  techniques  and  the  analyses  used  to  extract  the  data.  Next,  we  present  the  results 
of  tests  on  aluminum  nitride  (AIN),  two  grades  of  alumina  (AI2O3),  and  boron  carbide  (B4C), 
followed  by  a  comparison  of  these  results  with  those  obtained  by  others  using  different  testing 
techniques  and  with  results  obtained  for  silicon  carbide  (SiC).  In  Chapter  4,  we  describe  the 
other  candidate  tests  considered  and  present  simple  analyses  of  them.  We  also  describe  a  brief 
investigation  we  undertook  on  penetrator  erosion  phenomenology.  Conclusions  and 
recommendations  are  presented  in  Chapter  5.  The  Appendix  contains  copies  of  papers  we 
published  on  this  project,  along  with  a  list  of  presentations. 
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Chapter  2 

PENETRATION  PHENOMENOLOGY 


QUALITATIVE  PENETRATION  CONDITIONS 

When  an  approaching  penetrator  initially  contacts  a  block  of  ceramic,  a  dynamic  Hertzian 
contact  stress  wave  field  is  generated.  The  material’s  ability  to  resist  this  initial  pan  of  the  attack 
is  therefore  governed  by  its  dynamic  hardness.2  If  the  penetrator  can  apply  more  stress  than  the 
strength  associated  with  the  dynamic  hardness  of  the  target,  it  will  advance  into  the  ceramic. 
Further  advance  requires  that  the  ceramic  be  moved  out  of  the  penetrator  path.  In  metal  armor, 
most  of  this  movement  is  accommodated  by  dislocation  plasticity,  but  in  ceramics,  although  the 
initial  movement  may  be  accompanied  by  dislocation  plasticity,3  too  few  slip  systems  are 
available  for  sustained  plasticity.  Therefore,  the  large  deformations  required  to  move  material 
out  of  the  penetrator  path  are  accommodated  by  the  sliding  and  flow  of  microscopic  fragments.2 
Thus,  the  material's  ability  to  resist  granulation  is  an  important  material  property,  as  is  its  ability 
to  resist  granular  flow. 

For  continued  penetration  through  a  thick,  well-confined  block  of  ceramic,  the  material 
ahead  of  the  penetrator  must  be  moved  laterally,  then  opposite  the  direction  of  attack.  In  steady- 
state  penetration,  the  fine  ceramic  fragments  produced  in  the  "Mescall"  zone  at  the  tip  of  the 
penetrator4'6  (Figure  1)  flow  out  of  the  penetrator’s  path.  To  an  observer  traveling  with  the 
interface  between  the  penetrator  and  the  ceramic  (which  may  not  be  a  sharp  boundary),  a  stream 
of  comminuted  ceramic  appears  to  approach,  then  divide  and  byp.iss  the  observer's  position. 
Likewise,  penetrator  material  appears  to  approach  from  behind,  divide,  and  reverse  direction  as 
the  penetrator  is  eroded.  The  deformation  of  penetration  is  expected  to  occur  under  conditions  of 
high  mean  stress  and  high  strain  rates.  Any  tests  for  measuring  the  failure  criteria  of  ceramic 
under  ballistic  penetration  must  exercise  the  material  under  such  conditions.  In  this  work,  we 
concentrated  on  the  material  properties  that  govern  the  steady-state  phase  of  penetration,  namely, 
granulation  and  granular  flow  properties. 
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Figure  1 .  Cracking  pattern  in  ceramic  targets  during  the  steady-state  phase  of 
the  penetration  process. 
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QUANTITATIVE  PENETRATION  CONDITIONS 


To  fail  and  deform  well-confined  ceramic  requires  that  a  stress  on  the  order  of  that 
required  for  yielding  in  uniaxial  strain  compression  (i.e.,  the  Hugoniot  elastic  limit  (HEL)]  be 
exceeded.  To  obtain  a  quantitative  estimate  of  the  conditions  ahead  of  a  penetrator,  we  used  the 
model  of  Curran  et  al.  for  the  motion  of  granular  beds4  in  Cooper's  two-dimensional  Lagrangian 
hydrocode7  and  simulated  the  penetration  of  a  1.6-km/s  tungsten  rod  into  a  block  of  aluminum 
nitride.  The  resulting  contour  plots  of  mean  stress  and  effective  strain  rate  are  shown  in  Figure  2. 
Essentially,  a  material  element  ahead  of  a  penetrator  is  deformed  at  a  strain  rate  greater  than 
105  s'1  while  under  a  mean  stress  of  order  of  the  HEL,  9  GPa8  in  this  case.  The  results  in 
Figure  2  are  not  strongly  dependent  on  the  details  of  the  tungsten  and  ceramic  material  models. 

The  details  of  the  loading  path  for  a  material  element  are  complex,  and  complex  loading 
paths  induce  complex  responses,  such  as  anisotropy,  in  the  ceramic.  In  this  work,  the  effects  of 
complex  loading  paths  are  considered  secondary.  We  attempted  to  develop  tests  that  could 
duplicate  simply  the  mean  stress,  strain  rate,  and  strain  conditions  of  penetration. 
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Figure  2.  Simulation  of  1.6-km/s  tungsten  rod  penetrating  confined  A^N.  (concluded) 
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Chapter  3 

PRESSURE/SHEAR  TESTS 


Two  types  of  pressure/shear  tests  were  conducted:  high-strain-rate  pressure/shear 
(HSRP/S)  tests  and  symmetric  prcssure/shear  (SP/S)  tests.  The  HSRP/S  experiments  were 
performed  on  AIN,  and  the  SP/S  experiments  were  performed  on  AIN,  two  grades  of  AI2O3,  and 
B4C.  The  experimental  results  were  compared  with  those  obtained  for  SiC,  with  results  obtained 
by  other  workers,  and  with  ballistic  performance. 

GENERAL  EXPERIMENTAL  SETUP 

The  plate  impact  experiments  were  performed  at  SRrs  inch-bore  gas  gun  facility 
(Figure  3).  This  gun  is  capable  of  launching  230-g  projectiles  to  1  km/s.  It  has  two  breech 
mechanisms:  a  wraparound  type  for  lower  velocities  and  a  rapid-opening-valve  type  for  high 
velocities.  The  breech  mechanisms  are  described  by  Erlich9  in  a  publication  included  in  the 
Appendix.  A  key  on  the  side  of  the  projectile  engages  a  slot  broached  in  the  gun  barrel  and 
prevents  the  projectile  from  rotating,  thus  allowing  oblique/parallel  impact  experiments  to  be 
performed.  The  velocity  of  the  projectile  is  measured  just  before  impact  by  recording  the  contact 
times  of  four  wires  placed  in  its  path.  A  flyer  plate  in  the  form  of  a  disk  is  carried  on  the  nose  of 
the  projectile.  A  target  or  anvil  plate,  also  a  disk,  is  situated  in  an  evacuated  tank  at  the  gun’s 
muzzle.  The  tank  and  barrel  are  evacuated  to  below  100  mTorr  to  prevent  the  compression  of  a 
gas  column  ahead  of  the  projectile,  which  would  disturb  the  anvil  and  lubricate  the  impact 
interface.  The  HSRP/S  projectile  and  the  target  assemblies  are  shown  in  Figure  4. 

In  the  experiments  reported  here,  the  diameters  and  thicknesses  of  the  impacting  ceramic 
plates  were  chosen  so  that  the  First  unloading  wave  to  arrive  at  the  center  was  from  the  anvil  rear 
surface.  The  release  waves  from  the  periphery  arrived  slightly  later  and  that  from  the  flyer  free 
surface  still  later.  The  flyer,  anvil,  and  (in  the  case  of  HSRP/S)  specimen  plates  were  made 
flatter  than  0.6  p.m  over  a  40-mm  diameter  to  ensure  wave  planarity. 
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Four  pairs  of  electrical  contact  pins  were  located  near  the  edge  of  the  anvil  at  the  3,  6,  9, 
and  12  o'clock  positions  (Figure  4).  Single  gold  contact  pins  were  evaporated  onto  the  anvil; 
pairs  of  contact  pins  were  evaporated  onto  the  specimen  or  flyer  impact  face,  aligned  so  that  they 
would  short  across  the  paired  pins  and  pads  on  the  anvil  at  impact.  The  times  of  electrical 
contact  indicated  the  "tilt"  or  deviation  from  simultaneity  (lack  of  parallelism)  of  the  impact 
across  the  face  of  the  target  plate.  When  the  anvil  was  placed  in  the  gun  target  chamber,  it  was 
aligned  parallel  to  the  face  of  the  projectile  within  0.08  mrad,  by  using  either  a  newly  developed 
interferometric  technique10  or  the  optical  technique  of  Kumar  and  Clifton.11 

Signals  from  the  projectile  velocity  contact  pins,  the  tilt  pins,  and  the  interferometer 
photodetector  output  were  recorded  on  a  4-channel,  500-megasample/s  transient  digitizer  and 
analog  oscilloscopes.  The  analog  data  were  reduced  by  magnifying  and  digitizing  the 
oscilloscope  photographs  and  processing  the  resulting  computer  files. 

HIGH-STRAIN-RATE  PRESSGRE/SHEAR  TEST 

High-strain-rate  pressure/shear  plate  impact  testing  is  a  relatively  new  technique12'14  that 
has  been  applied  to  metals  to  obtain  stress-strain  curves  at  shear  strain  rates  greater  than  105  s1 
under  high  and  controllable  levels  of  nearly  hydrostatic  pressure.  The  method  is  based  on 
concepts  drawn  from  plate  impac  testing  and  Hopkinson  bar  experiments.  From  the  former  is 
taken  the  use  of  plane  waves  generated  by  parallel  impact,  so  that  a  one-dimensional  wave  theory 
is  applicable  and  the  interpretation  of  experimental  results  is  simplified.  From  the  latter  comes 
the  concept  of  using  elastic  "probes”  to  sandwich  and  deform  a  thin  specimen  and  to  condu  t 
specimen  loads  and  velocities  to  a  remote  location  for  measurement. 

High-strain-rate  pressure/shear  impact  testing  is  performed  by  impacting  a  thin  specimen 
plate  against  a  hard  anvil  plate  inclined  at  an  angle,  9,  as  shown  in  Figure  5.  This  thin  specimen 
plate  is  backed  up  by  a  thick,  hard  flyer  plate,  which  creates  a  state  of  high  pressure  and  high 
shear  strain  rate  in  the  specimen.  The  flyer  and  anvil  plates  are  designed  to  remain  elastic  while 
the  specimen  deforms;  during  the  early  stages  of  the  test,  before  release  waves  reach  the  center  of 
the  plate  from  the  periphery,  stress  components  in  the  specimen  can  be  inferred  from  the 
projectile  velocity  and  panicle  velocities  at  the  rear  surface  of  the  anvil.  Thus,  measuring  a 
history  of  particle  velocity  at  the  rear  surface  allows  us  to  construct  a  shear  stress-shear  strain 
curve. 
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Diffraction  Grating 
(on  surface) 
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Figure  5.  Configuration  for  HSRP/S  experiment. 
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In  the  early  stages  of  the  test,  waves  reverberate  in  the  specimen  (Figure  6).  Because 
these  waves  are  plastic,  they  are  quickly  dispersed  and  the  stress  state  in  the  specimen  becomes 
nominally  homogeneous.  The  out-of-plane  shear  and  normal  stress  components  are  the  same  as 
those  in  the  elastic  plates.  Because  the  specimen's  elastic  resistance  to  volume  change  will  not 
allow  a  normal  difference  in  velocity  to  be  maintained  across  the  specimen,  the  homogeneous 
normal  stress  in  the  specimen  quickly  "rings  up"  to  the  normal  stress  that  would  exist  if  the  flyer 
and  anvil  impacted  in  the  absence  of  the  specimen.  Thus,  if  the  flyer  and  anvil  are  made  from 
the  same  material,  the  stress  component  normal  to  the  specimen  surface,  <T,  is  given  by 
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a=2PC\uo 


(1) 


where  p  is  the  density  of  the  elastic  plates,  c\  is  the  elastic  longitudinal  wave  speed,  and  uq  is  the 
component  of  projectile  velocity  normal  to  the  plates.  Computations  verify  that  the  hydrostatic 
pressure  in  the  specimen  is  nearly  equal  to  the  negative  of  the  normal  stress. 

A  finite  transverse  difference  in  velocity  can  be  maintained  across  the  specimen  for  the 
duration  of  the  experiment.  Thus,  the  shear  stress  need  not  ring  up  to  the  level  that  would  exist 
in  the  absence  of  the  specimen.  The  shear  stress  component,  r,  at  the  surface  of  the  specimen  is 
the  same  as  that  in  the  elastic  anvil,  which  is  related  to  the  transverse  particle  velocity,  v,  in  the 
anvil  by 


T=pC2V  (2) 

where  C2  is  the  elastic  shear  wave  speed  in  the  flyer  and  anvil.  Because  the  anvil  remains  elastic, 
v  is  exactly  half  the  velocity  at  the  rear,  free  surface  of  the  anvil,  Vfs.  Thus, 

r=lpc2V[s  (3) 

The  nominal  shear  strain  rate  in  the  specimen,  y  is  given  by  the  transverse  velocity 
difference  across  the  specimen  divided  by  its  thickness;  that  is. 


14 


RM6795-39 


Figure  6.  HSRP/S  l-x  diagram  showing  waves  reverberating  in 
the  specimen. 

Elastic  wavefronls  only  are  shown  in  the  specimen. 
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where  vp  and  va  are  the  flyer  and  anvil  particle  velocities  adjacent  to  the  specimen  and  h  is  the 
thickness.  If  the  flyer  and  anvil  remain  elastic,  vp  -  Va  =  Vg  -  v^,  where  Vo  is  the  component  of 
projectile  velocity  transverse  to  the  plates.  Thus, 


‘  -VQ-Vfe 
r~  h 


(5) 


By  integrating  Eq.  (5)  over  time,  we  obtain  the  nominal  shear  strain  in  the  specimen,  and 
combining  the  strains  with  the  stresses  from  Eq.  (3)  allows  us  to  construct  a  stress-strain  curve. 

In  the  HSRP/S  experiments,  the  transverse  particle  velocity  at  the  rear  surface  of  the  anvil 
was  monitored  by  using  a  transverse  displacement  interferometer  (TDI).  This  device,  developed 
by  Kim  et  al.,15  uses  the  interference  of  light  beams  of  different  Doppler-shifted  frequencies  to 
indicate  the  transverse  velocity.  The  TDI  is  shown  schematically  in  Figure  7.  The 
interferometer  consists  of  a  beam  of  laser  light  directed  at  normal  incidence  onto  a  diffraction 
grating  applied  to  the  rear  surface  of  the  anvil.  The  grating  is  oriented  so  that  the  diffracted 
beams,  which  are  spaced  symmetrically  about  the  surface  normal,  are  in  the  plane  of  the 
transverse  velocity  component  to  be  measured.  A  pair  of  diffracted  beams  with  equal  angles  on 
each  side  of  the  normal  are  selected  and  mixed  with  a  system  of  mirrors  and  a  beam-splitter  and 
projected  onto  a  photodetector.  When  the  rear  surface  of  the  anvil  plate  moves  transversely,  the 
two  selected  diffracted  beams  are  Doppler-shifted,  one  to  a  higher  frequency  and  one  to  a  lower. 
When  these  beams  are  mixed,  beating  occurs  at  a  frequency  proportional  to  the  transverse 
velocity.  The  transverse  displacement,  S,  per  beat  cycle  is 


where  d  is  the  pitch  of  the  diffraction  grating  and  n  is  the  diffraction  order  of  the  chosen  beams. 
We  used  first-order  (n  =  1 )  diffracted  beams  and  a  5-pm-pitch  grating  in  the  experiments 
reported  here.  The  TDI  is  insensitive  to  the  normal  motion  of  the  plate. 
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Figure  7.  Transverse  displacement  interferometer. 
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In  the  first  HSRP/S  experiment  on  AIN,  the  flyer  and  the  specimen  were  disks  55  mm  in 
diameter  with  thicknesses  of  5.3  and  0.58  mm,  respectively;  in  the  second  experiment  the 
specimen  thickness  was  0.60  mm.  The  anvils  were  disks  41  mm  in  diameter  and  4.8  mm  thick. 
The  flyers  and  the  anvils  were  made  from  hot-pressed  B4C,  and  the  specimens  were  AIN. 

The  specimen  was  glued  to  the  flyer  at  the  periphery,  and  the  flyer-specimen  assembly 
was  then  glued  to  the  nose  of  a  projectile  cut  at  an  angle  of  20  degrees.  Small  relief  holes  were 
drilled  in  the  hollow  projectile  to  allow  it  and  the  specimen -flyer  interface  to  be  evacuated  at  the 
same  time  as  the  breech,  banel,  and  target  chamber  of  the  gun. 

The  TDI  diffraction  grating  was  applied  to  the  anvil  by  vapor-depositing  a  platinum 
mirror  and  copying  a  master  200-lines/mm  (5-pm-pitch)  Ronchi  ruling  onto  the  polished  rear 
surface  with  photoresist.  The  anvil  plate  was  glued  into  the  center  of  a  larger  polymethyl 
methacrylate  (PMMA)  disk,  which  was  in  turn  supported  by  the  gun's  target  holder  mechanism. 

SYMMETRIC  PRESSURE/SHEAR  TEST 

Experimental  Technique 

A  key  requirement  in  the  HSRP/S  test  is  that  the  flyer  and  anvil  plates  that  sandwich  the 
specimen  not  yield.  This  requirement  cannot  be  met  for  specimens  of  the  very  strongest 
ceramics,  such  as  SiC,  B4C,  and  T1B2,  unless  diamond  or  cubic  boron  nitride  flyers  and  anvils 
are  used.  An  alternative  technique  must  therefore  be  considered.  A  good  alternative  is  the 
symmetric  pressure/shear  (SP/S)  test,1 1  in  which  an  anvil  is  impacted  in  a  pressure/shear 
configuration  with  a  flyer  of  the  same  material  (Figure  8).  While  not  nearly  so  easy  to  analyze  as 
the  high- strain-rate  pressure/shear  test,  the  symmetric  pressure/shear  test  can  still  provide  similar 
conditions  of  high  confining  pressure  and  strain  rate.  In  this  test,  the  impact  velocity  and  angle 
and  the  histories  of  normal  and  transverse  particle  velocity  at  the  target  rear  surface  are 
measured.  These  histories  are  compared  with  those  from  a  numerical  simulation  using  a 
postulated  material  model,  which  is  iteratively  modified  until  the  experiment  and  the  simulation 
agree.  If  the  match  is  good,  the  model  is  at  least  valid  for  the  test  conditions. 

At  impact  in  the  SP/S  test,  plastic  longitudinal  and  shear  waves  propagate  forward  into 
the  anvil  and  backward  into  the  flyer  (Figure  9),  carrying  information  about  the  plates'  plastic 
response.  The  shear  wave  is  especially  telling  of  the  defonnational  behavior,  whereas  the 
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Ceramic  Flyer  Plate 


Figure  8.  The  symmetric  pressure/shear  experiment. 
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F.gure  9.  SP/S  t-x  diagram  showing  elastic  wavefronts  only. 


longitudinal  wave  is  largely  dominated  by  volumetric,  nearly  elastic  behavior.  We  chose  to 
measure  only  the  transverse  velocity  history  in  our  experiments. 

Although  the  laser-Doppler  velocimeter  (LDV)  has  been  in  use  for  nearly  30  years  to 
measure  the  velocity  of  fluid-entrained  particles  and  solid  surfaces,1®  its  use  for  plate  impact 
experiments  was  not  suggested  until  the  TD1  was  introduced.15  To  our  knowledge,  we  are  the 
first  group  to  implement  it.  When  the  TDI  system  is  used,  a  diffraction  grating  must  be  applied 
to  the  anvil  rear  surface  and  the  grating  must  survive  for  the  duration  of  the  experiment.  At  the 
high-impact  velocities  necessary  for  yielding  the  stronger  ceramics,  there  is  doubt  about  whether 
a  polymer  photoresist  grating  can  survive  and  the  expense  of  etching  the  grating  directly  into  the 
ceramic  is  currently  prohibitive.  One  solution  to  this  problem  is  to  confine  the  grating  behind  a 
window  of  sapphire,  as  described  by  Espinosa  and  Clifton.17  Another  solution  is  to  dispense 
with  the  grating  and  use  an  LDV  system. 

A  schematic  of  our  LDV  system,  TSl  Model  9100-3,  is  shown  in  Figure  10.  A  collimated 
beam  of  single-wavelength  (not  necessarily  single-frequency)  laser  light  is  split  by  the  beam¬ 
splitter  into  two  parallel  beams  that  are  made  to  pass  through  an  objective  lens,  which  focuses 
them  to  narrow  waists  and  causes  them  to  cross,  forming  an  ellipsoidal  measurement  volume. 
When  the  collimation  of  the  incoming  beam  is  correct,  the  location  of  the  waists  and  the  center  of 
the  intersection  region  coincide  and  a  group  of  planar,  parallel,  equispaced  interference  fringes  is 
formed  within  the  volume.  The  fringes  are  perpendicular  to  the  plane  of  the  intersecting  beams 
and  parallel  to  a  plane  containing  the  bisector  of  the  angle  between  the  beams. 

There  are  two  ways  to  illustrate  how  an  LDV  works.  Both  ways  lead  to  the  same 
equation  relating  light  modulation  frequency  to  target  transverse  velocity.  An  analogy  to  the  first 
view  is  to  think  of  a  flashlight  moving  behind  a  picket  fence.  Consider  a  light- scattering  object 
moving  through  the  measurement  volume.  Light  will  scatter  strongly  from  the  object  when  it  is 
in  a  bright  fringe  plane  and  weakly  (not  at  all,  ideally)  when  it  is  in  a  dark  fringe  plane.  The 
modulation  frequency  of  the  scattered  light  will  be  the  object  speed  divided  by  the  fringe 
spacing,  which  is  given  by 


d(  = 


A 

2  sin  k 


(7) 
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Moving  Target 


Figure  10.  The  laser  Doppler  velocimeter  (LDV)  system. 


where  X  is  the  wavelength  of  the  laser  light  and  2k*  is  the  beam  crossing  angle.  Thus, 

W\  (8) 

where  d  is  the  beam  spacing  entering  the  objective  and  F  is  the  objective  focal  length.  Thus,  the 
frequency,/,  for  a  given  velocity  is 


/= 


2v  sin  K 
— 


(9) 


The  other,  equivalent  way  to  understand  how  an  LDV  works  is  to  consider  the 
interference  of  light  scattered  from  each  intersecting  beam.  The  light  scattered  from  one  beam 
will  be  Doppler-shifted  fast,  the  other  will  be  slow,  and  die  mixed  light  will  interfere  with  a  beat 
frequency  given  by  Eq.  (9). 

Like  the  TDI,  the  LDV  is  insensitive  to  the  normal  motion  of  the  anvil.  Unlike  the  TD1, 
however,  the  LDV  is  also  insensitive  to  target  tilt  or  to  the  condition  of  the  anvil  rear  surface  as 
long  as  scattered  light  reaches  the  detector.  Thus,  the  LDV  can  be  applied  in  situations  where  the 
TDI  cannot. 

Analysis  Techniques 

The  extraction  of  stresses,  strains,  and  other  mechanical  variables  from  the  transverse 
velocity  history  recorded  in  an  SP/S  test  requires  that  the  test  be  modeled  with  a  postulated 
material  model  and  that  the  model  parameters  be  adjusted  until  the  computed  and  measured 
velocity  histories  match.  Then,  the  desired  mechanical  variables  can  be  extracted  from  the 
simulation.  We  have  used  two  models,  one  a  continuum  model  that  uses  an  elastic-viscoplasdc, 
Mohr-Coulomb,  dilatant  flow  law  and  the  other  a  model  that  views  the  flow  of  ceramic  as  the 
movement  of  blocks  within  a  fragment  bed.  Both  models  are  fully  described  elsewhere4-1®  and 
will  be  only  briefly  outlined  here;  the  description  by  Klopp  and  Shockey18  is  included  in  the 
Appendix. 

The  Mohr-Coulomb  model's  formulation  follows  that  of  Abou-Sayed  and  Clifton  for 
aluminum,19  which  is  in  turn  based  on  the  viscoplastic  formulation  of  Perzyna.20  We  have 
added  comminution,  dilatancy,  and  mean-stress-dependent  yielding  following  the  formulation 
used  by  Rudnicki  and  Rice.21  The  resulting  constitutive  equation  has  the  form 
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(10) 


f'j  “ 


1  +  v 


<*j  -  +  {<k  t, 


where  Cjj  is  the  total  strain  rate,  vis  Poisson's  ratio,  E  is  Young's  modulus,  <5jj  is  the  Kronecker 
delta,  <p  is  the  equivalent  plastic  strain  rate,  t  is  the  effective  shear  stress,  yP  is  the  equivalent 
plastic  shear  strain,  p  =  -Okk/3  is  the  negative  of  the  mean  stress  or  hydrostatic  pressure,  D  is  the 
amount  of  comminution,  and/is  a  flow  potential.  The  bracket  function  <$>  returns  the  value 
of  the  function  <p  when  it  is  positive,  and  zero  otherwise.  The  effective  stress  is  given  in  the  usual 
way  in  terms  of  the  deviatoric  stress,  S,j  =  cry  -  <$jOkk/3:  r  =  (jSjjSjj)1/2.  The  overstress,  ?  -  ry, 
is  proportional  to  the  logarithm  of  the  plastic  strain  rate.  Thus, 


<P 


(11) 


The  yield  stress  ry  depends  on  the  mean  stress,  amount  of  comminution,  and  strain  hardening. 
When  comminution  is  extensive,  D  »  1  and  Mohr-Coulomb  behavior  described  by  ry  ~pp 
dominates,  p  being  the  friction  angle  tangent.  The  amount  of  comminution  evolves  from  0 
toward  1  at  a  rate  given  by  the  amount  the  second  invariant  of  the  total  stress  exceeds  some 
constant  limit. 


.  _  {h  -  Ajm) 

U  2 

P<y0 


(1  -D) 


(12) 


Here  p  is  the  density,  ci  is  the  longitudinal  wave  speed,  and  tQ  is  a  time  constant.  The  choice  of 
12  as  the  driving  force  for  comminution  is  based  on  the  notion  that  even  purely  hydrostatic 
loadings  comminute  polycrystalline  ceramic.  Hydrostatic  loading  causes  comminution  because 
neighboring  grains  often  have  different  crystallographic  orientations  and,  thus,  different  strain 
fields  as  a  result  of  anisotropy.  Cracks  form  when  the  stresses  necessary  to  enforce  the 
compatibility  of  these  strain  fields  across  grain  boundaries  exceed  some  failure  stress.  Recent 
work  by  Lankford3  supplies  strong  microscopic  evidence  of  grain  boundary  cracking  caused  by 
incompatible  plastic  strains  within  neighboring  grains. 

The  constitutive  Eq.  (11),  together  with  the  momentum  and  compatibility  equations,  are 
rewritten  as  equations  along  characteristics.  We  solved  these  equations  by  using  a  modification 
of  the  method  of  Ranganath  and  Clifton.22  The  boundary  conditions  for  the  anvil  plates  are  that 
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the  rear  surface  is  traction-free  and  the  impact  face  has  an  imposed  velocity  of  half  the  impact 
velocity. 

The  fragment  bed  (FRAGBED)  model4  is  based  on  an  analogy  between  fragment  block 
sliding  and  dislocation  theory.  The  basic  relation  is 

djf/dr  =  (r,  -  T)  ( 1 3) 

where  is  the  plastic  strain  on  the  ith  slip  plane;  T/  is  the  shear  stress;  p  is  the  tangent  of  the 
friction  angle;  is  the  normal  stress  on  the  slip  plane;  ps  is  the  solid  density;  C2  is  the  speed  of 
shear  waves;  and  T  =  B/(2gfjC2),  a  time  constant.  B  is  the  block  dimension,  g  is  a  factor  on  the 
order  of  1  relating  the  active  slip  system  to  the  shear  strain  direction,  and.#  is  the  fraction  of 
blocks  that  have  an  adjacent  vacancy  large  enough  to  allow  slip. 

The  block  sliding  description  is  incorporated  into  a  multiple-plane  plasticity  model  23 
originally  developed  for  brittle  fracture  and  shear  band  problems.  In  addition  to  the  block  sliding 
description,  the  model  includes  a  treatment  of  the  creation  and  destruction  of  porosity  via  the 
competing  processes  of  dilatancy  and  compaction.  A  Mie-Griineisen  relation  is  used  for  the 
pressure.  Yield  and  tensile  strengths  are  reduced  by  the  developing  fragmentation.  The 
complete  model  was  incuiporated  into  a  one-dimensional  hydrocode,  called  COPS2,  designed  to 
model  pressure/shear  tests. 

RESULTS 

Pressure/Shear  Tests 

Selected  physical  properties  and  the  source  of  the  tested  ceramics  are  shown  in  Table  1 . 
The  density  and  wavespeed  values  were  provided  by  the  suppliers. 

The  impact  parameters  for  the  three  experiments  on  AIN  are  shown  in  Table  2.  Two 
HSRP/S  and  one  SP/S  experiment  were  performed.  The  histories  of  transverse  particle  velocity, 
shear  stress,  and  strain  rate  from  the  two  HSRP/S  experiments  are  shown  in  Figure  1 1 ,  and  the 
resulting  shear  stress-strain  curves  are  shown  in  Figure  1 2.  The  polymer  grating  on  the  anvil  in 
test  AN-2  deteriorated  so  rapidly  that  only  the  first  4  (of  perhaps  50)  TDI  interference  fringes 
were  captured.  Thus,  only  the  initial  response  is  known.  The  extrapolation  from  the  initial 
response,  assuming  no  strain  hardening,  is  shown  by  a  dotted  line.  The  history  of  transverse 
particle  velocity  from  the  subsequent  SP/S  experiment  is  shown  in  Figure  13. 
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Table  1 


MATERIAL  PROPERTIES  OF  TESTEO  CERAMICS 


P 

Ci 

c2 

Manufacturer’s 

Material 

_ (fl/cmg) _ 

(mm/us) 

(mm/us) 

Manufacturer  Designation 

AIN  3.26  10.8  6.2  Dow  Chemical  RokliteSOO 

(carbothermally 
reduced,  hot- 
pressed) 


85%  AI2O3 

3.42 

8.8 

5.1 

Coors 

AD-85 

99.5%  AI2O3 

3.81 

10.3 

6.2 

Coors 

AD-995 

B4C 

2.51 

13.5 

8.5 

Dow  Chemical 

Roklite  100 
(hot-pressed) 

SiC 

3.10 

11.2 

7.3 

Carborundum 

Sintered 
a  -  SiC 

Table  2 

IMPACT  PARAMETERS  FOR  AIN  EXPERIMENTS 


Nominal 

V0  Specimen  Angle  Pressure 

Experiment  Configuration*  (mm/us) _ Thickness _ (degree) _ (GPa) 


AN-1 

HSRP/s 

0.672 

0.577  mm 

20 

10.9 

AN-2 

HSRP/S 

0.680 

0.599  mm 

20 

11.1 

AN-3 

SP/S 

0.696 

N/Ac 

15 

6.6b 

* HSRP/S  -  high-strain-rate  pressure/shear  (Figures) 

SP/S  -  symmetric  pressure/shear  (Figure  8) 
b  Pressure  variee  significantly  in  SP/S. 
c  There  is  no  specimen  In  an  SP/S  test,  only  Dyer  end  anvil. 


26 


SHEAR  STRAIN  RATE  (ps1)  TRANSVERSE  VELOCITY  fnim/ps} 


TIME  (ns) 

(a)  Transverse  velocity  and  shear  stress  histories 


TIME  (ns) 

(b)  Strain  rate  history 

Figure  1 1 .  Results  from  HSRP/S  tests  on  AIN. 


RAM-6795-41 
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SHEAR  STRESS  (M  Pa) 


Shot  AN-2  (estimate) 


+  +  +  + 


+  +  ++  +  +  +  -»-  +  +  +  ++f+*  +  + 


n - 1 - 1 - 1 - 1 - 1 - r  . 

10  20  30  40  50  60  70x10' 
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Figure  12.  Stress-strain  curves  trom  HSRP/S  tests  on  AIN. 


TRANSVERSE  VELOCITY  (mm/ps) 


figure  13.  Transverse  velocity  history  from  SP/S  experiment  on  AIN. 
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The  impact  parameters  for  the  three  SP/S  experiments  on  AI2O3  are  shown  in  Table  3. 
One  Coors  AD-85  and  two  Coors  AD-995  (Coors  Porcelain  Co.,  Golden,  Colorado)  specimens 
were  tested.  The  resulting  transverse  velocity  histories  are  shown  in  Figure  14.  In  Shots  AO-1 
and  AO-2,  the  impact  faces  were  polished  to  a  «6-nm  RMS  finish  as  measured  by  a  profilometer. 
In  Shot  AO-3,  the  impact  faces  were  roughened  by  sandblasting  to  help  preclude  slip.  The  RMS 
roughness  of  the  sandblasted  surface  was  not  measured,  but  a  scanning-laser  microscope  profile 
of  the  surface  shows  that  sandblasting  created  pits  *=5  pm  in  diameter,  «1  pm  deep,  spaced  »15 
pm  apart  on  the  previously  nearly  pit-free  surface. 

The  impact  parameters  for  the  SP/S  experiment  on  B4C  are  shown  in  Table  4,  along  with 
the  impact  parameters  from  the  two  HSRP/S  experiments  on  AIN  in  which  the  flyer  and  anvil 
plates  sandwiching  the  specimen  were  B4C.  The  reason  for  including  these  results  will  become 
apparent  in  the  Discussion;  the  B4C  in  the  HSRP/S  tests  may  have  failed.  The  resulting 
transverse  velocity  history  from  the  SP/S  experiment  is  shown  in  Figure  15. 

Analyses 

The  stress-strain-pressure  behavior  measured  in  the  two  HSRP/S  experiments  on  AIN 
were  used  to  predict  the  results  of  the  SP/S  experiment  on  AIN.  Shot  AN-1  was  simulated  by 
using  the  Mohr-Coulomb  model,  and  parameters  were  adjusted  to  match  the  measured  response, 
as  shown  in  Figure  1 6.  The  final  parameters  for  AIN  and  the  other  ceramics  tested  are  shown  in 
Table  5.  The  result  of  modeling  the  AIN  SP/S  test  by  using  the  AIN  HSRP/S  parameters  is 
shown  in  Figure  17.  The  prediction  of  the  FRAGBED  model  is  also  shown.  Clearly,  both 
models  significantly  underpredict  the  measured  transverse  velocity.  The  reasons  for  this 
underprediction  will  be  presented  in  the  Discussion. 

The  Mohr-Coulomb  model  parameters  used  to  model  Shot  AO-1  on  AD- 8 5  AI2O3  are 
also  shown  in  Table  5,  and  the  calculated  transverse  velocity  is  shown  as  the  lower  dashed  curve 
in  Figure  18.  Those  same  parameters,  when  used  to  model  Shots  AO-2  and  AO-3  on  AD-995, 
give  the  upper  dashed  curve  in  Figure  18. 

Table  5  also  gives  the  Mohr-Coulomb  model  parameters  for  B4C  that  result  in  a  decaying 
velocity  history,  Figure  19.  The  velocity  history  is  qualitatively  similar  to  that  measured  in  Shot 
BC-1,  Figure  15,  but  the  amplitude  and  decay  rate  are  not  the  same.  We  expect  that  further 
iteration  would  produce  parameters  providing  a  better  match  to  the  experimental  history. 
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Table  3 


IMPACT  PARAMETERS  FOR  SP/S  EXPERIMENTS  ON  AI2O3 

Vo 

Experiment* _ Material _ (mm/ps) _ Impact  Face 


AO-1 

AD-85 

0.421 

Polished 

AO-2 

AD-995 

0.595 

Polished 

AO-3 

AD-995 

0.596 

Roughened 

*  All  at  0  -  15  degrees,  all  symmetric  pressure/shear, 


Table  4 

IMPACT  PARAMETERS  FOR  SP/S  EXPERIMENT  ON  B4C 

Vo 

Experiment* _ (mm/ps) _ Conflquratlonb 


BC-1 

0.939 

SP/S 

BC-2 

0.672 

BSP/S 

BC-3 

0.680 

BSP/S 

0  All  6  -  15  degrees. 

b  SP/S  -  symmetric  pressure/shear 

BSP/S  -  symmetric  pressure/shear  with  AIN  butter. 
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se  vetoc 


TRANSVERSE  VELOCITY  (mm/ps} 


Figure  15.  Transverse  velocity  history  Irom  SP/S  experiment  on  B4C. 
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SHEAR  STRESS  (MPa) 


Figure  16.  Mohr-Coulomb  and  FRAGBED  simulations  of  AIN  HSRP/S  results. 
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Table  5 


SIMULATION  PARAMETERS  FOR  VARIOUS  CERAMICS 


AD-65  & 
AD-995 


AIN 

AbOs 

BaC 

Density  (kg/m3) 

P 

3.26 

3.42 

2.51 

Longitudinal  wave  speed  (mnVps) 

Cl 

10.8 

8.8 

13.7 

Shear  wave  speed  (mm/ps) 

c2 

6.2 

5.1 

8.5 

initial  yield  stress  (GPa) 

"tyo 

1.6 

3.52 

10.0 

Friction  coefficient 

M 

0.057 

0.23 

0.45 

Shear  modulus  (GPa) 

G 

126.0 

88.0 

181.5 

Poisson's  ratio 

V 

0.25 

0.25 

0.19 

Mean  stress  at  HEL  (GPa) 

Phel 

5.0 

3.39 

6.8 

Shear  stress  at  HEL  (GPa) 

*HE1. 

3.0 

2.03 

6.4 

Square  root  of  second  stress  invariant 
at  damage  limit  (GPa) 

luM 

7.0 

4.8 

3.1 

Strain  rate  constant  (s'1) 

io 

0.11 

3.0 

1.0 

Shear  stress  constant  (MPa) 

*0 

1.2 

1.5 

0.5 

Strain  rate  power 

m 

2.0 

2.0 

2.0 

Pressure  constant  (GPa) 

Po 

1.0 

1.0 

1.0 

Time  constant  (ns) 

to 

0.1 

1.0 

0.1 

Hardening  modulus  (GPa) 

h 

1.14 

10.0 

10.0 

Dilatancy  ratio 

13 

0.2 

0.2 

0.2 
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TRANSVERSE  VELOCITY  (mm/ps) 
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Figure  17.  Result  ol  simulating  an  SP/S  experiment  on  AIN  using  parameters  inferred 
from  HSRP/S  experiments. 


36 


TRANSVERSE  VELOCITY  (mrrVps) 


Figure  18.  Results  o'  simulating  SP/S  experiments  on  AI2O3  using  the  Mohr-Coulomb 
model. 

AD-85  strength  parameters  were  used  to  simulate  the  AD-995  experiments. 
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RANSVERSE  VE 


DISCUSSION 

Tabic  6  summarizes  the  Mohr-Coulomb  properties  of  the  ceramics  in  terms  of  the  friction 
coefficient  and  the  shear  strength  at  the  applied  mean  stress,  i.e.,  a  slope  and  intercept  inferred 
from  the  simulations.  To  facilitate  comparison.  Table  6  also  presents  the  properties  in  terms  of 
the  shear  strength  at  a  mean  stress  of  5  GPa.  The  fact  that  the  strengths  regularized  to  5  GPa 
mean  stress  are  all  (but  one)  nearly  the  same  is  remarkable  but  possibly  misleading;  the  strengths 
were  not  measured  at  5  GPa  and  would  possibly  be  different  had  they  been  measured  at  5  GPA 
mean  stress.  Likewise,  the  fact  that  all  but  one  of  the  friction  coefficients  (slopes)  in  Table  6  are 
equal  to  0.23  is  misleading:  except  for  the  AIN  HSRP/S  simulation,  we  did  not  explore  other 
values.  There  is  nothing  to  suggest  that  the  value  0.23  is  unique.  Failure  surfaces  with  high 
slope  and  low  intercept  might  serve  just  as  well  as  failure  surfaces  with  low  slope  and  high 
intercept. 

The  pressure/shear  experiments  have  revealed  a  wide  range  of  mechanical  response 
among  the  ceramics  tested.  The  transverse  velocity  for  AIN  nearly  reached  the  elastic  level, 
whereas  the  velocity  for  B4C  was  only  35%  of  the  elastic  level.  For  AIN  and  B4C,  velocities 
decayed  with  time;  for  AD-85  and  AD-995  AI2O3,  velocities  increased  with  time.  The  former 
behavior  is  characteristic  of  faster  decay  waves  overtaking  the  initial  shear  wave  from  behind, 
whereas  the  latter  behavior  is  characteristic  of  slower  following  waves.  The  different  ceramics 
seem  to  have  different  mechanisms  for  failure  and  postfailure  flow,  and  these  mechanisms 
operate  at  different  threshold  loads. 

Figures  20  and  21  summarize  the  HSRP/S  and  SP/S  results  obtained  for  AIN  and  AD-85 
and  AD-995  AI2O3.  The  data  are  presented  in  terms  of  equivalent  shear  strength  as  a  function  of 
pressure,  deduced  from  the  Mohr-Coulomb  model  parameters,  assuming  105  S'1  shear  strain  rate, 
0.05  equivalent  plastic  strain,  full  damage,  and  a  mean  stress  around  the  average  attained  in 
simulating  a  particular  test.  Figure  20  includes  results  for  AIN  obtained  by  Heard  and  Cline24 
using  quasistatic  compression,  by  Lankford25  using  a  split-Hopkinson  pressure  bar,  and  by 
Grady26  and  Rosenberg  et  al.27  using  normal  plate  impact.  Figure  21  includes  the  results  for 
AD-85  obtained  by  Yeshurun  et  al.,28  Wilkins,29  and  Gust  and  Royce30  and  the  results  for 
AD-995  obtained  by  Gust  and  Royce. 3^  Figure  22  includes  the  results  for  B4C  obtained  by 
Grady,26  Wilkins,29  and  Rosenberg  and  Yeshurun.31  Figure  23,  from  Klopp  and  Shockey,18 
shows  similar  results  obtained  by  Johnson  et  al.,32  by  Grady,26  Kipp  and  Grady,8  and  Klopp  and 
Shockey18  on  silicon  carbide  (SiC). 
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Table  6 


AIN  (HSRP/S) 
AIN  (SP/S) 
AD-85 
AD-995 
SiC 


SUMMARY  OF  RESULTS— INFERRED  MOHR-COULOMB 
PARAMETERS  FOR  VARIOUS  CERAMICS 


Shear  Strength 

Friction  Shear  Strength  Mean  Stress  at  5  GPa  Mean 

Coefficient _ (GPa) _ (GPa) _ Stress 


0.057 

1.8 

0.23 

4.9 

0.23 

1.1 

0.23 

1.8 

0.23 

1.3 

11.0 

1.4 

6.6 

4.5 

3.5 

1.4 

6.5 

1.5 

4.0 

1.5 

40 


SHEAR  FLOW  STRESS  (GPa) 


MEAN  STRESS  (GPa) 


RAM-6795-51 


Figure  20.  Equivalent  shear  How  stress  as  a  Junction  of  mean  stress  Jor  AIN. 
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SHEAR  FLOW  STRESS  (GPa) 


Figure  21 .  Equivalent  shear  flow  stress  as  a  lunction  of  mean  stress  lor  AD-85  and 
AD-995  AI2O3. 
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Figure  22.  Equivalent  shear  flow  stress  as  a  function  of  mean  stress  for  B4C. 
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Figure  23.  Equivalent  shear  How  stress  as  a  function  of  mean  stress  tor  SiC 


From  Reference  18. 
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From  Figures  20  through  23,  wc  conclude  that,  with  the  exception  of  AIN  tested  in  SP/S, 
the  shear  strengths  inferred  from  the  pressure/shear  tests  are  lower  than  the  strengths  reported  by 
others.  As  suggested  by  Klopp  and  S  hockey, 18  we  think  the  reason  is  that  we  have  essentially 
measured  the  response  of  comminuted  ceramic.  The  longitudinal  wave,  which  is  strong  enough 
to  cauy  damage  even  in  compression,  runs  ahead  of  the  shear  wave.  The  shear  wave,  which  we 
measured,  is  always  running  through  damaged  material.  Thus,  the  mean-stress-strength  behavior 
of  SiC,  for  example,  measured  in  pressure/shear  is  closer  to  lying  along  a  line  through  the  points 
of  Johnson  et  al.28  for  ball-milled  SiC  than  along  the  line  connecting  data  points  for  intact 
ceramic  (Figure  23). 

We  established  at  the  outset  that  any  tests  for  measuring  the  ballistic  behavior  of  armor 
ceramics  must  measure  their  behavior  under  penetration  conditions  (at  least  until  such  tests  show 
that  ballistic  behavior  is  independent  of  certain  penetration  conditions).  Both  the  HSRP/S  and 
SP/S  tests  have  met  this  requirement,  because  mean  stresses  and  strain  rates  characteristic  of  the 
Mescall  zone  were  achieved  in  both.  For  example,  the  nominal  mean  stress  and  strain  rate  in  the 
HSRP/S  tests  on  AIN  were  1 1  GPa  and  1.3  x  105  s:,  respectively,  and  Figure  24  shows  that  the 
mean  stress  and  plastic  shear  strain  rate  were  4  GPa  and  1.0  x  104  s'1,  respectively,  obtained 
from  the  Mohr-Coulomb  simulation  of  the  SP/S  test  on  AD-85. 

The  Mohr-Coulomb  simulations  can  reproduce  the  transverse  velocity  histories  at  least 
qualitatively  and,  in  the  cases  of  AIN  tested  in  HSRP/S  and  AD-85  and  AD-995  AI2O3  tested  in 
SP/S,  quantitatively.  Matching  the  AIN  and  AD-85  velocity  histories  was  essentially  a  curve- 
fitting  exercise.  However,  the  inability  of  the  simulation  to  predict  the  SP/S  response  of  AIN 
from  the  HSRP/S  results  and  its  ability  to  predict  the  SP/S  response  of  AD-995  from  the  SP/S 
response  of  AD-85  are  significant. 

The  response  of  AIN  inferred  from  the  HSRP/S  tests  is  soft.  The  response  inferred  from 
the  SP/S  test  is  hard.  Since  it  is  unlikely  that  the  SP/5  test  could  indicate  a  response  that  is  too 
hard  (too  high  a  transverse  velocity),  the  HSRP/S  results  should  be  questioned.  There  are  several 
reasons  why  an  HSRP/S  test  might  indicate  too  low  a  transverse  velocity  and,  hence,  too  soft  a 
response:  (1)  the  impact  interface  may  have  slipped  and  transmitted  too  little  traction  into  the 
anvil,  (2)  the  B4C  flyer  and  anvil  may  have  failed  in  compression  by  fracturing  or  yielding  in  the 
interior,  and  (3)  the  rear  surface  of  the  anvil  under  the  TDI  diffraction  grating  or  the  grating  itself 
may  have  failed  in  shear.  The  evidence  supports  the  third  possibility  more  than  the  others. 
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Tlie  evidence  against  slip  is  that  both  polished  and  sandblasted  AD-995  gave  the  same 
response  (Figure  14).  We  would  expect  that,  if  slip  were  occurring,  changing  the  frictional 
characteristics  of  the  impact  interface  would  make  a  difference.  The  frictional  characteristics  of 
AIN  against  B4C  may  be  different  from  those  of  AD-995  against  itself.  The  fact  that  two 
HSRP/S  tests  on  AIN  gave  nearly  the  same  response  further  suggests,  however,  that  slip  was  not 
occurring.  Slip  is  probably  highly  sensitive  to  cleanliness,  roughness,  gun  chamber  vacuum,  and 
projectile  blow-by,  and  these  conditions  were  not  identical  in  the  two  tests. 

The  second  possibility,  that  the  B4C  failed  in  compression  in  the  interior,  would  require 
that  it  had  been  loaded  beyond  the  failure  limit  in  compression/shear.  The  failure  limit  of  B4C  in 
uniaxial  strain  compression,  the  HEL,  is  reportedly  greater  than  14  GPa.29  The  maximum 
equivalent  shear  stress  that  could  have  been  applied  to  the  B.tC,  assuming  the  AIN  specimen  did 
not  yield,  was  around  5.4  GPa,  whereas  the  failure  limit,  assumed  to  be  the  equivalent  shear 
stress  at  the  HEL,  was  around  6.2  GPa.  Thus,  the  B4C  probably  did  not  fail  in  compression  in 
the  interior. 

The  third  possibility,  that  the  B4C  failed  near  the  diffraction  grating,  seems  the  most 
likely.  The  B4C  probably  failed  in  tension  on  45-degree  planes  when  the  initial  shear  wave  ran 
into  the  unloaded  region  near  the  back  surface,17  as  indicated  in  Figure  9.  In  all  the  experiments, 
the  45-degree  tensile  stresses  far  exceeded  the  preshocked  spall  strengths  reported  by  Grady.26 
In  addition,  the  FRAGBED  model  includes  a  simple  tension  cutoff  failure  model.  In  a 
simulation  of  AIN  loaded  in  SP/S,  the  model  predicts  tensile  failure  near  the  location  indicated  in 
Figure  9  and  the  anvil  rear  surface  ends  up  moving  at  a  greater  normal  velocity  than  the  normal 
component  of  the  projectile  velocity.  Nevertheless,  the  results  show  that  all  the  ceramics 
maintained  at  least  some  shear  strength  after  high-amplitude  shock  and  release. 

The  fact  that  the  same  model  parameters  can  match  both  die  AD-85  and  AD-995 
responses  does  not  mean  they  have  the  same  mechanical  properties.  AD-85  is  well  known  to  be 
weaker  than  AD-995.  If  the  AD-85  had  been  tested  at  the  same  impact  velocity  (595  instead  of 
421  m/s)  as  the  AD-995,  presumably  the  AD-85  rear-surface  velocity  would  have  been  lower. 
The  match  between  the  results  shows  the  fallacy  of  using  the  models  to  extrapolate  away  from 
the  conditions  of  a  given  test.  The  results  are  valid  at  least — perhaps  only — for  the  conditions  of 
a  given  test.  Extrapolation  using  the  current  model  is  not  warranted. 
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Chapter  4 

OTHER  TECHNIQUES 


ROD  IMPACT 

The  rod  impact  technique,  or  Taylor  test,  has  been  applied  to  metals,33  polymers,34  and 
ceramics.35  In  the  test,  a  circular  rod  about  4  to  6  diameters  long  is  impacted  axially  against  a 
"rigid"  wall  or  another  like  rod.  The  deformed  profile  of  the  rod  is  recorded  with  a  framing 
camera33  or  determined  posttest,  and  this  profile  is  compared  with  that  predicted  from  analysis. 
The  material  model  in  the  analysis  is  adjusted  until  the  analytical  and  experimental  profiles 
match.  The  main  advantage  of  the  rod  impact  test  over  plate  impact  tests  is  that  less 
sophisticated  experimental  facilities  are  required.  The  main  disadvantage  is  that  the  loading  on 
the  rod  is  unknown  and  complicated  and  must  be  inferred  from  computations,  making  the 
quantitative  determination  of  material  properties  difficult. 

Another  disadvantage  is  that  the  stress  state  in  the  rod  impact  test  is  nearly  uniaxial,  and 
therefore  the  mean  stress  can  never  significantly  exceed  the  deviatoric  stress,  unlike  the  stress 
state  ahead  of  a  penetrator.  Thus,  to  be  relevant  to  armor  penetration  measurements,  the  test 
must  be  modified  to  increase  the  mean  stress.  The  modification  we  used  was  to  contain  the  rod 
within  a  cylinder  of  another  material,  thereby  constraining  radial  expansion  of  the  rod. 

We  simulated  a  contained  rod  impact  test  with  a  hydrocode  and  deduced  that  the  stress 
and  strain  paths  were  similar  to  those  occurring  ahead  of  a  penetrator.  We  then  performed 
several  rod  impact  experiments  in  which  a  ceramic  rod  was  surrounded  by  a  metal  or  PMMA  and 
metal  containment.  Two  examples  are  shown  in  Figure  25.  The  ends  of  the  assemblies  were 
impacted  by  explosively  accelerated  plates.  Because  the  ceramic  rod  was  unlikely  to  be 
recoverable  for  examination  posttest,  we  hoped  that  the  profile  would  be  recorded  in  the  bore  of 
the  containment. 

We  tried  the  first  configuration  shown  in  Figure  25  (a)  and  found  that,  indeed,  the  rod 
could  not  be  recovered;  the  containment  also  was  not  recoverable.  In  an  attempt  to  capture  and 
preserve  the  ceramic,  we  designed  a  blunt  tapered  flyer  that  was  supposed  to  explosively  weld 
itself  to  the  containment  cylinder  (Figure  25b).  However,  either  the  flyer  velocity  was  too  low 
relative  to  its  strength  or  the  taper  angle  was  wrong;  in  any  case,  welding  did  not  occur. 
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Figure  25.  Rod  'mpact  configurations. 

(All  dimensions  in  centimeters) 
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Although  some  contained-rod  impact  configuration  could  be  made  to  work,  we  did  not 
pursue  this  problem  any  further. 

CAVITY  EXPANSION 

The  success  of  cavity  expansion  analyses^ 38  of  modeling  penetration  suggests  that 
cavity  expansion  tests  would  exercise  a  material  under  penetration  conditions.  We  have  analyzed 
the  expansion  of  a  Mohr-Coulomb  (presumably  ceramic)  cylinder  confined  within  an  elastic/ 
plastic  (presumably  metal)  cylinder  (Figure  26),  subjected  to  monotonically  increasing  internal 
pressurization.  Our  analysis  is  based  on  Hill's39  analysis  of  the  expansion  of  a  monolithic  tube 
and  is  valid  from  the  initial  yielding  of  the  ceramic  to  complete  yielding  of  the  metal 
confinement.  We  provide  relations  for  the  radius  of  the  elastic-plastic  boundary,  when  in  either 
the  ceramic  or  the  metal  confinement,  and  the  outer  radius  of  the  confinement  as  a  function  of  the 
internal  pressure.  The  loading  is  assumed  to  be  quasi-static,  plane  strain,  and  the  ceramic  and 
metal  are  modeled  as  elastic-perfectly-plastic.  The  Lam6  elastic  constants  for  the  ceramic  and 
metal  are,  respectively,  A,  fi  and  X$,  11$.  Postyield,  the  ceramic  is  assumed  to  obey  a  Mohr- 
Coulomb  flow  law  with  a  tension  cutoff  and  the  metal  is  assumed  to  obey  Tresca's  flow  law. 

In  Figure  26,  a  and  b  are  the  inside  and  outside  radii,  respectively,  of  the  ceramic 
cylinder,  and  b  and  d  are  the  inside  and  outside  radii,  respectively,  of  the  metal  confinement. 

The  radius  of  the  elastic-plastic  boundary  is  c.  The  yield  condition  of  the  ceramic  is  shown 
graphically  in  Figure  27.  For  monotonic  pressurization,  stress  states  are  on  the  branch 


(14) 


<*00  =  mOrr  +  Yj 


where  m  is  the  slope  defined  in  Figure  27,  Yj  is  the  tension  cutoff,  and  and  c&e  are  the  radial 
and  tangential  stress  components,  respectively.  The  yield  condition  is  equivalent  to 


(15) 


xZ/P  +  T, 


where  r  is  the  maximum  in-plane  shear  stress,  P  is  the  negative  of  the  in-plane  mean  stress, /is 
the  tangent  of  the  friction  angle,  and  T  is  the  cohesion  in  shear. 


(16) 
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and 


*T  = 


2  T 

1  +/ 


(17) 


The  boundary  conditions  are  that  Oh  =  -Pq  at  r  =  a  and  Oh  =  0  at  r  =  d.  Further  conditions 
necessary  for  solving  the  problem  are  that  the  material  just  inside  r  =  c  is  at  incipient  yield  and 
that  traction-  are  continuous  across  r-c.  When  the  elastic-plastic  boundary  is  in  the  ceramic, 
the  values  of  c  are  the  roots  of 

~(Pq  +  i“^)(l  -  ~  (Po  +  +  m)a]-mcm-  ]d2  +  YTdP^~  =  0  (18) 


where  &■  =  fi/2  and  d  is  the  outside  radius  of  a  ceramic  containment  with  elastic  constants  A,  fi 
equivalent  to  a  metal  containment  of  outside  radius  d  with  elastic  constants  As,  fis  and 


r  _  { [(A  +  /i)  -  +  ps )]/ish2  +  (A  +  £  +  /isXAs  +  /isM2  }/ih2 

t^s(As  +  /is  +  p)tp-  +  (M  -  /is)(As  +  /is)^](A  +  /iV2 


(19) 


The  pressure  for  initial  yielding  of  the  ceramic  (i.e.,  when  c  -  a)  is  given  by 


Pva  — 


[(m  -  1  )a2  ~(m  + 
and  the  pressure  for  final  yielding  of  the  ceramic  (c  -  bm)  is  given  by 


^yb’  = 


^ j^[(l  +mjd 2  +  (1  - m)bt\  -  2^YT 

-(1  t(l  +  mV/2  +  (1  -  m)b2]  | 


(20) 


(21) 


The  relation  between  the  internal  pressure  and  the  change  in  the  outside  radius  of  the 
confinement,  through  final  yielding  of  the  ceramic,  is  given  by 

(A  +  2/i)(As+  2fis)YTlP-d 
u~  D 


(22) 
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where 


D  =  2(1  -m)(A  +  /i)[/4(As  +  Ms  +  M)^2  +  (M“  Ms)(As  +  Ms )^]  + 
2(1  +  m)^  {  f (A  +  //)  -  (As  +  Ms^Ms*4  +  (A  +  m  +  MsXAs  +Ms)^  d2  } 


(23) 


When  the  elastic-plastic  boundary  is  in  the  metal,  its  radius  is  given  implicitly  by 

1/2 


c  - 


(24) 


The  pressure  for  initial  yielding  of  the  confinement  (c  =  b+ )  is  given  by 

Yx 

m 


sr-i5 


(25) 


and  the  pressure  for  complete  yielding  of  the  confinement  (c  =  d)  is  given  by 

xWniw.xI  -  m  yT 


pyi^[ri 


m 


+  nu,(5H5 


1  -m 


(26) 


Finally,  the  relation  between  internal  pressure  and  the  change  in  the  outside  radius  of  the 
confinement,  through  final  yielding  of  the  confinement,  is  given  by 


u(d)  = 


hmf)+ ^-('o + + ^ 


2(A  +  m)M(c2  ~  d2) 


(27) 


Figure  28  shows  the  plastic  rone  radius  and  the  confinement  radial  displacement  as  a 
function  of  the  internal  pressure  for  an  example  set  of  parameters,  given  in  Table  7.  Results  are 
plotted  for /  =  0.25, 0.5,  and  0.75.  The  plots  show  that,  provided  the  metal  confinement  is  well 
characterized,  the  frictional  behavior  of  the  ceramic  can  be  inferred  from  measuring  the  change 
in  the  radius  of  the  confinement  or  the  plastic  zone  radius.  Large  relative  changes  in / lead  to 
large  relative  changes  in  u(d)  for  a  given  pressure. 


CONFINEMENT  RADIAL  DISPLACEMENT  u(d)  (pm)  PLASTIC  ZONE  RADIUS  c  (mm) 


0.0  0.5  1.0  1.5  2.0  2.5 


INTERNAL  PRESSURE  P0  (GPa) 

(b)  Conlinement  radial  displacement 

RAM-6795  59 

Figure  28.  Plastic  zone  radius  and  confinement  radial  displacement  versus 
internal  pressure. 
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Table  7 


PARAMETERS  FOR  CYLINDRICAL  CAVITY 
EXPANSION  EXAMPLE  PROBLEM 


Parameter  Value 


a  (mm) 

10 

b  (mm) 

50 

d(mm) 

70 

X  (GPa) 

115 

(GPa) 

147 

Xs  (GPa) 

119 

Ms  (GPa) 

80 

Yj  (MPa) 

500 

Ys  (MPa) 

350 
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The  above  analysis  is  strictly  quasi-static.  We  would  expect,  however,  that  a  dynamic 
analysis  taking  wave  propagation  into  account  would  also  show  that  u(d,t )  is  sensitive  to  the 
internal  pressure  history.  The  dynamic  analysis  would  apply  when  the  internal  pressure  is 
supplied,  for  example,  by  the  detonation  of  an  explosive. 

PENETRATION  ARREST  EXPERIMENTS 

To  understand  ceramic  penetration,  we  must  also  understand  penetrator  erosion.  Much  of 
the  phenomenology  of  ceramic  penetration  can  be  obtained  by  sectioning  ballistically  tested 
targets.  However,  penetrators  are  often  completely  eroded  or  otherwise  badly  deformed  after 
steady-state  penetration  and  little  information  is  available  on  erosion  phenomenology.  To  better 
understand  erosion,  we  designed  1/4-scale  penetration  experiments  in  which  penetration  would 
be  arrested  partway  through  the  targets. 

We  shot  6-ram -diameter,  60-mm-long  tungsten-nickel-iron  (W-Ni-Fe)  alloy  rods  with 
hemisphencal  tips  into  highly  confined  Coors  AD-995  alumina  ceramic  targets,  as  shown  in 
Figure  29.  Targets  were  constructed  by  shrink-fitting  steel  rings  around  25-mm-thick  ceramic 
disks,  stacking  two  such  assemblies,  and  attaching  front  and  back  plates  with  bolts.  The  ceramic 
disks  were  made  1 3  pm  thicker  than  the  rings  and  100  pm  larger  in  diameter  than  the  inside 
diameter  of  the  rings,  so  that  the  disks  were  under  an  axial  prestress  of  around  20  MPa  and  a 
radial  prestress  of  around  500  MPa  (both  prestresses  compressive).  The  steel  parts  were  heat- 
treated  to  HRC  40  before  assembly.  Impact  velocities  were  around  850  m/s. 

Penetration  was  arrested  partway  through  the  targets  as  desired.  The  tail  of  the  penetrator 
was  visible  through  the  hole  in  the  cover  plate.  The  ceramic  was  significantly  dilated  as 
evidenced  by  the  bulging  of  the  front  plate,  opening  of  gaps  between  the  rings,  and  stretching  of 
the  bolts.  Unfortunately,  we  did  not  have  sufficient  time  or  funds  to  section  the  targets  and 
examine  the  penetrators.  We  recommend  that  sectioning  be  pursued  in  a  future  program.  We 
have  available  several  more  AD-995  targets  and  targets  of  TiB2.  We  also  have  six  uranium- 
3/4%  titanium  (DU-3/4Ti)  penetrators.  By  performing  further  experiments,  we  can  compare  the 
erosive  properties  of  AD-995  and  TiB2  and  the  erosion  resistance  of  DU-3/4  Ti  and  W-Ni-Fe 
alloys. 
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Figure  29.  Apparatus  for  penetration  arrest  experiments. 
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Chapter  S 

CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

From  the  results  of  the  research  performed  in  this  program,  we  reached  the  following 
conclusions: 

•  Penetration  of  a  confined  ceramic  block  by  an  impacting  long  rod  projectile 
occurs  by  the  comminution  of  a  zone  of  the  ceramic  at  the  leading  edge  of  the 
projectile  and  the  flow  of  the  fine  fragments  out  of  the  projectile  path. 

•  The  material  properties  governing  the  penetration  resistance  of  a  confined 
ceramic  block  include  the  dynamic  compressive  strength,  the  comminution  and 
granular  flow  resistance  under  high  rate  and  high  pressure  conditions,  and  the 
abrasiveness  of  the  ceramic  fragments. 

•  Pressurc/shear  tests  in  two  configurations  with  rear  surface  laser  diagnostics 
and  with  appropriate  constitutive  analysis  are  suitable  for  measuring  the  high 
pressure,  high  strain  rate  comminution  and  granular  flow  properties  of  armor 
ceramics. 


recommendations  for  future  work 

Our  research  on  the  phenomenology  of  penetration,  the  development  of  tests  for 
measuring  governing  properties,  and  the  application  of  those  tests  to  candidate  armor  ceramics 
suggest  many  avenues  for  future  research,  in  terms  of  both  specific  test  techniques  and  armor 
design  in  general. 

One  issue  that  needs  to  be  addressed  is  whether  the  pressure/shear  anvils  fail  in  45-degree 
tension  when  the  shear  wave  runs  into  uncompressed  material.  One  way  to  remove  doubt  in 
future  experiments  is  to  simply  prevent  shear  failure.  A  sapphire  window  placed  behind  the 
anvil,  as  implemented  by  Espinosa  and  Clifton,!7  delays  longitudinal  unloading  and  ensures  that 
the  shear  wave  always  runs  into  compressed  specimen  material.  Another  way  is  to  use  in- 
material  gages  to  make  measurements  nearer  the  impact  plane  than  the  potential  shear  failure 
plane.  For  nonconductive  ceramics,  the  magnetic  particle  velocity  technique  developed  by  Gupta 
et  al.4*!  is  a  proven  method  for  making  such  measurements.  Both  types  of  experiments  are 
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recommended,  because  the  latter  technique  cannot  be  applied  to  conductive  ceramics  such  as 
TiB2,  whereas  the  former  technique  is  restricted  to  cases  in  which  the  sapphire  does  not  yield. 

A  second  issue  that  should  be  addressed  is  whether  slip  occurs  between  the  impacting 
plates.  The  magnetic  particle  velocity  technique  can  also  be  applied  to  resolve  this  question.  If 
slip  is  not  occurring  in  a  symmetric  pressure/shear  experiment,  the  particle  velocity  detected  by  a 
gage  located  on  the  impact  face  of  the  anvil  will  jump  to  half  the  impact  velocity  and  remain 
there  for  the  duration  of  the  one-dimensional  part  of  the  experiment  Any  variation  in  velocity, 
in  an  otherwise  perfectly  symmetric,  planar  impact,  indicates  that  slip  is  occurring. 

To  reduce  ambiguity  in  the  measured  stress-strain  response  of  very  hand  ceramics  like 
TiB2,  SiC,  and  B4C,  we  recommend  HSR  pressure/shear  experiments  using  flyer  and  anvil  plates 
of  cubic  boron  nitride  (CBN).  Since  CBN  has  a  higher  Hugoniot  elastic  limit  than  these 
candidate  armor  ceramics,  the  flyer  and  anvil  plates  will  remain  elastic  at  the  high  impact 
velocities  required  to  fail  the  specimen  materials,  allowing  a  direct  determination  of  the  stress- 
strain  response.  CBN  is  available  relatively  cheaply  in  disks  25  mm  in  diameter  and  2  mm  thick 
as  a  cutting  tool  material. 

The  fragmentation  and  flow  models  should  be  further  developed.  The  Mohr-Coulomb 
model  we  used  is  largely  empirical;  how  the  parameters  in  the  model  relate  to  physical  properties 
and  processes  in  the  material  is  unclear.  The  block-sliding  model,  FRAGBED,  is  more 
physically  based,  but  is  incomplete.  We  recommend  that  FRAGBED  be  modified  to  include 
details  of  the  fracture  and  fragmentation  process.  The  result  will  be  improved  understanding  of 
fragment  size  and  shape  evolution,  fragment  interaction,  and  the  possible  interplay  of  plastic 
(dislocation-based)  deformation  and  fragmentation.  The  relationship  between  model  parameters 
and  physics  would  be  more  clear.  The  model  would  be  used  for  simulating  both  SP/S 
experiments  and  ballistic  penetration. 

Our  data  base  of  pressure/shear  results  is  limited  in  the  number  of  ceramics  tested  and  the 
range  of  mean  stress  and  strain  rates  covered.  We  recommend  that  additional  pressure/shear  tests 
be  performed,  particularly  on  ceramics  of  common  interest  to  other  investigators,  such  as  AD- 
995,  and  to  designers,  such  as  TiB2. 

An  investigation  into  penetrator  erosion  is  recommended.  The  following  questions 
needing  to  be  addressed:  What  ceramic  properties  or  ceramic  fragment  characteristics  control 
penetrator  erosion  rates?  How  do  abrasive  wear  and  plastic  flow  relate  to  erosion?  If  sharp, 
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angular  fragments  erode  pcnetrators  better,  what  properties  control  fragment  morphology,  and 
can  these  properties  be  enhanced  without  degrading  other  important  properties? 

We  still  need  a  simple,  inexpensive  test  that  exercises  ceramics  under  penetration 
conditions,  because  the  pressure/shear  tests  described  here  require  a  sophisticated  gas  gun  facility 
with  accompanying  interferometers  and  high-speed  recording  equipment.  We  recommend  that 
the  search  for  a  simpler  test  begin  with  continued  investigation  of  the  cylindrical  cavity 
expansion  test.  That  test  requires  an  explosives  test  capability  and  the  capability  to  machine 
cylinders  of  armor  ceramics,  but  no  more  sophisticated  instrumentation  than  a  micrometer  is 
required  :o  measure  the  diameter  of  the  confinement  before  and  after  the  test.  Our  results  show 
that  the  granulation  and  flow  properties  of  armor  ceramics  are  relatively  insensitive  (compared 
with  metals,  say)  to  loading  rate,  but  are  sensitive  to  mean  stress.  Thus,  a  simpler  static  test  may 
be  possible,  but  only  if  sufficiently  large  mean  stresses  can  be  applieJ.  A  dynamic  test  may  be 
necessary,  if  only  to  attain  sufficiently  large  mean  stress. 


63 


REFERENCES 


1.  D.  A.  Shockey,  A.  H.  Marchand,  J.  E.  Crocker,  and  R.  W.  Klopp,  “Dynamic  Failure 
Behavior  of  Ceramic  Targets  and  Textured  Tungsten  Alloy  Rods,"  Final  Technical  Report 
on  Contract  No.  9-X69-3295Z-1,  Los  Alamos  National  Laboratory,  SRI  Project  7525 
(November  1991). 

2.  D.  A.  Shockey,  A.  H.  Marchand,  S.  R.  Skaggs,  G.  E.  Cort,  M.  W.  Burkett,  and  R.  Parker, 
“Failure  Phenomenology  of  Confined  Ceramic  Targets  and  Impacting  Rods,” /rtf.  J.  Impact 
Eng.  9(3),  1  (1989). 

3.  J.  Lankford,  “The  Compressive  Strength  of  Strong  Ceramics:  Microplasticity  versus 
Microfracture,”  J.  Hard  Mater.  2, 55-77  (1991). 

4.  D.  R.  Curran,  L.  Seaman,  T.  Cooper,  and  D.  A.  Shockey,  “Micromechanical  Model  for 
Comminution  and  Granular  Flow  of  Brittle  Material  under  High  Pressure  and  High  Strain 
Rate:  Application  to  Penetration  of  Ceramic  Armor,”  submitted  to  Int.  J.  Impact  Eng. 
(1992). 

5.  J.  Mescall  and  C.  Tracy,  “Improved  Modeling  of  Fracture  in  Ceramic  Armor,”  Proceedings 
of  the  1986  Army  Science  Conference,  U.S.  Military  Academy,  West  Point,  17-20  June 
1986. 

6.  J.  Mescall  and  W.  Weiss,  “Materials  Behavior  under  High  Stress  and  Ultrahigh  Loading 
Rates-Part  II,”  Proceedings  of  the  29th  Sagamore  Army  Conference,  Army  Materials  and 
Mechanics  Research  Center,  Watertown,  MA  (1984). 

7.  T.  Cooper,  “A  Computer  Code  for  Numerical  Simulation  of  Shock  Waves  in  Fluids  and 
Solids,”  Report  DS  1980: 16,  Swedish  Detonic  Research  Foundation,  Stockholm,  Sweden 
(1980). 

8.  M.  E.  Kipp  and  D.  E.  Grady,  "Shock  Compression  and  Release  in  High-Strength 
Ceramics,"  Sandia  Report  SAND89-146l*UC-704  (July  1989). 

9.  D.  C.  Erlich,  “Adaptations  to  Existing  Gas  Guns  to  Increase  Impact  Velocity  and  Flyer 
Plate  Diameter,”  Proceedings  of  the  42nd  Aeroballisric  Range  Association  Meeting,  21-25 
October  1991,  Adelaide,  South  Australia. 

10.  R.  W.  Klopp,  "Fiber-Optic  Interferometers,"  Poulter  Laboratory  Technical  Report  001-89, 
SRI  International,  Menlo  Park,  California  (1989). 

11.  P.  Kumar  and  R.  J.  Clifton,  “Optical  Alignment  of  Impact  Faces  for  Plate  Impact 
Experiments,”  J.  Appl.  Phys.  48,  1366-1367  (1977). 

12.  A.  S.  Abou-Sayed,  R.  J.  Clifton,  and  L.  Hermann,  “The  Oblique  Plate  Impact  Experiment,” 
Exp.  Mech.  16, 127-132  (1976). 


65 


13.  C.  H.  Li,  "A  Pressure-Shear  Experiment  for  Studying  the  Dynamic  Plastic  Response  of 
Metals  at  Shear  Strain  Rates  of  105  s'1,"  Ph.D.  Thesis,  Brown  University,  Providence, 
Rhode  Island  (1981). 

14.  R.  J.  Clifton  and  R.  W.  Klopp,  "Pressure-Shear  Plate  Impact  Testing,"  Metals  Handbook, 
Vol.  8.,  9th  ed.  (ASM  International,  Metals  Park,  Ohio  1985),  pp.  230-239. 

15.  K.  S.  Kim,  R.  J.  Clifton,  and  P.  Kumar,  "A  Combined  Normal-  and  Transverse- 
Displacement  Interferometer  with  an  Application  to  Impact  of  y-Cut  Quartz,"  J.  Appl.  Phys. 
48(10),  4132-4139  (1977). 

16.  M.  J.  Rudd,  “A  New  Theoretical  Model  for  the  Laser  Dopplermeter,”  J.  Phys.  E.  2(2),  55- 
58(1969). 

17.  H.  D.  Espinosa  and  R.  J.  Clifton,  “Plate  Impact  Experiments  for  Investigating  Inelastic 
Deformation  and  Damage  in  Advanced  Materials,”  in  Experiments  in  Micromechanics  of 
Failure-Resistant  Materials,  K.  S.  Kim,  Ed.  (ASME,  New  York,  1991). 

18.  R.  W.  Klopp  and  D.  A.  Shockey,  “The  Strength  Behavior  of  Granulated  Silicon  Carbide  at 
High  Strain  Rates  and  Confining  Pressure,”  J.  Appl.  Phys.  70(12),  7318-7326  (1991). 

19.  A.  S.  Abou-Sayed  and  R.  J.  Clifton,  “Analysis  of  Combined  Pressure-Shear  Waves  in  an 
Elastic/Viscoplastic  Material,”  J.  Appl.  Mech.  44,  79-84  (1977). 

20.  P.  Perzyna,  “Fundamental  Problems  in  Viscoplasticity,”  Adv.  Appl.  Mech.  9, 243-377 
(1966). 

21.  J.  W.  Rudnicki  and  .T.  R.  Rice,  “Conditions  for  the  Localization  of  Deformation  in  Pressure- 
Sensitive  Dilatant  Materials,”  J.  Mech.  Phys.  Solids  23, 371-394  (1975). 

22.  S.  Ranganath  and  R.  J.  Clifton,  “A  Second-Order  Accurate  Difference  Method  for  Systems 
of  Hyperbolic  Partial  Differential  Equations,”  Comput.  Methods  Appl.  Mech.  Eng.  1,  173- 
188(1972). 

23.  L.  Seaman,  “Development  of  a  Model  for  Shear  Banding:  SHEAR  3,”  Final  Report,  Vol. 

2,  Contract  DAAKl  1-78-C-0115  to  Ballistics  Research  Laboratory,  SRI  Project  No.  PYU- 
7893  (April  1983). 

24.  H.  C.  Heard  and  C.  F.  Cline,  "Mechanical  Behavior  of  Polycrystalline  BeO,  AI2O3,  and 
AIN  at  High  Pressure,"  J.  Mater.  Sci.  15(8),  1889-1897  (1980). 

25.  J.  Lankford,  personal  communication  (November  1989). 

26.  D.  E.  Grady,  “Dynamic  Material  Properties  of  Armor  Ceramics,”  Sandia  Report  SAND91- 
0147*UC-704  (March  199 1 ). 

27.  Z.  Rosenberg,  N.  S.  Brar,  and  S.  J.  Bless,  “Dynamic  High  Pressure  Properties  of  AIN 
Ceramic  as  Determined  by  Flyer  Plate  Impact,"  submitted  to  J.  Appl.  Phys.  (1992). 


66 


28.  Y.  Yeshurun,  Z.  Rosenberg,  D.  G.  Brandon,  and  N.  A.  Travitzky,  “Effect  of  Heat 
Treatment  on  the  Dynamic  Response  of  Commercial  85  wt%  AI2O3,”  Mat.  Sci.  and  Engng. 
71,71-75(1985). 

29.  M.  L.  Wilkins,  “Use  of  Boron  Compounds  in  Lightweight  Armor,”  in  Boron  and 
Refractory  Borides,  V.  I.  Matkovich,  Ed.  (Springer- Verlag,  Berlin,  1977). 

30.  W.  H.  Gust  and  E.  B.  Royce,  “Dynamic  Yield  Strengths  of  B4C,  BeO,  and  AI2O3 
Ceramics,”  J.  Appl.  Phys.  42(1),  276-295  (1971). 

31.  Z.  Rosenberg  and  Y.  Yeshurun,  “The  Relation  Between  Ballistic  Efficiency  and 
Compressive  Strength  of  Ceramic  Tiles,”  Int.  J.  Impact  Engng.  7(3),  357-362  (1988). 

32.  G.  R.  Johnson,  T.  J.  Holmquist,  J.  Lankford,  C.  E.  Anderson,  and  J.  Walker,  “A 
Computational  Constitutive  Model  and  Test  Data  for  Ceramics  Subjected  to  Large  Strains, 
High  Strain  Rates,  and  High  Pressures,”  Final  Technical  Report  for  Optional  Task  1, 
Department  of  Energy  Contract  DE-AC04-87AL-42550,  Advanced  Technology 
Assessment  Center,  Los  Alamos  National  Laboratory  (August  1990). 

33.  D.  C.  Erlich  and  P.  Chartagnac,  “Determination  of  Dynamic  Flow  Curve  of  Metals  at 
Ambient  and  Elevated  Temperatures  by  Rod  Impact  Techniques,”  J.  de  Physique,  Colloque 
C5,  supplement  au  no.  8,  Tome  46,  C5-455-462  (1985). 

34.  D.  C  Erlich  and  L.  Seaman,  “Dynamic  Failure  of  Energetic  Propellants,”  SRI  International 
unpublished  results  (1990). 

35.  N.  S.  Brar  and  S.  J.  Bless,  in  Proceedings  of  the  International  Conference  on  Shock-Wave 
and  High-Strain-Rate  Phenomena  in  Materials  (EXPLOMET  90),  M.  A.  Meyers,  L.  E. 
Mutt,  and  K.  P.  Staudhammer  Eds.  (Dekker,  New  York,  1991). 

36.  R.  F.  Bishop,  R.  Hill,  and  N.  F.  Mon,  “The  Theory  of  Indentation  and  Hardness  Tests,” 
Proc.  Phys.  Soc.  57,  147-159  (1945). 

37.  M.  J.  Forrestal,  D.  B.  Longcope,  and  F.  R.  Norwood,  “A  Model  to  Estimate  Forces  on 
Conical  Penctrators  into  Dry  Porous  Rock,”  J.  Appl.  Mech.  48,  25-29  (1981). 

38.  A.  L.  Florence,  P.  R.  Gefken,  L.  Seaman,  D.  R.  Curran,  and  D.  A.  Shockey, 

“Computational  Models  for  Armor  Penetration,”  SRI  International  Technical  Report  on 
Project  8521,  prepared  for  Alliant  Tcchsystems,  Inc.,  on  Contract  No.  635848-FN  (January, 
1992). 

39.  R.  Hill,  The  Mathematical  Theory  of  Plasticity  (Oxford  University  Press,  London,  1950). 

40.  Y.  M.  Gupta,  D.  D.  Keough,  D.  F.  Walter,  K.  C.  Dao,  D.  Henley,  and  A.  Urweider, 
“Experimental  Facility  to  Produce  and  Measure  Compression  and  Shear  Waves  in  Impacted 
Solids,”  Rev.  Sci.  Instrum.  51(2),  183-194  (1980). 


67 


ADAPTATIONS  TO  EXISTING  GAS  GUNS  TO  INCREASE 
IMPACT  VELOCITY  AND  FLYER  PLATE  DIAMETER 


ADAPTATIONS  TO  EXISTING  GAS  GUNS  TO  INCREASE 
IMPACT  VELOCITY  AND  FLYER  PLATE  DIAMETER 


by 

DAVID  C.  ERLICH 

Poulter  Laboratory 
SRI  International 
333  Ravenswood  Ave. 
Menlo  Park,  Ca  94025 

for 


41st  Meeting  of  the  Aeroballlstlc  Range  Association 

October  22,  1990 
San  Diego,  California 


ABSTRACT 


Adaptations  were  recently  wJe  to  the  two  light  gas  guns  at  SRI  International's 
Poulter  Laboratory  to  improve  their  versatility  as  research  tools.  An  alternative  breech  plug 
configuration  was  created  foi  the  smaller  (2.5-in. -diameter)  gas  gun,  resulting  in  a  nearly 
50%  increase  in  impact  velocity,  and  a  new  projectile/ target/recovery  system  was  designed 
and  built  for  use  in  the  larger  (4-in.-diameter)  gas  gun  that  allows  planar  impact  of  plates  up 
to  12  in.  in  diameter  at  relatively  low  velocities.  This  paper  will  briefly  describe  the  two 
gas  guns,  discuss  the  motivations  for  making  the  changes,  then  show  in  detail  the  design 
and  operation  of  the  new  configurations. 
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2.5-IN.  GAS  GUN:  BREECH  MODIFICATION 
FOR  HIGHER  VELOCITIES 


DESCRIPTION 

SRI  International's  2.5-in.  gas  gun  has  been  in  operation  for  about  30  years  and 
was  last  modified  in  1977.  This  modification  involved  slotting  the  barrel  and  adding  to 
each  projectile  a  key  that  rides  in  the  slot,  thus  preventing  projectile  rotation  and  allowing 
planar  impact  of  a  slanted  flyer  and  target  with  negligible  tilt. 

The  standard  wraparound  breech  configuration  (shown  in  Figures  1  and  2a)  has  the 
projectile  acting  as  the  rapid-opening  breech  valve.  With  the  projectile  latched  in  place,  the 
barrel  ahead  of  and  behind  it  is  evacuated  and  the  reservoir  is  pressurized  to  a  maximum  of 
=6000  psi.  To  activate  a  launch,  the  projectile  is  unlatched  and  high-pressure  gas  from  the 
reservoir  is  allowed  through  the  evacuation  tube  into  the  region  behind  the  projectile, 
pushing  it  forward  and  unblocking  the  reservoir  opening. 

The  projectile  walls  must  be  relatively  strong  (a  minimum  wall  thickness  of  »  5  mm 
of  7075-T6  aluminum  in  the  region  of  the  reservoir  opening)  to  withstand  the  high  pressure 
of  the  gas  prior  to  launch;  and  therefore,  the  projectile  is  relatively  heavy  (from  -500  to 
1200  g).  The  results  of  the  549  tests  performed  from  1977  to  1989  are  shown  in  Figure  2, 
plotted  in  terms  of  impact  velocity  (V)  versus  gas  pressure  divided  hy  projectile  mass 
(P/M).  The  highest  impact  velocity  has  been  677  m/s,  using  a  custom  projectile  (with 
outrigger  for  slanted  impact)  designed  for  minimal  mass.  The  projectile  weighed  535  g, 
and  the  tests  had  a  gas  pressure  of  5980  psi  (P/M  of  1 1 .2). 
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Standard  wraparound  breech  configuration  for  2.5  in.  gas  gun. 
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Figure  2.  Schematic  of  alternative  2.5  in.  gas  gun  breech  configurations 


MOTIVATION  FOR  CHANGE 

Recently  there  has  been  interest  in  gas  gun  tests  of  high-strength  ceramics,  such  as 
silicon  carbide  (SiC),  boron  carbide  (B4C)  and  titanium  diboride  (TiBj),  that  have 
Hugoniot  elastic  limits  (HEL)  in  the  150-kbar  range.  One  such  test,  called  the  pressure- 
shear  test,*  involves  slanted  impact  to  produce  shear  stresses  strong  enough  to  yield  the 
material.  Back-surface  interferometry  is  used  to  characterize  the  dynamic  flow  properties 
of  the  ceramic . 

Impact  velocities  on  the  order  of  900-950  m/s  are  required  to  obtain  the  desired 
stresses.  Figure  3  shows  an  extrapolation  of  the  V-P/M  curves  to  1000  m/s,  using  a  power 
law  fit  to  the  lower-velocity  data.  The  P/M  value  would  need  to  increase  by  =150%  (from 
=1 1  to  =27)  to  obtain  the  50%  increase  in  V  over  the  highest  previously  attained  value  of 
677  m/s.  Because  we  are  already  at  the  upper  limit  of  gas  pressure  (=6000  psi),  we  must 
reduce  the  projectile  mass  by  60%  (from  =550  to  =220  g). 

After  we  examined  many  possibilities  for  redesigning  the  projectile  (including 
changing  the  composition  of  the  walls,  reducing  the  wall  thickness  and  filling  the  inside 
with  a  light  foam,  and  ribbing  the  walls),  it  became  clear  that  nothing  would  effect  such  a 
major  reduction  in  the  projectile  mass  as  long  as  the  projectile  needed  to  withstand  the 
6000-psi  gas  pressure.  We  concluded  that  we  had  to  redesign  the  breech  configuration  to 
eliminate  the  use  of  the  projectile  as  a  breech  valve.  We  strived  to  minimize  (1)  the  costs  by 
using  as  many  existing  breech  components  as  possible  and  (2)  the  time  required  to  switch 
back  and  forth  between  the  new  and  the  old  breech  configurations. 


*R.  J.  Clifton  and  R.  W.  Klopp,  "  Pressure -Shear  Plate  Impact  Testing,"  in  Metals  Handbook A  Vol.  8, 9th 
Edition,  ASM'  International,  Metals  Park,  Ohio  (1985). 
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NEW  HIGHER-VELOCITY  BREECH  CONFIGURATION 

The  newly  designed  and  implemented  higher- velocity  breech  configuration  is 
shown  in  Figures  2b  and  4.  A  solid  polycarbonate  cylinder  replaces  the  projectile  as  the 
breech  valve,  while  the  projectile,  a  new  ultralight,  thin-walled  design,  is  positioned  in  the 
barrel  ahead  of  the  reservoir  opening.  The  plastic  slug  is  attached  by  a  tube  to  a  piston  in  a 
cylinder  outside  the  barrel.  The  tube,  which  is  also  used  to  evacuate  the  region  behind  the 
projectile,  passes  through  a  tapered  high-strength  steel  breech  plug. 

The  modifications  were  relatively  inexpensive  to  implement,  because  only  two 
components  of  the  standard  breech  configuration  needed  to  be  replaced  by  machined  parts. 
The  barrel  insen  was  redesigned  to  include  a  shoulder  (against  which  the  plastic  slug  and 
the  projectile  are  positioned)  and  holes  for  the  latching  pins.  A  new  tapered  hreech  plug 
replaces  the  original  plug. 

When  the  gun  is  fired,  compressed  air  pushes  the  piston,  which  pulls  the  front  end 
of  the  plastic  slug  past  the  reservoir  opening.  At  that  point,  high-pressure  gas  from  the 
reservoir  flows  into  the  barrel,  driving  the  plastic  slug  backward  (funher  opening  the 
reservoir)  and,  after  enough  pressure  is  attained  to  shear  off  the  latching  pins  (a  few 
hundred  psi),  propelling  the  projectile  forward.  The  reservoir  opening  time  is  expected  to 
be  about  the  same  for  this  configuration  as  for  the  standard  wraparound  breech.  The  plastic 
slug  is  severely  deformed  and  is  eventually  stopped  by  the  taper  of  the  breech  plug. 

The  ultralight  projectile  design  is  shown  in  Figure  5.  The  material  is  still  7075-T6 
aluminum,  but  the  wall  thickness  has  been  reduced  to  a  mere  1.6  mm.  The  design  allows 
the  gas  to  flow  to  the  outside  as  well  as  the  inside  of  the  wall,  which  serves  simply  as  a 
guide  to  keep  the  front  of  the  projectile  (=5  mm  thick,  to  withstand  the  high  pressure)  in 
alignment.  The  mass  of  the  projectile  shown  in  Figure  5  is  about  235  g,  depending  on  the 
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Figure  4.  New  higher-velocity  breech  configuration  Id*  2.5  b.  ga$  gun 
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density  and  thickness  of  the  flyer  plate.  Further  reductions  can  be  made  by  reducing  the 
projectile's  length  and  cutting  holes  in  its  side  walls. 

Thus  far,  four  experiments  have  been  performed  with  the  new  breech 
configuration,  shown  in  Figure  6.  The  highest  velocity  attained  has  been  939  m/s,  using  a 
projectile  with  a  mass  of  223  g  (breech  pressure  was  6090  psi,  resulting  in  a  P/M  of 
27.2  g).  The  velocity  was  a  only  a  few  percent  below  that  predicted  for  the  same  P/M  by 
the  power  law  extrapolation  from  lower-velocity  data.  The  measured  impact  tilt  has  been  as 
low  as  0.6  mRad,  which  is  comparable  to  the  best  attained  using  the  standard  wraparound 
breech  configuration. 
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4-IN.  GAS  GUN:  NEW  SYSTEM  FOR  IMPACT 
OF  12-IN.-DIAMETER  PLATES 


DESCRIPTION 

SRI  International's  4-in.-diameter  light  gas  gun  has  been  in  operation  for  about  20 
years.  Its  =14-m-long  barrel  is  unslotted,  so  the  gun  is  suitable  for  nonnal  impacts  only. 
The  gun  produces  impacts  with  uniaxial  strain  pulse  durations  (before  impingement  of 
radial  release  waves  from  the  target  or  projectile  periphery)  as  high  as  10-15  ps,  depending 
on  the  target  material,  which  is  enough  for  most  laboratory  shock  wave  experiments. 

The  gun  has  three  different  breech  configurations  to  span  a  wide  range  in  impact 
velocity.  A  standard  wraparound  breech  system  (Figure  7a),  similar  to  that  described 
above  for  the  2.5-in.  gun  but  with  a  maximum  reservoir  pressure  of  3000  psi,  has  yielded  a 
peak  impact  velocity  of  approximately  875  m/s  with  standard  thick-walled  aluminum 
projectiles.  A  double-diaphragm  breech  (Figure  7b)  with  a  maximum  reservoir  pressure  of 
6000  psi  has  yielded  a  peak  impact  velocity  of  nearly  1 .5  km/s  using  custom  lightweight 
projectiles  of  nylon  and  foam. 

A  third  configuration,  dubbed  the  "Monkey’s  Fist”  (Figure  8),  is  used  for  impact 
velocities  below  about  150  m/s.  It  uses  the  wraparound  breech  reservoir,  but  the  projectile 
is  held  near  the  muzzle  end  of  the  barrel  at  a  distance  of  only  1  -  5  ft  from  the  target. 
Shortening  the  distance  of  projectile  travel  allows  us  to  use  a  higher  gas  pressure  to 
accelerate  the  same  mass  to  the  same  velocity,  thus  reducing  the  great  effect  of  frictional 
losses. 


The  Monkey’s  Fist  is  attached  by  a  strong  flexible  steel  cable  to  the  breech  plug.  Its 
hardened  steel  fingers  grip  the  back  end  of  the  projectile  while  gas  fills  the  reservoir  and  the 
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Figure  8.  4-in.  gas  gun  'Monkey's  Fist"  configuration  (or  lower  velocity  impacts. 


barrel  region  behind  the  projectile.  A  gas-driven  piston  is  activated  to  launch  the  projectile, 
loosening  the  grip  of  the  steel  fingers.  The  Monkey’s  Fist  system  allows  us  to  predict 
impact  velocities  as  low  as  50  m/s  with  less  than  10%  uncertainty  and  as  low  as  15  m/s 
with  less  than  20%  uncertainty.  In  comparison,  the  wraparound  breech  system  had 
velocity  prediction  uncertainties  of  50%  at  50  m/s  and  75%  at  15  m/s. 

MOTIVATION  FOR  CHANGE 

There  has  been  significant  interest  over  the  last  few  years  in  stress  gage  validation 
for  both  in  situ  and  underwater  stress  gage  packages.  These  packages  consist  of  one  or 
more  ytterbium  or  other  stress  transducers  encapsulated  within  a  metal  jacket.  A  uniaxial 
strain  pulse  of  30-50  pis  duration  may  be  required  to  bring  the  gage  and  its  jacket  into 
stress  equilibrium  with  the  surrounding  material.  Because  these  gage  packages  are 
designed  to  operate  in  regions  of  relatively  low  stress,  peak  stress  pulse  magnitudes  of 
interest  for  gage  validation  are  from  hundreds  of  bars  to  several  kilobars. 

Gage  validation  tests  at  SRI  International  have  been  performed  using  a  16-in.- 
diameter  barrel  fy-C^-gas-explosion  gun  located  at  the  remote  test  facility.  Although 
providing  more  than  enough  uniaxial  strain  pulse  durations  (flyer  and  target  plates  could  be 
as  large  as  24  in.  in  diameter)  for  gage  validation  tests,  the  facility  is  relatively  expensive  to 
operate  (the  16-in. -diameter  by  24-in. -long  projectiles  are  not  reusable).  Our  goal  therefore 
was  to  modify  the  4-in. -gas  gun  to  allow  the  planar  impact  of  1 2- in  .-diameter  plates.  A 
12-in. -diameter  is  enough  to  provide  adequate  uniaxial  strain  pulse  durations  for  many 
applications.  To  minimize  costs  for  both  implementing  these  modifications  and  using  the 
new  configuration  in  testing,  we  planned  to  use  as  many  of  the  existing  4-in.  gas  gun 
components  as  possible  and  to  design  the  system  so  that  the  projectile  could  be  reused. 
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LARGE  DIAMETER  PLATE  IMPACT  SCHEME 


The  new  4-in.  gas  gun  configuration  for  impacting  1 2-in. -diameter  plates  is  shown 
in  Figures  9  and  10.  The  projectile  consists  of  an  18-in.-long  shaft  that  fits  into  the  muzzle 
end  of  the  barrel,  a  12-in, -diameter  by  l-in.-thick  aluminum  pusher  plate  bolted  to  the 
shaft,  and  a  12-in. -diameter  flyer  plate.  Both  the  pusher  and  flyer  plate  are  positioned  in 
the  target  tank,  with  the  impact  end  of  the  flyer  plate  located  9  in.  from  the  target  face.  The 
projectile  shaft  and  pusher  plate  are  reusable  from  shot  to  shot,  while  a  new  flyer  plate  is 
bolted  or  glued  onto  the  pusher  plate. 

The  target,  whose  diameter  is  slightly  larger  than  that  of  the  flyer  plate,  is  attached 
to  a  target  holding  and  alignment  fixture  bolted  to  the  floor.  Velocity  and  tilt  pin  blocks 
may  be  attached  to  the  front  of  the  target. 

The  Monkey’s  Fist  is  used  to  hold  and  release  the  projectile  in  the  same  manner  as 
described  above  for  lower  velocity  tests.  Because  the  projectile  assembly  is  much  heavier 
(as  high  as  =*20  kg)  and  the  distance  to  impact  shorter  than  in  the  usual  Monkey’s  Fist 
arrangement,  higher  gas  pressures  are  required  to  attain  the  same  impact  velocity  in  the  new 
configuration.  So  the  steel  cable  between  the  Monkey’s  Fist  and  the  breech  plug  was 
replaced  by  a  series  of  10-ft-long  by  1 -in. -diameter  high-strength  steel  rods,  joined  together 
by  connectors  and  pretensioned  to  remove  the  slack  before  pressurization.  The  system 
was  designed  to  operate  up  to  a  pressure  of  1500  psi  (the  original  Monkey’s  Fist  system 
operated  up  to  about  500  psi). 

The  projectile  must  be  prevented  from  completely  exiting  the  barrel  after  impact  to 
ensure  that  the  projectile  shaft  and  pusher  plate  are  undamaged  after  eaeh  shot.  For  this 
purpose,  several  aluminum  honeycomb  disks,  approximately  1 1  in.  in  diameter  and  3  in. 
thick,  are  positioned  behind  the  target  (far  enough  away  to  avoid  affecting  the  wave 
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Figure  9.  4-In.  gas  gun  configurallon  for  planar  Impact  of  12-in.-diameter  plates:  breech  section. 


1 1 _ CnjahabEa  AhjmSnum  Honeycomb 

Velocity  Shock  Absorbers 


e  10.  4-In.  gss  cun  configuration  for  planar  Impact  of  12in.-diameter  plates:  target  section. 


propagation  in  the  target  during  times  of  interest)  to  act  as  shock  absorbers.  The  first  layer 
or  two  of  the  honeycomb  crushes  on  impact  and  stops  the  projcctilc/targct  mass  after  it  has 
traveled  several  inches,  leaving  a  substantial  length  of  the  projectile  shaft  within  the  barrel. 

Honeycomb  layers  are  available  in  several  different  strengths,  depending  on  the 
mass  of  the  projectile  and  the  intended  impact  velocity.  The  honeycomb  shock  absorbers 
are  attached  to  a  flange  on  the  target  tank  door,  and  the  momentum  from  the  stopped  impact 
is  transmitted  through  a  massive  reaction  frame  to  a  momentum  trap  bolted  to  a  concrete 
pad  in  the  floor. 

Thus  far,  seven  tests  have  been  performed  using  the  new  12-in. -diameter  plate 
impact  system.  Impact  velocities  have  spanned  a  range  from  1 3  to  42  m/s,  as  shown  in 
Figure  1 1.  The  impact  tilt  has  been  .bout  0.5  mRad.  The  tests  so  far  have  involved  gage 
package  validations  in  both  dry  and  saturated  limestone.  Future  tests  may  call  for  the  use 
ol  large  Helmholtz-typc  coils  on  both  sides  of  the  target  to  create  a  uniform  magnetic  field 
and  allow  the  use  of  particle  velocity  gages. 
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BREECH  PR ES U R E/PROJFO TILE  MASS  (psi/g) 


THE  STRENGTH  BEHAVIOR  OF  GRANULATED  SILICON  CARBIDE 
AT  HIGH  STRAIN  RATES  AND  CONFINING  PRESSURE 


The  strength  behavior  of  granulated  silicon  carbide  at  high  strain  rates 
and  confining  pressure 

R.  W.  Kfopp  and  D.  A.  Shockey 

SKI  International,  333  Pavenswood  Avenue .  Menlo  Park,  California  94025-3493 
(Received  12  August  1991;  accepted  for  publication  lb  September  1991) 

The  dynamic  Mohr-Coulomb  behavior  of  silicon  carbide  (SiC)  was  inferred  from  symmetnc 
pressure/shear  plate-impact  experiments  which  entail  planar  impact  of  two  SiC  plates 
inclined  at  15*  to  the  impact  direction.  The  transverse  velocity  of  the  free  rear  surface  of  the 
target  plate  was  recorded  using  a  laser  Doppler  veiocimeter  system,  and  the  experiments 
were  simulated  using  a  postulated  viscopiastic  constitutive  model  that  accounts  for 
comminution  and  dilatancy.  Model  parameters  were  vaned  until  the  computed  and  measured 
velocity  histones  agreed.  The  results  indicate  that  comminution  occurred  soon  after 
loading,  and  thus  the  experiment  measures  the  behavior  of  granulated  material  at  shear  strain 
rates  of  -105  s"1  and  mean  stress  ranging  from  1  to  9  GPa.  A  friction  coefficient  of 
0.23  was  obtained,  which  is  about  half  the  value  for  quasistatic  compression  of  precomminuted 
ceramic  reported  in  the  literature.  The  simulation  results  were  strongly  affected  by  the 
values  chosen  for  the  friction  coefficient  and  yield  strength  parameters;  changes  in  the 
dilatancy,  rate-sensitivity,  and  strain- hardening  parameters  had  a  lesser  effect. 


INTRODUCTION 

Examination  of  confined  ceramic  targets  that  have 
been  impacted  by  long  tungsten  rods  suggests  that  the 
comminution  and  granular  flow  behavior  of  the  ceramic 
governs  the  resistance  to  penetration. '  A  penetrator  can 
advance  in  a  ceramic  target  that  is  confined  on  all  sides  hv 
steel  plates  only  if  ceramic  material  at  the  leading  penetra¬ 
tor  surface  is  moved  out  of  the  penetrator  path.  Fracto- 
graphic  observations  of  targets  sectioned  after  ballistic  test¬ 
ing  indicate  that  this  occurs  by  the  crushing  of  a  small 
quantity  of  the  ceramic  within  about  two  penetrator  radii 
ahead  of  the  advancing  penetrator  and  the  subsequent  flow 
of  the  comminuted  material  lateral  to  and  then  opposite 
the  direction  of  attack.  Thus  the  shear  strength  of  commi¬ 
nuted  as  well  as  intact  material  contributes  to  penetration 
resistance.  Hydrocode  simulations  show  that  the  mean 
stress  in  this  zone,  which  we  shall  call  the  Mescal)  zone 
after  John  Mescih  who  first  deduced  its  existence,2*3  is  on 
the  order  of  the  confined  compressive  strength  of  the  ce¬ 
ramic.  For  silicon  carbide  the  confined  compressive 
strength,  i.e.,  the  Hugoniot  elastic  limit  (HEL),  is  around 
13  GPa.°  Simulations  also  indicate  that  equivalent  shear 
strain  rates  in  the  zone  are  on  the  order  of  104-105  s~*. 
Thus  the  penetrability  of  armor  ceramics  is  controlled  by 
their  resistance  to  comminution  and  fragment  flow  at  high 
loading  rates  and  at  high  compressive  mean  stress  levels. 

Computational  simulations  can  make  armor  material 
and  package  design  more  efficient  by  eliminating  tnal-and- 
error  ballistic  testing.  These  simulations  require  material 
models  that  reflect  the  micro  mechanisms  of  deformation 
and  failure  and  experimental  measurements  of  properties 
controlling  these  mechanisms.  However,  relatively  few 
techniques  are  available,  particularly  for  ceramics,  to  make 
measurements  of  comminution  and  granular  flow  proper¬ 
ties  under  penetration  conditions.  The  objective  of  this 
work  was  to  perform  pressure/shear  plate -impact  experi¬ 


ments  on  an  armor  ceramic  to  measure  the  comminution 
and  flow  properties  under  conditions  of  confining  pressure 
and  strain  rate  representative  of  those  at  the  tip  of  an  ad¬ 
vancing  penetrator. 

Up  to  now,  the  properties  of  armor  ceramics  have  been 
measured  using  quasistatic  and  dynamic  compression  of 
small  cylinders,  with  or  without  confining  stress,  or  using 
uniaxial  strain  impact  of  plates.  Heara  and  Cline6  investi¬ 
gated  the  quasistatic  mechanical  behavior  of  BeO,  AL03, 
and  AIN  cylinders  in  axial  compression  with  a  radial  con¬ 
fining  pressure  and  constructed  stress-strain  curves  and 
Stress-difference  versus  confinmg-stress  curves.  Lankford7 
and  Johnson  er  j/.5  performed  similar  tests  on  armor  ce¬ 
ramics,  but  used  a  Hopkinson  bar  apparatus  to  achieve 
strain  rates  near  1C3  s  ”  \  However,  neither  of  the  cylinder 
compression  techniques  can  reach  the  mean-stress  levels 
thought  to  exist  ahead  of  an  advancing  penetrator.  The 
radial  confining  stress  is  always  well  below  the  axial  failure 
stress,  and  thus  the  mean  stress  is  dominated  by  the  axial 
stress.  Posnest  examination  of  the  specimens  showed  fail¬ 
ure  occurred  by  faulting,  which  on  a  macroscopic  level  is  a 
different  failure  mechanism  than  comminution  and  flow  of 
fragments.4  In  the  Hopkinson  bar  tests,  reported  strengths 
are  for  the  onset  of  failure,  i.e.,  plastic  strains  are  negligi¬ 
ble. 

Johnson  ei  al.s  performed  quasistatic  compression  tests 
on  comminuted  SiC  that  was  prepared  by  ball  milling.  The 
data  are  reported  til  terms  of  Von  Mises  equivalent  com¬ 
pressive  stress  a Y  as  a  function  of  mean  stress  p .  and  when 
converted  to  equivalent  shear  stress  rY  =  cY/  ^3,  as  a  func¬ 
tion  of  mean  stress  and  fitted  with  a  straight  line,  a  friction 
coefficient  of  0.46  is  indicated. 

Brar  and  Bless0  reported  a  bar-impact  test  for  measur¬ 
ing  the  high-strain-rate  compressive  strength  of  ceramics 
and  glass.  Two  slender,  constant-diameter  bars  were  im¬ 
pacted  axially  at  a  relative  velocity  of  100-300  m/s,  and 
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the  impact  was  photographed  using  a  high-speed  framing 
camera.  The  photographs  indicated  the  velocity  of  defor¬ 
mation  and  fracture  waves  traveling  along  the  target  bar.  A 
piezoresistanr  stress  gauge  sandwiched  in  the  target  bar 
several  diameters  from  the  impact  end  provided  a  measure 
of  the  dynamic  compressive  strength.  Unfortunately,  mean 
stresses  and  strain  rates  in  these  tests  are  not  accurately 
known,  and  so  the  results  are  difficult  to  interpret. 

Gust,  Holt,  and  Roycc,10  Kipp  and  Grady,4  and 
Johnson  ex  al ?  performed  uniaxial-strain  shock  compres¬ 
sion  and  release  experiments  on  several  ceramics,  in  which 
they  recorded  free  surface  or  window  interface  particle  ve¬ 
locities.  Kipp  and  Grady  constructed  constitutive  models 
by  iteratively  varying  a  proposed  stress-strain  path  until 
computer  simulations  using  the  path  reproduced  the  mea¬ 
sured  velocities.  Kipp  and  Grady  were  able  to  simulate 
their  measurements  on  SiC  by  using  an  eiasuc-plastic 
model  adapted  from  metals.  However,  because  theirs  were 
uniaxial  strain  experiments,  plastic  strains  were  small,  be¬ 
ing  on  the  order  of  and  because  velocities  were  high, 
mean  stresses  were  higher  than  those  in  the  zone  ahead  of 
a  penetrator,  which  are  on  the  order  of  30  GPa. 

Klcpp  and  Shockey11  reported  shear  stress-strain 
curves  for  aluminum  nitnde  (AIN)  at  a  mean  stress  of  11 
GPa  and  a  shear  rate  of  1.3 X  105  s“  1  for  shear  strains  to 
These  values  were  obtained  using  the  high -strain -rate 
pressure/shear  plate-impact  technique, 12  in  winch  an  in* 

'  clined  flyer  plate  carrying  a  thin  Inver  of  specimen  material 
impacts  a  similarly  inclined  anvil  plate  such  that  the  spec¬ 
imen  is  compressed  and  sheared  between  the  plates.  Before 
waves  arrive  at  the  middle  of  the  plates  from  the  periph¬ 
eries,  and  if  the  flyer  and  anvil  plates  remain  dasne  while 
the  specimen  yields  and  flows,  an  analysis  very  similar  to 
that  for  the  Hopkinson  bar  applies,  and  a  stress-strain 
curve  can  be  directly  obtained  for  conditions  like  those  at 
the  tip  of  a  penetrator.  Unfommately,  this  "sandwich’*  test 
cannot  be  used  for  the  strongest  ceramics,  since  no  mate¬ 
rials  that  will  remain  elastic  (except  perhaps  diamond)  are 
available  for  the  flyer  and  anvil  plates. 

To  enable  evaluation  of  very  strong  ceramics,  the  high- 
strain-ratc  pressurc/shear  plate-impact  experiment  must 
be  modified.  By  impacting  one  plate  of  the  specimen  ma¬ 
terial  directly  with  a  second  plate,  i.c.,  a  symmetric 
pressure/shear  experiment  (Fig.  1),  both  longitudinal  and 
shear  waves  axe  produced  in  the  plate?.,  and  material  be¬ 
havior  can  be  inferred  from  velocity-ome  profiles  recorded 
in  the  intenor  of  the  anvil  or  at  its  rear  surface.  The  trans¬ 
verse  velocity  profile  due  to  shear  waves  is  especially  indic¬ 
ative  of  the  resistance  of  a  material  to  deformation  at  high 
mean-stress  levels  and  loading  rates.  Conditions  such  as 
those  at  the  tip  of  a  penetrator  can  be  produced.  A  draw¬ 
back  of  the  symmetric  pressure/shear  experiment  is  that 
computer  simulations  are  required  to  infer  material  prop¬ 
erties  from  the  recorded  velocity  histones  ft  r  all  but  the 
very  simplest  constitutive  benaviors  n 

This  paper  desenbes  two  symmetric  pressure/shcar 
plate- impact  experiments  on  SiC  ceramic.  We  first  desenbe 
the  experimental  procedure  and  the  results.  The  interpre¬ 
tation  of  the  results  by  numerical  simulations  is  then  pre- 


s*c  Py*  PKt* 


FIG.  1.  The  symmetric  pressure/ shear  experiment. 

sented  and  a  material  modei  is  proposed.  Finally,  the  re¬ 
sults  and  model  are  discussed  in  relation  to  data  obtained 
by  other  researchers. 

PRESSURE/SHEAR  EXPERIMENTS 

Two  experiments  were  performed  usi^s  plates  pre¬ 
pared  from  a  single  tile  of  sintered  a-SiC  supplied  by  Car¬ 
borundum.  The  density  of  the  SiC  was  3.2  g/cm3  and  the 
grain  size  approximately  3  fxm.  Both  flyer  and  anvil  plates 
were  disks  oominully  55.4  mm  in  diameter,  and  the  flyers 
were  3  mm  thick.  The  anvil  plates  for  experiments  l  and 
2  were  5,06  an*.<  7.51  mm  thick,  respectively.  The  impact 
and  rear  faces  of  rhe  anvils  were  initially  lapped  to  a  flat¬ 
ness  of  2.5  /im,  with  a  roughness  of  about  20  nm  root  mean 
square  (rms).  Then  the  rear  face  of  the  anvils  was  further 
polished  to  a  roughness  of  about  6  nm  rms  to  improve 
reflectivity. 

The  flyers  were  mounted  at  an  angle  of  15#  oo  the 
leading  end  of  a  6-in.-long  aluminum  projectile  and 
launched  from  SRI’s  2.5-in.-bore  slotted  gas  gun  at  0.913 
and  0.715  mmy/is,  A  key  on  the  projectile  engages  the  slot 
to  prevent  rotation  and  loss  of  parallelism.  Before  launch, 
the  anvil  plates  were  inclined  at  an  angle  of  15*  and  aligned 
to  be  parallel  with  the  flyers  within  0.2  mrud.  Four  gold 
contacts  that  had  been  vapor -deposited  onto  the  periphery 
of  the  plates  indicated  parallelism  and  triggered  recording 
instruments.  The  impact  velocity  was  measured  using 
spaced  shorung  pins  located  at  the  gun’s  muzzle.  Impact 
occurred  in  an  evacuated  tank  to  reduce  the  lubricating 
effect  of  gas  trapped  between  the  flyer  and  anvil. 

In  previous  pressure/sbear  experiments  performed  on 
metals,  the  transverse  motion  of  the  anvil  rear  surface  has 
been  measured  using  either  a  transverse  displacement 
interferometer14  (TDI)  or  two  VISARs.15  The  TDI  re¬ 
quires  a  diffraction  grating  on  the  anvil  rear  surface,  but 
because  of  the  fragility  of  polymer  diffraction  gratings  and 
the  difficulty  of  etching  gratings  directly  into  the  ceramic, 
the  TDI  could  not  be  used  for  these  experiments.  Although 
the  VISAR  technique  could  probably  be  used,  we  chose 
instead  to  adapt  a  laser  Doppler  velocimeter  (LDV)  sys¬ 
tem,  commonly  used  in  fluid -mechanics  experiments,  to 
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FIG.  2.  The  laser  Doppler  veioameter  (LDV)  system. 


these  shock-physics  experiments  (Fig,  2).  In  this  system 
two  parallel  laser  beams  pass  through  an  objective  lens  that 
focuses  them  and  causes  them  to  intersect  and  form  an 
ellipsoidal  measurement  volume.  Within  the  measurement 
volume,  the  two  intersecting  beams  interfere  and  form  a 
pattern  of  equispaccd  bright  and  dark  planar  fringes.  As  a 
light-scattering  object  moves  transversely  through  the  mea¬ 
surement  volume,  the  scattering  intensity  is  modulated  hy 
the  fringes.  When  the  object  moves  at  a  constant  velocity, 
the  scattered  light  intensity  history  is  sinusoidal.  Since  the 
fringe  spacing  is  fixed  by  the  laser  wavelength  and  beam 
intersection  angle,  the  modulation  frequency  gives  the 
speed  of  the  scattering  object.  The  fringes  in  the  measure¬ 
ment  volume  arc  parallel  to  a  plane  that  contains  the  bi¬ 
sector  of  the  beam  intersection  angle  and  is  perpendicular 
to  the  plane  of  the  intersecting  beams.  Thus,  if  the  bisector 
is  made  to  be  perpendicular  to  the  surface  of  a  solid  target, 
the  modulation  frequency  will  be  insensitive  to  the  normal 
motion  of  the  target.  Another  attractive  feature  of  the  sys¬ 
tem  is  that  the  results  are  independent  of  the  condition  of 
the  anvil  rear  surface,  provided  that  this  surface  scatters 
light  toward  the  detector  and  the  reflectivity  changes  rela¬ 
tively  slowly  so  as  not  to  be  confused  with  intensity 
changes  due  to  motion  through  the  fringes. 

In  our  system  the  514.4-nm  argon-laser  beams  were 
spaced  52  mm  apart  before  entering  the  objective.  The  ob¬ 
jective  had  a  focal  length  of  800  mm,  giving  a  fringe  spac¬ 
ing  of  7.9  /im.  The  measurement  volume  was  about  150  /im 
in  diameter.  The  scattered  light  was  collected  through  the 
center  of  the  objective  lens  and  convened  to  an  electrical 
signal  with  a  photomultiplier  tube  (PMT).  We  found  that 
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EApenmcm  1 

Experiment  2 

Impact  velocity  lmm/^s) 

0913 

o.7i ; 

Skew  angle  (deg) 

15 

Tilt  (rarad) 

36 

06 

Anvil  thickness  (mm) 

5.06 

2.51 

the  rear  surface  of  the  target  (which  had  a  6-nm-rms  fin¬ 
ish)  had  enough  scratches  and  pits  to  scatter  light  and 
provide  measurable  signals. 

The  electrical  signal  from  the  PMT  was  amplified  and 
split  so  that  half  could  be  recorded  directly  and  the  other 
half  processed  with  a  10-MHz  high-pass  filter  to  remove 
meaningless  large,  lower-frequency  excursions.  Signals 
were  recorded  on  250-MHi-bandwidth  analog  oscillo¬ 
scopes.  The  transverse  velocity  was  determined  by  measur¬ 
ing  the  times  of  the  peaks  and  valleys  of  the  sinusoidal 
signal.  We  divided  the  displacement  herweea  consecutive 
peaks  (7.9  pm)  by  the  time  between  those  peaks  to  deter¬ 
mine  the  velocity  history. 

Parameters  of  the  two  experiments  are  given  in  Table 
I,  and  the  oscilloscope  traces  and  transverse  velocity  his¬ 
tories  are  shown  in  Fig.  3.  Where  the  sinusoid  is  clear  and 
has  relatively  constant  frequency,  such  as  in  experiment  2, 
the  velocity  is  known  to  within  5%.  Where  it  is  not 
clear — such  as  near  wave  arrival  and  near  times  when  scat¬ 
tering  objects  enter  or  leave  the  measurement  volume — the 
uncertainty  is  greater.  For  example,  the  uncertainty  is 
nearly  100%  between  0.80  and  0.90  p$  in  Fig.  3(a)  (ex¬ 
periment  1),  where  the  LDV  sinusoid  contains  a  segment 
of  high-frequency,  low-amplitude  noise,  presumably 
caused  hy  the  absence  of  light  scattercrs  in  the  measure¬ 
ment  volume  Similarly,  the  uncertainty  is  at  least  50%  at 
1,03  /is,  where  the  phase  of  the  LDV  signal  suddenly 
jumps,  presumably  as  a  result  of  scatterers  entering  or  leav¬ 
ing  the  measurement  volume.  We  believe  that  in  experi¬ 
ment  1  the  transverse  velocity  was  around  0.17 ±0.03 
mm /fis  and  probably  steady  between  0.4  and  1.2  /is. 

ANALYSIS 

We  used  an  iterative  computational  procedure  to  de¬ 
duce  the  stresses  and  strains  within  the  anvil  plates  from 
the  measured  rear -surface  transverse  velocities.  In  this  pro¬ 
cedure  we  propose  a  material  model  with  associated  pa¬ 
rameters  and  simulate  the  experiments  with  it.  If  the  cal¬ 
culated  transverse  particle  velocity  histories  match  those 
measured  in  the  experiments,  then  we  suggest  that  the 
model  is  appropriate”— at  least  for  the  conditions  in  the 
experiments — and  correctly  infers  the  histories  of  stress 
and  strain  within  the  anvils.  This  iterative  procedure  has 
been  often  applied.  Examples  include  the  Taylor  test,  in 
which  a  constitutive  model  is  varied  until  the  computed 
deformed  rod  profiles  match  those  observed;16  shock-com¬ 
pression  and  release  experiments,  in  which  a  loading  path 
is  varied  until  computed  and  experimental  particle  veloci¬ 
ties  agree/  and  previous  symmetric  pressurc/shear  exper¬ 
iments  on  metals.17 
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FIG.  3.  LDV  oscilloscope  recordi  and  meaju red  and  computed  transverse 
velocity  histories. 


For  this  project  we  simulated  the  pressure/shear  ex¬ 
periments  following  the  analysis  of  Abou-Sayed  and  Clif¬ 
ton  for  aluminum.17  However,  we  expanded  their  analysis 
to  include  dilatancy,  comminution,  and  mean-stress -depen¬ 
dent  yielding,  because  these  phenomena  are  thought  to  be 
important  in  the  flow  of  ceramics.  We  consider  the  speci¬ 
men  as  an  infinite  piate  normal  to  the  x  ajus,  with  x  =  0 
being  the  impact  face,  and  the  y  axis  extending  in  the  di- 
rccnon  of  the  imposed  transverse  velocity.  Considering 
symmetry  and  that  the  motion  is  independent  of  the  y  and 
z  coordinates,  there  is  no  motion  in  the  z  direction  In¬ 
direction  velocity  w  =  0),  and  the  shear  stresses  and 

must  he  zero.  Furthermore,  if  the  dilatancy  is  isotropic, 
then  we  can  show  that  the  in-plane  normal  stresses  are 
equal  ~  o^).  Thus  there  are  five  unknowns:  the  nor¬ 
mal  velocity  ti,  the  y  transverse  velocity  vt  the  normal 
stresses  and  and  the  shear  stress  (TXjr  The  governing 
equations  include  the  momentum  equation,  the  compati¬ 
bility  equation,  and  a  constitutive  relation.  These  equations 
can  be  combined  to  yield  the  five  necessary  to  find  the 
unknowns. 

We  use  the  constitutive  framework  developed  by 


Perzyoa,11  in  which  the  total  strain  rate  is  additively  de¬ 
composed  into  elastic  and  plastic  parts  and  the  plastic  part 
is  determined  by  the  oveistress  (the  difference  between  the 
applied  equivalent  stress  and  the  current  yield  equivalent 
stress).  The  resulting  constitutive  equation  has  the  form 

Uv  v  __  <3/ 

£>j=  *tj  -  £  Si/Tkk  +  WT'ripJ))  )  ^T,  (  1  ) 

where  itJ  is  the  total  strain  rate,  v  is  Poisson’s  ratio,  E  is 
Young's  modulus,  6 is  the  KLronecker  delta,  p 
=  —  0**/3  is  the  negative  of  the  mean  stress  or  hydro¬ 
static  pressure,  d  is  the  equivalent  piasuc  strain  rate,  r  is 
the  effective  shear  stress,  y*  is  the  equivalent  plastic  shear 
strain,  D  is  the  amount  of  comminution,  and  /  is  a  flow 
potential.  The  effective  stress  is  given  in  the  usual  way  in 
terms  of  the  de  via  cone  stress,  StJ  =  a 

f- (2) 


and  the  equivalent  plastic  shear  strain  is  defined  as 


{<P)dt 


The  bracket  function  (<f>)  is  defined  such  that 


0,  for  <£<0, 
<f>t  for  <£>0. 


(3) 


(4) 


The  first  two  terms  on  the  right-hand  side  of  the  constitu¬ 
tive  equation  ( 1 )  determine  the  elastic  part  of  the  strain 
rate  and  thr  third  term  gives  the  plastic  strain  rate. 

Following  Rudmcld  and  Rice,19  we  define  our  flow 
potential  to  be 


/= 


\n  ,  i  _ 


(5) 


where  {3  is  the  so-called  dilatancy,  which  is  defined  as  the 
increment  in  volumetric  strain  resulting  from  a  unit  incre¬ 
ment  of  plastic  shear  strain. 

The  experimental  results  of  Heard  and  Cline6  and 
Lankford7  for  AiN  and  A12Oj  and  of  Johnson  et  ai 5  for 
Si C  show  rather  weak  elevation  of  the  overstress  with  in¬ 
creasing  strain  rate.  Thus  we  choose  the  overstress  to  he 
proportional  to  the  logarithm  of  the  plastic  strain  rate, 
such  that 


(6) 


where  ry  is  the  equivalent  shear  yield  stress  and  y 5,  r0,  and 
m  are  constants. 

We  choose  the  equivalent  shear  yield  stress  to  depend 
on  the  mean  stress,  amount  of  comminution,  and  strain 
hardening: 


h  , 

Ty=F{pJ))Ty o-r  j 


(7) 


where  F  is  the  mean  stress-damage  function,  h  is  the  tan¬ 
gent  modulus  of  a  uniaxial  stress-strain  curve,  and  r ^  is 
the  initial  shear  yield  stress  at  atmospheric  pressure.  As  the 
material  changes  from  an  intact  to  a  tulJy  comminuted 
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state,  D  goes  from  zero  to  one  and  F  simultaneously 
changes  from  a  form  appropriate  to  undamaged  material 
( Fu )  to  a  form  appropriate  to  comminuted  material 

{FDy. 


F=(l  -D)Fu~  DFo-  (8) 

The  function  Fu  is  cho.sen  so  that  the  yield  stress 
smoothly  increases  with  mean  stress  up  through  the  HEL 
point: 


lt  for  p<p& 


for  p>p& 

(9) 


where  and  ^HEL  are  the  shear  yield  stress  ana  mean 
stress,  respectively,  at  the  Hugoniot  elastic  limit  and  p0  is  a 
constant. 

The  function  FD  is  chosen  so  that,  below  a  limiting 
mean  stress,  yield  stress  varies  linearly  with  mean  stress 
according  to  Mohr-Couiomb  behavior.  Above  the  limiting 
mean  stress,  the  yield  stress  is  assumed  to  be  constant.  The 
limiting  mean  stress  is  chosen  such  that  the  shear  strength 
cannot  exceed  the  theoretical  strength  (chosen  to  be  G/2Q, 
where  G  is  the  shear  modulus): 


PI* 


20  V 


fOr  P  P mJl> 

for 


(10) 


where  p  is  the  Mohr-Coulomb  friction  coefficient  and 
Pm*i  ~  G/20fl. 

The  amount  of  comminution,  Dy  is  initially  zero  and 
increases  toward  one  if  the  square  root  of  the  second  in¬ 
variant  of  the  total  stress  I2  ~  exceeds  some  con¬ 

stant  limit  7^.  The  rate  of  increase  of  D  is  proportional  to 
the  amount  by  which  I2  exceeds  7^.  Thus 


D^—  .M(i.j),  (11) 

Pc  T:o 

where  p  is  the  density,  c,  is  the  longitudinal -wave  speed, 
and  r0  is  a  time  constant.  The  choice  of  Iz  as  the  driving 
force  for  comminution  ls  based  on  the  notion  that  even 
purely  hydrostatic  loadings  comminute  ceramic.  Hydro 
static  loadrng  causes  comminution  because  neighboring 
grains  often  have  different  crystallographic  orieutauons 
and  thus  have  different  strain  fields  as  a  result  of  amsot- 
ropy.  Cracks  form  when  stresses  necessary  to  enforce  com¬ 
patibility  of  these  strain  fields  across  grain  boundaries  ex¬ 
ceed  some  failure  stress. 

Constitutive  equations  ( 1 )-( 1 1 )  are  combined  and,  to¬ 
gether  with  the  momentum  and  compatibility  equations, 
arc  rewritten  as  equations  along  characteristics.  vV'e  solved 
the  equations  using  a  modificanon  of  the  method  of  Ran- 
ganath  and  Clifton  70  The  boundary  conditions  for  the  an¬ 
vil  plates  art  that  ..ic  rear  surface  is  tractiou  free  and  the 
impact  face  has  an  imposed  velocity  of  one-half  the  impact 
velocity. 


RESULTS 

The  solid  lines  in  Figs.  3(a)  and  3(b)  are  the  free- 
surface  transverse  velocity  histories  for  the  two  experi¬ 
ments.  computed  using  the  parameters  in  Table  II.  The 
computed  histones  follow  the  measured  histories  reason¬ 
ably  well.  Figure  4  shows  the  mean  stress  p  as  a  function  of 
position  through  the  anvil  (as  a  percentage  of  the  anvil 
thickness)  at  various  times.  Similar  plots  for  equivalent 
plastic  strain  rate  6,  equivalent  plastic  strain  p7,  and  com¬ 
minution  D  are  shown  in  Figs.  5-7.  The  travel  of  the  lon¬ 
gitudinal  waves  through  the  anvil  to  the  rear  surface  and 
back  is  readily  apparent.  At  0.913  mm/ps9  a  large  portion 
of  the  specimen  sees  a  mean  stress  of  about  9.2  GPa  before 
dropping  to  roughly  half  that  level.  At  0.715  nun /us,  the 
mean-stress  wave  front  is  less  sharp,  rising  to  a  peak  level 
of  about  7  GPa  before  dropping  by  half.  Less  of  the  spec¬ 
imen  reaches  the  peak  mean-stress  levels  in  the  0."15- 
mm/ps  simulation.  From  Fig.  5  it  is  clear  that  pLasnc  shear 
strain  rates  of  the  order  of  3xl05  $~l  arc  achieved,  al¬ 
though  these  rates  are  sustained  for  only  a  short  duration 
at  any  given  location.  Figure  6  shows  that  the  cumulative 
effect  of  those  strain  rates  are  plastic  shear  strains  of  3.3% 
and  2.6%,  respectively,  at  the  impact  faces,  and  6.5%  and 
4.5%  peak  levels,  respectively,  for  0.913-  and  0.715- 
mm/ps  impacts.  Figure  7  shows  that  the  comminution 
level  D  rapidly  approaches  1  in  the  interiors  of  the  anvils 
for  hoth  impacts,  and  comparison  with  Fig.  5  shows  that 
this  occurs  before  the  peak  in  strain  rate.  Thus  the  majority 
of  the  deformation  takes  place  in  highly  comminuted  ma¬ 
terial. 

DISCUSSION 

The  results  of  experiment  2  (0.715  mm/ps)  are  rela¬ 
tively  unambiguous.  The  transverse  velocity  time  record 
[Fig.  3(b)]  looks  similar  to  those  reported  previously  for 
rocks  and  metals,  although  the  initial  plateau  followed  by  a 
later  rise  is  more  reminiscent  of  the  double-wave  structure 
seen  in  longitudinal  waves.  This  double-wave  structure 


TABLE  II.  Simulation  parameters 


Density  (kg/m') 

P 

3  .ZO* 

Lon ptudmaJ- wave  speed  (mrn/Ats) 

'l 

u.r 

Shear-wave  speed  (mm/jis) 

c: 

«* 

Initial  yield  stress  (GPa) 

r>o 

10.0 

Fncuon  coefficient 

p 

0.23 

Shear  modulus  (GPa) 

G 

164.0* 

Poisson  s  ratio 

V 

G.U* 

Mean  stress  at  HEL  (GPaj 

Phel 

6.3h 

Shear  stress  at  HEL  (GPaJ 

EL 

64 

Square  root  of  second  stress  invariant 

7  iim 

3.1 

at  damage  limn  (GPa) 

Strain -rate  constant  (s-1) 

Vo 

3.0 

Shear  stress  constant  (MPa) 

rD 

1.5 

Strajmratc  power 

m 

2.0 

Pressure  constant  (GPa) 

Pn 

1.0 

Time  constant  (ns) 

h 

1.0 

Hardening  modulus  (GPa) 

.h 

10.0 

•Reference  21. 
''Reference  4 
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FIG.  4.  Computed  mean  stress  as  a  foe  cron  of  position  within  the  anvil. 

may  result  from  strong  coupling  between  the  shear  and 
longitudinal  waves,  which  would  provide  evidence  that  the 
flow  stress  of  comminuted  SIC  is  controlled  by  the  mean 
stress.  We  could  investigate  this  issue  further  by  perform¬ 
ing  experiments  in  which  the  anvil  thickness  varies  widely. 
If  true  double  waves  are  present,  the  lengths  of  the  plateaus 
separating  the  waves  will  vary  with  the  anvil  thickness. 

Because  the  transverse  velocity  record  of  experiment  2 
has  a  plateau  structure  and  the  velocity  is  well  below  that 
expected  if  the  SiC  were  clastic,  we  infer  that  the  SiC 
yielded.  If  the  SiC  had  remained  elastic,  the  transverse 
velocity  would  have  jumped  to  a  level  of  0.185  mm//*s 
upon  shear-wave  arrival.  The  normal  stress  would  have 
been  12.0  GPa,  the  in-plane  normal  stress  1.95  GPa,  and 
the  shear  stress  2.09  GPa,  The  equivalent  normal  stress 
would  have  been  10.4  GPa.  This  value  is  below  the  value  at 
the  reported  HHL,°  i.e.,  e=  13  GPa.  Thus  we  can  conclude 
independently  of  the  simulations  that  the  SiC  appears  to 
have  yielded  at  an  equivalent  stress  near,  hut  below,  that  of 
the  HEL.  We  do  not  know  if  the  SiC  would  yield  at  even 
lower  stress  levels.  Rosenberg,  Yeshurun.  and  Brandon22 
have  reported  evidence  of  plasticity  in  Coons  AD-85  alu¬ 
mina  loaded  below  the  HEL.  whereas  Rosenberg,  Brar, 
and  Bless23  found  no  evidence  of  plasticity  in  porous  or 
noriporous  B4C  loaded  below  the  HEL. 

Plots  of  shear  stress  as  a  function  of  mean  stress  from 
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FTG.  5.  Computed  equivalent  plisoc  strain  rate  u  t  function  of  position 
within  the  tnvil. 

the  simulations,  that  is,  loading  paths  for  material  elements 
in  the  anvil*  do  not  present  a  very  clear  picture  of  the 
response  of  SiC  to  hugh-strain-rate,  high-mean -stress  load¬ 
ing,  because  the  plastic  strain,  strain  rate,  and  comminu¬ 
tion  ah  vary  simultaneously.  A  clearer  picture  emerges  if 
we  plot  the  response  of  the  model  when  constant  values  for 
plastic  strain,  strain  rate,  and  comminution  are  assumed. 
The  response  of  the  model  at  two  extremes,  ( 1 )  a  strain 
rate  of  10J  s  “ 1  and  strain  of  0  and  (2)  a  strain  rate  of  I05 
s  * ]  and  strain  of  0.04,  both  for  a  comminution  of  1.0,  is 
shown  in  Fig.  8.  Also  included  in  this  plot  are  the  data  of 
Johnson  et  al 5  for  intact  and  powdered  material  and  that 
of  Kipp  and  Grady  for  intact  ceramic.4  From  this  plot  we 
conclude  that  the  inferred  response,  i.e.,  the  flow  strength 
at  a  given  pressure,  is  softer  than  even  the  response  of 
powdered  material  at  quasistatic  loading  rates. 

In  the  simulations  the  response  is  only  transiently  af¬ 
fected  by  the  properties  of  the  intact  ceramic,  as  can  he 
seen  by  the  steep  front  and  high  speed  of  the  comminution 
wave  as  it  traverses  the  anvil  (Fig.  7).  The  effects  of  the 
intact  ceramic  are  transient  because  the  comminution  rime 
constant  (r0)  is  only  1.0  ns.  The  response  is  dominated  by 
the  friction  coefficient  (/x  =  0.23)  and  the  combination  of 
parameters  that  determine  the  overall  flow  stress  (r0 
=  L3  MPa,  y0  =  3.0  s~ l,  and  m  =  2).  Varying  the  dila- 
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FIG.  6.  Computed  equivalent  pLtuoc  shear  strain  as  a  function  of  position 
within  the  snviL 

tan.  •  (£),  rate-sensitivity  (y^m),  and  strain-hardening 
( h )  parameters  while  keeping  the  friction  coefficient  and 
overall  flow  stress  constant  has  a  lesser  effect  on  the  results, 
A  short  relaxation  time  seems  reasonable  when  one  con¬ 
siders  the  high  wave  speeds  and  short  crack  propagation 
distances  involved  in  comminution.  The  value  fi  =  0.23  is 
below  the  range  0.4<iu<0.9  quoted  by  Rudnicki  and  Rice19 
for  mineral-on-nuneral  sliding  and  is  one-half  the  value 
obtained  by  Johnson  et  at }  on  precomminuted  Si C. 

We  suggest  that  the  reason  the  flow  strength  at  a  given 
pressure  is  less  than  that  measured  in  cylinder-compression 
and  uniaxial  strain  plate-impact  tests  is  that  plastic  strains 
in  the  latter  tests  are  of  the  order  of  1  %  or  less,  whereas 
plastic  strains  in  our  tests  were  4. 5%— 6,59c.  Such  large 
differences  in  plastic  strains  imply  large  differences  in  the 
amount  of  fragment  sliding  that  has  occurred  and  the 
amount  of  concomitant  interffagment  debris  that  has  been 
generated.  Interffagment  debris  might  tend  to  act  as  rollers 
between  the  fragments.  Shear  localization  in  the  pres- 
sure/shear  experiments  might  also  account  for  the  low  in¬ 
ferred  strength.  Indeed,  numerical  overflow  of  the  plastic 
strain  rate  occurred  in  the  simulations  when  impact  veloc¬ 
ities  much  above  0.913  mm//ts  were  input. 

The  reason  that  the  flow  strength  we  inferred  is  less 
than  that  measured  quasi  statically  on  hall-milled  SiC  hy 
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FIG.  7.  Computed  comminution  *s  t  function  of  position  within  the  anviL 

Johnson  et  al5  might  be  the  differences  in  fragment  sizes, 
shear  localization,  strain-rate  effects,  or  the  fact  that  our 
mean  stresses  were  an  order  of  magnitude  higher.  At  high 
strain  rates,  fragment  interfaces  could  melt  if  frictional 
heat  builds  up  there.  At  high  mean  stresses,  interfragment 
debris  would  form  more  readily  than  at  low  stresses. 


FIG.  8.  Equivalent  flow  stress  is  i  (unction  of  mean  stress  for  SiC. 
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Our  simulations  show  clearly  that  the  SiC  specimen 
experiences  a  range  of  mean  stresses,  strain  rates,  and  shear 
stresses,  all  of  which  contribute  to  the  wave  profile  re¬ 
corded  at  the  free  surface.  Thus  a  simulation  that  ade¬ 
quately  reproduces  experimental  results  requires  a  material 
model  chat  is  effective  over  a  relatively  broad  range  of 
conditions.  We  have  not  vaned  the  model  parameters 
widely  and  thus  do  not  know  whether  other  combinations 
of  model  constants  would  give  the  same  results  or  whether 
a  Mchr-Coulomb  modei  is  the  only  appropriate  one.  To 
reduce  the  uncertainty,  more  experiments  with  different 
anvil  thicknesses,  impact  velocities,  and  skew  angles  arc 
required.  However,  we  believe  the  model  reasonably  rep¬ 
resents  the  stresses  within  the  anvil,  because  the  model 
reasonably  reproduces  velocity  measurements  at  the  anvil 
rear  surface.  However,  extrapolation  from  the  test  condi¬ 
tions  should  be  done  with  caution. 

The  results  of  experiment  1  (0.913  mm  //is )  show  rwo 
disturbing  features.  First,  the  shear  wave  arrived  late  so 
that  the  free  surface  did  not  attain  a  high  transverse  veloc¬ 
ity  until  well  after  it  was  expected  to.  The  lace  arrival 
suggests  that  the  impact  faces  initially  slipped  hecause  they 
were  iubneated  by  gas  that  leaked  past  the  projectile.  Two 
other  pieces  of  evidence  that  suggest  gas  leaked  ahead  of 
the  projectile  are  the  large  (36  mrad)  hit  and  the  fact  that 
the  target  was  moving,  albeit  slowly,  before  impact.  In  the 
absence  of  gas  leakage  (presumably  the  case  for  experi¬ 
ment  2),  we  believe  that  the  impact  faces  would  not  slip. 
One  reason  is  that  the  maximum  possible  ratio  of  shear 
stress  to  normal  stress  for  a  15*  impact,  obtained  assuming 
the  plates  remained  elastic,  is  0,17,  which  is  below  our 
inferred  fnenon  coefficient  of  0.23.  Furthermore,  in  exper¬ 
iments  on  metals,  slippage  only  occurred  for  impact  angles 
above  23"  (Ref  12)  and  was  associated  with  a  sudden  drop 
in  transmitted  shear  stress. 

The  second  disturbing  feature  of  experiment  1  is  the 
oscillation  of  the  transverse  velocity.  If  the  flyer  and  anvil 
remained  elastic  or  i I  the  longitudinal  and  shear  waves 
were  otherwise  uncoupled,  the  velocities  normal  and  trans¬ 
verse  to  the  anvil  rear  surface  would  remain  constant  after 
the  initial  transient  arrived,  because  the  anvil  is  thinner 
than  the  flyer.  We  have  not  seen  unsteady  behavior  re¬ 
ported  for  similar  tests  on  metals  or  rocks.  It  is  probably  an 
artifact  of  scattering  objects  moving  in  and  out  of  the  LDV 
measurement  volume.  This  effect  is  known  to  give  false 
indications  of  tnrbulence  in  fluid  ancmomecry.24’  The  be¬ 
havior  may  be  an  indication  of  suck -slip  of  the  impact  faces 
or  of  heterogeneous  cracking  of  the  anvil.  If  moving  cracks 
emitted  waves  that  traversed  the  measurement  volume  of 
the  LDV,  then  the  transverse  velocity  would  have  been 
perturbed.  If  we  had  measured  the  transverse  velocity  at 
two  locations,  wc  could  tell  whether  the  oscillations  were 
due  to  a  phenomenon  that  was  uniformly  disenbuted  over 
a  plane  parallel  to  the  plates,  such  as  stick-shp,  or  to  phe¬ 
nomena  that  were  heterogeneously  distributed,  such  as 
fracturing  and  LDV  artifacts. 
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CONCLUSIONS 

Comminution  of  sintered  silicon  carbide  under  shock 
loading  occurs  at  a  stress  level  slightly  below  the  HEL.  The 
strength  of  the  highly'  comminuted  material  increases 
monotonically  with  mean  stress  (confining  pressure)  in  the 
1-9  GPa  range  ana  is  three  to  six  times  less  than  the 
strength  of  uncomminuted  material.  The  strength  of  com¬ 
minuted  material  at  high  shear  strain  rates  (10*  $~l)  is 
significantly  lower  than  at  quasisutic  rates  and  may  result 
from  rate  dependence  of  interfragment  friction.  The  fric¬ 
tion  coefficient  of  0.23  inferred  from  the  pressure/shear 
experiments  reported  here  is  smaller  than  that  reported  for 
quasistatic  mineral-on-mineral  sliding  by  a  factor  of  2—*. 

A  laser  Doppler  velocuneter  is  a  viable  cost-effective 
alternative  to  the  transverse  displacement  interferometer 
and  the  V1SAR  for  measunng  the  transverse  morion  of  the 
rear  surface  of  shock-loaded  ceramic  targets. 

A  Simple  nonassociated  continuum  modei  that  in¬ 
cludes  dtfatancy.  rate  dependence,  and  a  one-parameter 
measure  of  the  extent  of  comminution  reasonably  describes 
the  measured  results. 

Using  the  experimental  measurements  and  model,  we 
have  demonstrated  a  method  for  deducing  flow  properties 
of  a  fragment  bed  under  high  strain  rate  and  confining 
pressure. 
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Summary— The  mechanism  by  which  a  long  rod  penetrates  a  steel -encased  ceramic  block  was 
sought  by  performing  impact  experiments  at  a  range  of  vclocihes  and  examining  the  fracture  and 
dctormation  in  the  recovered  largcis  and  impaaors  The  ko>  processes  are  the  crashing  of  a  small 
volume  of  ceramic  adjacent  to  ihe  leading  surlace  ul  the  advancing  penetramr.  and  the  subsequent 
flow  of  I  he  fine  fragments  lateral  to  and  then  opposite  the  direction  of  atiack.  The  results  suggest 
that  nonconvcniionul  material  properties  such  as  the  dynamic  compressive  failure  energy  and  the 
friction,  flow  and  abrasive  properties  of  the  finely  fragmented  material  govern  the  penetration 
resistance  of  confined  ceramics.  This  understanding  of  pcnciraiion  mechanisms  can  be  used  to 
guide  development  of  specialized  tests  and  failure  models  lo  measure  pcrlincnt  material  properties 
and  to  predict  penetration  behavior,  respective!). 


INTRODUCTION 

When  ceramic  plaies  arc  used  as  overlays  or  incorporated  as  a  layer  within  conventional 
monolithic  steel  armor,  the  ballistic  protection  is  significantly  enhanced  [I],  This  finding 
has  encouraged  the  use  of  ceramics  as  a  component  of  advanced  armoi  structures  and  lias 
motivated  researchers  to  identify  materials  and  structural  configurations  that  maximize 
ballistic  performance. 

The  design  of  current  ceramic  armors  is  based  predominantly  on  empirical  ballistic 
performance  data.  Test  firings  using  the  threat  of  interest  (long  rod,  shaped  charge  jet.  or 
small  caliber  ammunition)  are  conducted  against  armor  specimen."  in  which  ceramic 
component  parameters,  such  as  type  of  ceramic,  thickness,  and  spacing  of  plates,  arc  varied 
independently  and  systematically,  The  combination  of  parameters  that  produce  maximum 
hallistic  protection  is  determined  from  test  results  and  used  to  design  the  armor  package. 
This  procedure  is  lengthy  and  expensive,  and  because  the  number  of  potentially  influential 
material  and  geometry  variables  is  large,  a  comprehensive  test  matrix  is  not  practical  and 
so  it  is  doubtful  that  armoi  packages  affording  optimal  protection  are  attained. 

A  more  efficient  procedure  is  to  combine  experiments  with  computational  simulations 
of  experiments.  Computations  of  the  ballistic  behavior  for  various  impact  conditions  and 
target  geometries  indicate  which  target  parameters  are  important  and  suggest  com  btnations 
of  target  parameters  that  will  give  favorable  performance.  A  limited  matrix  of  ballistic 
experiments  based  on  these  guidelines  are  performed.  The  results  are  used  to  modify  the 
models  and  the  code;  then  a  second  generation  of  computational  simulations  is  conducted 
and  used  lo  design  a  second  set  of  test  filings.  This  procedure  is  repeated  until  an  armor 
package  with  acceptable  performance  isohtamed.  The  reduction  in  the  number  of  required 
test  firings  reduces  substantially  the  expense  and  time  required  to  attain  a  suitable  armor, 
and  the  understanding  gained  promotes  optimum  armor  design. 

Implementation  of  the  iterative  computational,  experimental  procedure,  however,  re¬ 
quires  reliable  models  for  the  microfailurc  behavior  of  ceramics  and  penctrators  under 
penetration  conditions  These  models  should  be  based  on  an  understanding  of  the  failure 
phenomenology  during  penetration.  In  particular,  the  material  properties  governing 
penetration  resistance  must  be  known.  Unfortunately,  conventional  material  properties 
such  as  fracture  toughness,  strength,  and  hardness  correlate  poorly  with  penetration 
behavior  (_2.  3],  suggesting  that  under  the  complex,  high-rate,  multiaxi.il  load  produced 
by  the  penetrator  one  or  inoic  noncoiiventional  material  properties  control  penetration. 


DosAi.n  A  Shotkiy  or  nl 


Tlic  goal  of  the  work  reported  here  was  to  establish  the  failure  phenomenology  of 
confined  ceramic  targets  and  impacting  long  rods  during  penetration,  and  to  deduce  the 
ceramic  properties  governing  penetration  resistance.  Our  approach  was  to  perform  impact 
experiments  oil  confined  ceramic  specimens  at  several  velocities  to  produce  damage  ranging 
from  incipient  to  severe.  Very  early  stages  of  damage  were  studied  by  performing 
expel orients  with  low  velocity  spherical  particles.  From  fractographic  and  metal logruplue 
examination  of  the  targets  and  rods  after  impact,  we  inferred  the  failure  mechanisms  and 
speculate  on  the  properties  controlling  penetration  behavior.  The  study  was  aimed  at 
ceramics  as  a  class  of  materials  and  sought  a  qualitative  understanding  of  penetration 
phenomenology  to  provide  the  basis  for  a  computational  model:  therefore,  experiments 
were  performed  on  a  variety  of  ceramics  and  details  of  microstruclure  and  mechanical 
properties  of  the  individual  materials  are  not  presented. 

ROD  IMPACT  EXPERIMENTS 

Tungsten -nickel-iron  rods,  0.770  mm  in  diameter  and  7.70  mm  long,  having  hemi¬ 
spherical  noses  were  accelerated  in  a  powder  gun  to  desired  velocities  in  the  0.8- K4  km  s“  * 
range  and  caused  to  impact  a  stcel-encascd  block  of  ceramic  at  approximately  zero  degrees 
obliquity  (big.  1  Counties  investigated  included  AUO*.  $iC\  BaC  and  TiB^ 

After  the  impact  event,  the  target  assemblies  were  carefully  removed  from  the  mounting 
fixture  and  placed  on  the  floor  of  the  bunker  with  the  impact  surface  facing  up.  A 
two-component  epoxy  was  poured  into  the  hole  on  the  front  cover  plate  to  fill  the  crater 
and  run  into  the  cracks  in  the  ceramic.  This  procedure  was  intended  to  strengthen  the 
fractured  ceramic  block  so  that  the  confining  siccl  could  he  removed  and  the  ceramic  block 
could  be  scetione<l  without  crumbling  of  the  fractured  ceramic.  We  found  in  later 
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Fig.  2.  Rear  surface  crack  pattern  (a)  and  crack  patterns  on  cross-sections  (b)  in  targets  of  B4C 

and  TiB,. 

experiments  that  impacted  targets  were  often  strong  enough  to  be  disassembled,  sectioned, 
and  even  wafered  without  infiltration  of  epoxy. 

The  front  surface  of  the  ceramic  block  and  the  inside  surface  of  the  steel  cover  plate 
showed  a  starburst  pattern  of  linear  markings  radiating  outward  from  the  impact  site. 
These  markings  were  produced  by  ceramic  and  tungsten  rod  fragments  originating  near 
the  leading  edge  of  the  penetrator. 

All  ceramics  tested  exhibited  qualitatively  similar  cracking  patterns;  however,  the 
numbers  and  sizes  of  cracks  generally  differed.  The  crack  patterns  on  the  rear  surfaces  of 
a  B4C  and  TiB2  specimen  impacted  at  0.8  km  s“l  are  shown  in  Fig.  2(a).  Two  main  types 
of  cracks  are  evident:  circular  cracks  and  radial  cracks.  Three  to  four  dominant  circular 
cracks  were  observed  in  both  ceramics;  however,  approximately  50%  more  radial  cracks 
were  produced  in  the  B4C  (23  cracks)  than  in  the  TiB2  (15  cracks).  Fracture  damage  was 
heaviest  in  both  materials  in  a  region  directly  beneath  the  impact  site. 
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A  diamond  saw  was  used  to  eul  through  the  emter  centers  aiii  it  plane  containing  the 
impact  direction.  Cutting  the  ceramic  block  was  difficult.  The  diamond  wheel  wore  out 
quickly  and  slopped  cutting  about  halfway  through  the  hloek.  The  wheel  needed  to  he 
redressed  several  times  hefore  the  sectioning  was  completed. 

Orthogonal  views  of  the  craters  and  crack  patterns  were  obtained  on  the  as-sawn, 
unpolished  section  surfaces.  These  surfaces  for  a  B AC  and  a  TiB:  specimen  are  shown  in 
Fig.  2(b).  Crater  size  was  greater  in  the  BaC.  Material  directly  heneath  the  crater  in  the 
TiBc  specimen  remained  intact,  prohahly  because  of  the  Irtaxial  compressive  stress  state 
and  the  higher  compressive  yield  strength  of  Til3:.  At  higher  velocities  resulting  in  deeper 
rod  penetration,  the  material  m  advance  of  the  tip  of  the  penetrator  was  crushed  to  a  fine 
powder.  Similar  damage  patterns  were  ohserved  m  ALO*  and  SiC  targets. 

The  views  in  F  ig.  2{b|  show  that  the  circular  cracks  in  Fig.  2(a)  arc  traces  of  cracks  that 
ran  outward  at  an  angle  from  the  impact  site.  The  resulting  cone  configuration  corresponds 
to  the  I  iertzian  cracks  observed  in  the  particle  impact  experiments  reported  in  the  following 
section  and  ohserved  by  others  in  ceramics  and  glasses  under  static  indentation  and  particle 
impact  [4,5].  The  lateral  cracks  lying  roughly  parallel  to  tiie  surface  apparently  formed 
after  the  cone  cracks,  since  they  are  discontinuous  across  the  cone  cracks.  Thus,  lateral 
cracks  may  he  produced  hv  tensile  stress  waves  reflecting  from  the  specimen  houndarics 
or  hv  later  unloading  of  the  target.  Radial  cracks  are  not  revealed  on  cross-sections 
containing  the  impact  direction. 

In  several  instances,  tungsten  fragments  were  observed  lodged  between  faces  of  Hertzian 
cone  cracks.  That  these  fragments  were  not  moved  into  the  cracks  hy  the  sectioning 
operation  was  confirmed  tv  computed  tomograph^  results  that  showed  fragments  in  cracks 
in  unseetioned  specimens.  Such  observations  suggest  that  dehrisemanatingfrom  the  eroding 
end  of  the  penetrating  rod  can  have  a  significant  forward  velocity  component.  These 
observations  support  the  premise  of  Hauver  [6],  who  ohserved  tungsten  fragments  in 
advance  of  the  penetrator  in  X-radiographs  of  ceramic  blocks  during  impact  by  tungsten 
rods. 

Loose  fragments  produced  in  the  impact  experiments  were  collected  for  examination. 
The  size  distribution  of  the  collected  fragments  was  determined  by  a  sieve  analysis  (placing 
the  recovered  fragments  on  the  topmost  of  a  stack  of  successively  finer  screens  and  vibrating 
the  stack  for  1  h).  We  separated  ceramic  fragments  from  penetrator  fragments  by  passing 
a  strong  magnet  over  the  sieved  fragments  and  extracting  the  slightly  magnetic  tungsten 
alloy  fragments  from  the  ceramic  dchris. 

Figure  3  shows  the  fragments  of  SiC'  retained  on  screens  with  various  mesh  openings. 
Fragment  shape  did  not  vaiv  substantially  with  fragment  size:  aspect  ratios  ranged  from 
I  -ai  3.  Fracture  was  predominantly  transgranular  rather  than  intergranular.  Differences 
in  fragment  size  distributions  for  the  four  ceramic  materials  were  small  over  the  entire  si/e 
range.  Figure  4  compares  the  distributions  for  three  ceramics  in  the  2  40  am  size  range. 

In  higher  velocity  experiments,  in  which  the  rod  penetrated  60-100  mm.  the  crushed 
ceramic  material  produced  at  the  leading  edge  of  the  rod  flowed  around  and  behind  the 
rod  closing  the  hole  made  hy  the  rod  (Fig.  5).  So  well  consolidated  were  these  fine  fragments 
that  no  fragments  were  loosened  during  sectioning  ami  individual  fragments  were  not  easily 
discernible  by  high  magnification  examination.  Hardness  and  scratch  tests  indicated 
strong*  h*  of  the  compacted  powder  comparahle  to  that  of  the  umnipacted  material.  In  I  ig. 
5,  the  penetrator  has  stopped  just  short  of  the  rear  confinement  plate. 

The  debris  from  several  targets  was  searched  for  distal  portions  of  the  penetrator.  Distal 
portions  ranging  in  length  from  3  20  mm  were  found  [Fig.  6(a)].  Proximal  ends  had  cither 
a  mushroom  shape  or  a  sharpened  pencil  shape.  We  speculate  that  a  penetrator  tip  may 
alternate  shapes  hoi  ween  that  of  a  muslnoom  and  a  pencil  point  several  times  during 
the  penetration  process.  Initially,  we  expect  plastic  deformation  of  the  leading  edge  to 
produce  a  mushroom  shape.  The  mushroom  zone  then  shears  away  on  roughly  a  45 
conus,  producing  a  pencil  point.  T  he  pencil  point  then  deforms  plastically  and  the  tip 
acquires  a  mushroom  shape  again:  this  mushroom  shape  becomes  unwieldy  and  .shears  lo 
a  pencil  point.  This  alternating  shape  change  continues  until  penetration  ceases.  The 
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Fig.  3.  Fragments  from  a  SiC  target  retained  on  screens  with  different  mesh  openings. 


proximal  end  surfaces  of  all  recovered  penetrators  were  faceted  and  gouged,  suggestive  of 
shear  failure.  The  lateral  surfaces  were  unscored. 

Polished  and  etched  cross-sections  on  planes  containing  the  rod  axis  revealed  that  the 
tungsten  particles  in  this  sintered  alloy  retained  their  original  roughly  spherical  shape 
everywhere  except  near  the  proximal  failure  surface.  Adjacent  to  the  failure  surface,  the 
tungsten  particles  were  greatly  elongated,  often  to  aspect  ratios  of  5  or  more.  Particle 
distortion  decreased  with  distance  from  the  failure  surface,  rather  gradually  (over  a  distance 
of  about  4  mm)  in  the  mushroomed  region.  Thus,  the  distribution  of  deformed  tungsten 
particles  provides  a  map  of  the  plastic  strain  field  in  a  penetrator. 

Tungsten  fragments  extracted  magnetically  from  the  debris  ejected  from  the  impact 
surface  were  also  examined  with  a  scanning  electron  microscope.  Failure  surfaces  and 
etched  cross-sections  suggested  that  fragment  formation  was  by  localized  shearing  of  the 
microstructure,  in  accord  with  observations  on  distal  penetrator  ends.  Tungsten  particle 
distortion  in  the  fragments,  however,  decreased  abruptly  (usually  within  about  300  /*m) 
from  the  surface  of  fragments  [Fig.  6(b)]. 
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We  computationally  simulated  the  experiment  depicted  in  Fig.  1  using  a  two-dimensional. 
Lagrangian  finite  difference  code,  l  uc  results  provided  an  estimate  of  the  distribution  and 
time  variation  of  the  stresses  and  strains  produced  in  the  ceramic  target  by  the  impacting 
long  rod  before  failure  occurred,  and  assisted  in  the  interpretation  o!  the  fritctographic 
observations. 


PART  1C  LL  IMPACT  EXPERIMENTS 

l.n\s  velocity  particle  impact  experiments  were  performed  to  study  incipient  stages  of 
impact  damage.  The  evolution  of  fracture  damage  was  established  in  hot  pressed  (HP) 
silicon  nitride  by  accelerating  single  solid  spheres  of  tungsten  carbide  (WC|  or  steel  onto 
the  polished  surfaces  of  small  plate  specimens  of  HP  SixNa  at  a  90  angle  [7].  Particles 
were  accelerated  to  velocities  from  Jb-30N  m  s'1  by  filling  the  gun  chamber  with  nitrogen 
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Fig  6.  Distal  portions  of  tungsten  alloy  rods  recovered  from  ballistic  experiments  (a)  and  polished 
and  etched  cross-scction  showing  deformation  of  the  microstructurc  (b). 


gas  to  various  pressures,  then  suddenly  releasing  the  nitrogen  by  rupturing  a  disk.  The 
diameters  of  the  WC  spheres  were  1.6  and  2.4mm;  the  steel  spheres  were  2.4mm  in 
diameter.  Impact  and  rebound  velocities  were  recorded  with  photomultipliers.  Photo¬ 
multiplier  records  also  showed  whether  particles  remained  intact  or  fragmented  after  impact. 
The  specimen  fracture  damage  was  studied  by  optical  and  scanning  electron  microscopy 
of  impact  surfaces  and  polished  cross-sections  normal  to  the  impacted  surfaces. 

The  impact  tests  caused  several  kinds  of  cracks,  small  craters  and  fragmentation  in  the 
target  plates  and  eventually  plastic  deformation  or  fracture  of  the  impacting  spheres.  Targets 
sustained  no  damage  at  impact  velocities  below  17  ms'1,  at  which  point  ring  cracks 
appeared.  As  impact  velocity  increased,  the  damage  progressed  to  cone  cracks,  an  inelastic 
impression,  radial  cracks,  lateral  cracks,  and  median-vent  cracks.  Ring  cracks,  as  shown 
in  Fig.  7.  arc  circumferential  cracks  that  extend  less  than  a  millimeter  beneath  the  surface. 
As  the  impact  velocity  increased,  more  and  longer  ring  cracks  formed  [Fig.  7(b)].  The 
ring  cracks  are  similar  to  the  Hertzian  ring  cracks  lormcd  under  quasi-static  loading  [8], 
The  surface  ring  cracks  that  start  approximately  normal  to  the  specimen  surface  veer 
outward  at  various  angles  up  to  about  75  from  the  vertical  to  become  Hertzian  cone 
cracks  [9],  As  the  velocity  increase*,  additional  cone  cracks  form  both  inside  and  outside 
the  cxistingdamage  umbrella, and  the  innermost  cone  grows  several  millimeters  in  depth. 

An  inelastic  impression  and  radial  cracks  [F ig.  7(c)]  seemed  to  form  at  the  same  time 
in  the  failure  sequence.  As  the  impression  deepened  with  increasing  velocity,  the  radial 
cracks  grew  in  both  si/e  and  number,  although  only  a  small  number  {8  or  9)  of  the  radial 
cracks  grew  to  scvcial  millimeters  [Fig.  7(d)j. 
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Fig.  7.  Cracks  on  the  surface  of  HP  SisN4  caused  by  impact  of  2.4  mm  diameter  tungsten  carbide 
spheres  at  velocities  of  (a)  19.5  ms-1  (b)  46.2  ms-’  (c)  97.7  ms-’  and  (d)  159  ms-’. 


Figure  8  shows  the  internal  damage  and  the  extent  of  growth  of  the  various  cracks  below 
the  specimen  surface.  The  nucleation  and  growth  sequence  of  the  ring/cone  cracks  is 
illustrated  in  Figs  8(a)  and  8(b).  Under  increasingly  severe  impacts,  cone  cracks  seemed 
to  cease  growing;  instead,  two  new  types  of  cracks  were  created,  as  shown  in  Fig.  8(c). 
Lateral  cracks  nucleated  internally  near  the  contact  center  and  ran  approximately  parallel 
to,  and  eventually  veered  toward,  the  impact  surface  of  the  specimen.  Vertical  cracks 
initiated  internally  in  the  region  within  the  innermost  cone  crack.  These  latter  penny-shaped 
cracks  are  similar  to  the  median-vent  cracks  observed  by  Evans  and  Wilshaw  in  quasi-static 
indentation  experiments  on  ZnS  [8],  Observations  with  polarized  light  showed  that  a  zone 
of  densely  microcracked  material,  approximately  spherical  in  shape,  was  formed  beneath 
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Fig.  8.  Seciional  views  of  subsurface  cracking  pallern  in  HP  Si*N4  impacted  by  2  4  mm  diameter 
sled  spheres  at  velocities  of  56.4  ms" 1  (a)  and  231  ms" 1  fbi  and  by  a  2.4  mm  diameter  tungsten 

carbide  sphere  at  231  ms"  1  icy 


the  contact  area.  Zinc  sulfide  impacted  by  0.4  and  0.8  mm  WC  spheres  exhibited  a  similar 
microcracked  zone  [10]. 

Impacting  steel  spheres,  which  are  softer  than  tungsten  carbide,  caused  only  ring  and 
cone  cracks  and  introduced  little  additional  damage  above  300  m  s'1,  at  which  velocity 
the  particle  failed  by  plastic  deformation.  This  limit  on  the  damage  inflicted  on  the  ceramic 
occurs  because  the  particle  cannot  exert  a  pressure  on  the  ceramic  greater  than  the  particle's 
yield  strength.  Since  the  yield  strength  of  the  steel  is  less  than  the  pressure  required  for 
inelastic  deformation  ol  the  ceramic  surface,  higher  velocity  impacts  only  result  in  more 
deformation  of  the  steel  sphere. 

FAILURE  PHI’NOMKNOLOGV  Oh  THE  PENETRATION  PROCESS 

The  picture  of  the  penetration  process  that  begins  to  emerge  from  these  observations 
and  consideration  of  the  initial  stress  history  is  as  follows. 

Caleulalional  simulations  of  a  tungsten  alloy  rod  impacting  a  target  as  in  Fig.  I  at 
1600  m  s~  1  show  that  at  the  instant  of  impact,  a  shock  wave  with  an  amplitude  of  several 
hundred  kbars  is  generated  at  the  impact  site,  Radial  divergence  and  plastic  flow  and 
fracture  in  the  steel  cover  plate  quickly  and  drastically  reduce  the  stress  so  that  the  strength 
of  the  shock  that  passes  into  the  ceramic  is  below  the  Hugomot  elastic  limit.  Thus,  the 


Donald  A.  Shockey  ef  al. 


212 

initial  shock  wave  is  not  expected  to  condition  the  ceramic.  A  steady-state  ramp  wave 
follows  the  shock,  loading  the  ceramic  material  at  the  tip  of  the  penetrator  to  a  maximum 
compressive  stress  of  about  50-60  kbars.  The  ceramic  initially  resists  the  stress  in  the  ramp 
wave  and  exerts  large  stresses  on  the  tungsten  rod,  which  may  deform,  fracture  or  be 
deflected. 

Ceramics  are  substantially  stronger  in  compression  than  in  tension,  and  consequently, 
the  tensile  strength  of  the  ceramic  is  exceeded  at  the  impact  surface  near  the  rod  periphery 
and  tensile  fracture  begins  to  occur  soon  after  impact.  The  stress  fields  in  the  ceramic  are 
initially  elastic,  and  the  largest  tensile  stresses  are  in  the  radial  direction  (the  Boussinesq 
stress  field).  Therefore,  the  cracks  that  form  (normal  to  the  direction  of  maximum  principal 
stressj  arc  ring  cracks  concentric  about  the  impact  site.  These  cracks  are  shallow  cracks, 
extending  initially  only  a  millimeter  or  so  beneath  the  ceramic  surface.  Upon  continued 
loading,  however,  several  ring  cracks  continue  to  grow  and,  following  the  paths  normal 
to  the  direction  of  the  principal  tensile  stress,  assume  angled  trajectories  25-75T  outward 
from  the  initial  normal-to-thc-surfacc  direction.  Thus,  several  large  Hertzian  cone  cracks 
extend  through  the  ceramic  block,  intersect  the  specimen  surfaces,  and  cause  structural 
failure  of  the  target. 

Up  to  this  point,  the  stress  fields  and  the  fracture  response  arc  elastic.  But  as  the  rod 
continues  to  advance,  the  compressive  strength  is  exceeded  in  material  directly  beneath 
the  penetrator.  Microcracking  occurs  in  a  shallow  zone  near  the  penetrator  tip,  and  the 
stress  field  changes  in  character.  The  principal  tensile  stresses  are  now  in  the  circumferential 
direction,  and  a  new  type  of  tensile  crack  is  invoked — six  to  twelve  large  radial  cracks  run 
outward  from  the  impact  site  like  spokes  from  a  hub  [9].  These  cracks  intersect  the  impact 
surface  and  may  extend  to  all  specimen  boundaries,  resulting  in  strength  degradation  and 
eventual  structural  failure  of  the  target, 

A  third  crack  type,  lateral  cracks,  form  beneath  the  im  pact  surface  and  propagate  roughly 
parallel  to  it,  probably  during  i  loading.  These  cracks  intersect  cone  cracks  and  radial 
cracks,  thereby  providing  the  oithogonal  surfaces  necessary  for  fragment  formation. 
Cratering  results  when  these  large  fragments  are  liberated  from  the  vicinity  of  the  impact  site. 

The  tensile  cone,  radial,  and  lateral  cracks  do  not  provide  an  easy  path  for  penetrator 
advance  and  hence  do  not  assist  the  penetration  process  directly.  Intact  and  laterally 
confined  ceramic  remains  in  the  penetrator  path  despite  the  presence  of  these  tensile  cracks, 
and  this  material  must  be  moved  out  of  the  way  for  the  penetrator  to  advance.  This  occurs 
by  pulverization  of  the  ceramic  material  in  a  shallow  zone  immediately  ahead  of  the 
penetrator  and  the  subsequent  flow  of  this  material  laterally  and  opposite  to  the  impact 
direction,  processes  that  occur  under  large  compressive  and  shear  stresses. 

Thus,  the  development  of  a  densely  inicrocracked  zone  in  a  ceramic  directly  ahead  of 
the  impactor  is  a  prerequisite  for  penetration.  Insight  into  how  this  zone  forms  can  be 
gleaned  from  the  observations  of  Hagan  and  co-workers  [11,12]  of  damage  zones  in 
soda-lime  glass  and  zinc  sulfide  produced  by  quasistatic  indentation.  These  workers 
observed  a  curvilinear  grid  of  narrow,  fault-like  flow  lines  beneath  indentations,  and  voids 
and  microcracks  at  many  of  the  nodes  in  the  grid.  Flow  in  the  polycrystallinc  ZnS  occurred 
by  slip  and  twinning  within  the  grains  and  by  grain  boundary  sliding;  voids  formed  when 
grain  boundary  displacements  became  large  either  along  the  flow  lines  or  at  flow  line 
intersections. 

The  finely  fragmented  material  at  the  leading  edge  of  the  penetrator  wants  to  occupy  a 
larger  volume  (i.e.  dilation),  but  expansion  is  resisted  by  the  confinement  of  the  slccl 
encased  ceramic  block.  The  resulting  increase  in  pressure  makes  fragment  flow  more  difficult 
and  adds  to  the  resistance  exerted  on  the  penetrator.  The  tensile  cracks  may  assist 
penetration  indirectly  by  reducing  the  level  of  constraint  on  the  pulverized  material,  thereby 
allowing  easier  flow  of  the  material  out  and  away  from  the  advancing  penetrator,  but  the 
main  resistance  to  penetration  is  probably  coupled  to  the  flow  characteristics  of  highly 
comminuted  ceramic  powder.  The  cracking  pattern  in  the  ceramic  targei  envisioned  during 
the  steady-state  phase  of  the  penetration  process  is  depicted  in  Fig.  9. 

As  the  ceramic  particles  flow  across  the  leading  surface  of  the  penetrator,  they  erode 


Failure  phenomenology  of  confined  ceramic  targets 


27) 


impaci  Surface  Ejec;a 


RA- 7525-0 

Fiti.  9.  Cracking  pattern  inceramic  targetsdurine  thcsieadv-slatc  phase  of  the  penetration  process, 


the  rod,  shortening  and  eventually  consuming  it  as  the  rod  moves  through  the  ceramic. 
No  scoring  or  erosion  of  the  sides  of  the  penetrator  results  from  particles  flowing  opposite 
the  direction  of  penetration.  Fragments  of  the  penetrator  fretted  from  the  leading  surface 
generally  have  an  initial  forward  velocity  component  and  may  travel  into  open  cone  and 
radial  cracks  ahead  of  the  tip  of  the  penetrator.  Other  penetrator  fragments  mix  with  and 
flow  with  the  ceramic  powder,  becoming  part  of  the  front  surface  ejecta.  The  eroded 
tungsten  fragments  exhibit  greatly  elongated  grains  close  to  the  fragment  surfaces,  indicative 
of  heavy  localized  plastic  flow. 


CONCLUSIONS 

According  to  this  concept  of  penetration  phenomenology,  the  properties  of  a  ceramic 
that  govern  penetration  resistance  include  the  compressive  strength  and  hardness,  the 
pulverization  characteristics,  the  frictional  flow  characteristics  of  fine  fragments,  and 
fragment  abrasiveness.  These  properties  are  consistent  with  those  suggested  by  Meseall 
[13,14]. 

Initial  resistance  to  penetration  is  provided  by  the  compressive  strength  or  hardness  of 
a  ceramic.  High  compressive  strength  is  desirable  to  deform,  fracture  and  deflect  an 
impacting  body.  Projectiles  with  low  aspect  ratios  can  be  defeated  if  the  strength  of  the 
ceramic  exceeds  the  strength  of  the  projectile.  High  aspect  ratio  projectiles  such  as  long 
rods,  may  suffer  heavy  deformation  and  fracture  damage  to  the  proximal  end,  but  the  intact 
distal  portion  will  continue  to  advance  and  penetrate  the  ceramic.  Thus,  a  high  ceramic 
compressive  strength  can  resist  penetration  only  to  a  certain  extent. 

The  stresses  exerted  by  a  long  impacting  rod  will  eventually  pulverize  the  ceramic  material 
in  a  small  zone  immediately  ahead  of  the  leading  surface  of  the  penetrator.  As  explained 
in  the  following  paragraph,  production  of  a  pulverized  zone  is  a  necessary  condition  for 
a  penetrator  to  advance  in  a  confined  ceramic  target.  Thus,  resistance  to  comminution 
is  desirable  for  penetration  resistance.  Although  the  specific  fracture  surface  energy  for 
most  ceramics  is  small  in  tension,  the  energy  required  to  produce  a  unit  of  failure  surface 
area  under  large  dynamic  compressive  and  shear  forces  may  be  significantly  greater.  Thus, 
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the  energy  absorbed  in  dealing  th?  surface  area  of  the  powder  may  be  a  significant  ceramic 
property  for  penetration  resistance. 

A  penctrator  can  only  advance  if  the  material  in  its  path  is  pushed  ahead  of  it  or  to  the 
side.  Because  of  heavy  rear  confinement,  the  crushed  ceramte  cannot  he  pushed  ahead  and 
out  the  rear  surface  in  the  way  that  metallic  armor  plates  fail  by  plugging.  And  if  the 
ceramic  is  nonporous  and  snugly  confined  laterally,  the  pulverized  material  cannot  be 
pushed  to  the  side.  Indeed,  the  only  recourse  is  for  the  powder  to  flow  opposite  the 
penetration  direction  along  the  cavity  being  produced  by  the  penctrator.  Thus,  the 
pulverized  and  diluted  ceramic  must  flow  under  high  pressure,  and  so  the  frictional  flow 
property  of  the  comminuted  ceramic  should  influence  penci  ration  resistance.  We  expect 
this  property  to  depend  on  pressure,  strain  rate  and  size  and  shape  of  the  fragments,  and 
to  be  describable  bv  a  Mohr-Coulomb  type  curve.  For  thick,  highly  confined  ceramic 
blocks,  the  friction-flow  property  of  ceramic  fragments  may  be  the  most  important  material 
property  for  penetration  resistance. 

Finally,  the  ahilitv  of  a  ceramic  to  erode  a  penetrating  rod  is  a  desirable  property  for 
penetration  resistance.  Whereas  erosion  may  be  by  gross  local  plastic  flow  of  the  leading 
rod  surface,  ceramic  fragments  that  gouge,  score,  shear  or  otherwise  abrade  the  rod 
material  may  reduce  the  incoming  mass  and  terminate  the  penetration  earlier  than 
non-abrasive  target  materials.  Wear  and  erosion  can  occur  by  a  number  of  mechanisms 
depending  on  penetrator  and  target  material,  fragment  geometry  and  size,  pressure, 
temperature,  and  flow  rate.  Thus,  optimal  erosive  behavior  might  be  achievaole  by  matching 
the  abrasive  characteristics  of  a  ceramic  material  to  the  threat. 

The  fractographic  observations  and  the  deduced  penetration  phenomenology  reported 
here  can  also  be  used  to  identify  properties  governing  the  penetration  capability  of  rods. 
To  be  effective  as  a  penetrator,  a  material  should  have  high  density  to  produce  high  stresses 
in  the  target,  a  high  yield  strength  to  resist  mushrooming  at  the  leading  edge,  a  high  work 
hardening  rate  to  suppress  the  tendency  to  shear  band  and  fret,  a  high  fracture  toughness 
to  resist  the  propensity  for  rod  shaft  failure,  and  high  abrasion  resistance  to  resist  scoring 
and  erosion  by  ceramic  particles. 

In  future  work,  this  understanding  of  penetration  phenomenology  will  be  used  to  develop 
tests  that  measure  dynamic  shear  strength  and  flow  resistance  of  intact  and  fragmented 
ceramic  material  under  high  confining  pressure,  and  to  develop  computational  models  of 
penetration  that  can  be  used  to  assist  in  designing  ceramic  armor. 
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Ceramics  and  brittle  geologic  materials  can  deform  nonelastically  under  compression  and  shear  by  sliding  and  ride-up  of 
fragments,  grains,  or  blocks,  with  accompanying  competition  between  dilatancy  and  pore  compaction.  This  paper 
describes  a  micromcchanical  model  of  such  deformation  based  on  an  analogy  to  the  dynamics  of  atomic  dislocations  and 
using  the  slipping  processes  in  a  multiple- plane  plasticity  model.  In  the  analogy,  atomic  lattice  slip  planes  are  replaced  by 
slip  on  interfaces  between  blocks,  and  atoms  arc  replaced  by  fragments,  grains,  or  blocks.  The  model  treats  intact 
material,  a  gradual  transition  to  fragments,  and  then  the  behavior  of  the  fragments. 

The  coefficient  of  friction  between  the  grains  (tan  $  =  0. 1 1)  was  determined  for  aluminum  nitride  through  simulations  of  a 
pressure-shear  impact  experiment  The  calibrated  model  was  then  used  with  a  finite-element  computation  of  the 
penetration  of  a  long  tungsten  rod  into  an  aluminum  nitride  target  The  calculations  successfully  exhibited  the  key  features 
observed  experimentally,  such  as  the  rubblized  rone  just  ahead  of  the  penetrator  and  the  dilatant  rebound  of  the  ceramic 
behind  the  penetrator.  The  most  important  microscopic  material  property  governing  depth  of  penetration  was  found  to  be 
the  friction  between  granules  in  the  rubble  ahead  of  the  penetrator. 


1.  INTRODUCTION 

Hard,  brittle  materials  often  deform  nonelastically  by 
fracture,  fragmentation,  and  subsequent  frictional  sliding  and 
rotation  of  fragments.  The  sliding  produces  additional  ride- 
up  friction  and  associated  dilatancy.  The  dilatancy  ts 
counteracted  by  pore  compaction  under  confining  pressure. 
These  processes  are  especially  evident  in  the  eroding 
’process  zone”  or  ’’Mcscall  zone”  at  the  nose  of  a  long 
metallic  rod  penetrating  a  ceramic  armor.  Another  example 
is  the  ’’block  motion” in  faulted  geologic  materials  exposed 
to  ground  shocks  produced  by  earthquakes  or  nuclear 
weapons.  Hard  sands  also  deform  by  grain  sliding  or 
rolling,  ride-up,  and  compaction. 

We  developed  a  micromcchanical  model  for  motion  of 
granular  beds  and  calibrated  it  by  computationally  simulating 
a  pressure-shear  gas-gun  impact.  We  then  used  this  model 
to  compute  the  penetration  of  a  Jong  rod  into  ceramic. 

Attempts  to  develop  a  micromechanical  model  for  such 
behavior  usually  focus  on  building  continuum  constitutive 
relations  that  average  microscopic  fracture  or  grain 
movement  processes.1'5  For  our  model,  we  described  the 
motion  of  lines  ot  voids  between  the  moving  objects,  as  in 
atomic  dislocation  theory.  Our  model  combines  fracture  and 
slip-induced  softening  processes,  compressive  behavior 


(including  pore  compaction),  and  dilatancy. 

2.  MODEL  DEVELOPMENT 

In  developing  the  model,  we  visualized  an  aggregate  of 
blocks  or  grains  that  arc  initially  tightly  packed.  When  a 
large  shear  strain  is  applied  to  this  aggregate,  grains  may 
slide  or  roll  over  each  other,  but  as  they  move,  they  develop 
voids  because  the  material  cannot  remain  tightly  packed 
during  the  motion.  The  slipping  can  occur  only  on  the 
available  planes  of  the  blocks  and  not  on  the  orientation'*  of 
maximum  shear  stress.  These  considerations  led  us  to  a 
model  that  is  based  on  the  Orowan6  equation  relating  plastic 
strain  rate  to  the  number,  size,  and  velocity  of  dislocations. 

In  our  case  the  plastic  strain  is  the  inelastic  shearing 
deformation  associated  with  sliding  (rolling  is  not  treated), 
and  the  dislocations  are  the  blocks  or  grains.  The  basic 
relation,  which  is  derived  elsewhere,7  is 

dyf/dt  =  (Tj  -  HOnj)/(PsCiT)  (1) 

P 

where  Yi  is  the  plastic  strain  on  the  ith  plane,  tj  is  the  shear 
stress,  p  is  the  tangent  of  the  friction  angle,  C™  is  the  normal 
stress  on  the  plane,  ps  is  the  solid  density.  Q  is  the  speed  of 
transverse  (shear)  waves,  and  T  =  B/(2gfj  Q),  a  time 
constant.  B  is  the  block  dimension,  g  is  a  factor  on  the  order 


of  1  relating  die  active  slip  system  to  the  shear  strain 
direction*  and  f*  is  the  fraction  of  blocks  that  have  an 
adjacent  hole  large  enough  to  allow  slip  into  it 

Next  we  note  that  the  porosity  ($)  increases  in  proportion 
to  the  plastic  shearing  strain  (a  dilatant  process)  and 
decreases  by  compaction.  We  write  the  dilatancy  expression 
as 

d$/dt  =  a  dyf/dt  +  d^dt  (2) 

where  d<(>c/dt  is  a  compaction  rate  and  a  is  a  dilatancy 
coefficient  on  the  order  of  1 .  To  define  the  compaction 
process,  we  used  the  compaction  curve  developed  by  Carroll 
and  Holt 8  The  curve  is  given  by  an  implicit  relation 
between  pressure  and  density.  With  the  formulation 
represented  by  Eq.  (2)  and  the  compaction  curve,  when  we 
shear  the  material  at  constant  volume  and  cause  plastic 
shearing  strain,  the  pressure  rises  in  accordance  with 
Eq.  (2).  When  this  pressure  reaches  the  compaction  curve, 
the  pressure  ceases  to  rise  and  the  material  maintains  that 
porosity. 

Eqs.  (1)  and  (2)  are  solved  with  die  elastic-plastic 
deviator  stress  relations  and  a  Mie-Griineisen  relation  for 
pressure  to  determine  the  stress  state. 

The  preceding  constitutive  model  for  granular  flow  was 
incorporated  into  a  multiple -plane  plasticity  model  for 
computational  simulations  of  dynamic  loading  problems. 

The  computational  model  (termed  FRAG  BED)  has  the 
following  features: 

♦  Nine  planes  on  which  ad  slip  and  fracture  occur,  as 
seen  in  Figure  1, 

♦  The  granular  flow  and  dilatancy  processes,  porosity 
model,  and  Mic<rriineisen  relation  mentioned  above. 

♦  Yield  and  tensile  strengths  that  are  reduced  by  the 
developing  fragmentation. 

3 .  CORRELATION  OF  THE  MODEL  PREDICTIONS 

WITH  EXPERIMENTAL  DATA 

To  test  and  calibrate  the  FRAGBED  model  for  ceramic 
material,  we  made  computational  simulations  to  represent 
several  laboratory  experiments  for  AIN.  The  static  tests  of 
Heard  and  Cline9  and  the  dynamic  tests  of  Lankford10  were 


FIGURE  1 

Initial  orientations  of  the  slip  planes 
in  the  FRAGBED  model. 

simulated  using  the  intact  properties  of  the  AIN  with  a 
Coulomb  yield  strength  of  p  *  0.33,  density  of  3.26  Mg/m3, 
bulk  modulus  of  210  GPa,  initial  yield  of  3.0  GPa,  and  a 
shear  modulus  of  1 26  GPa. 

A  more  detailed  study  was  made  of  the  pressure-shear 
impact  test  of  S hockey  and  Klopp. 1 1  The  geometry  of  the 
experiment  is  shown  in  Figure  2.  A  laser  interferometric 
transverse  velocity  measurement  at  the  rear  surface  of  the 
boron  carbide  anvil  was  interpreted  to  show  a  constant 
normal  stress  of  10.9  GPa  and  a  nearly  constant  axy  in  the 
AIN  of  L78  GPa.  Several  one-dimensional  wave 
propagation  simulations  (with  the  FRAGBED  model)  were 
made  of  this  experiment  to  obtain  appropriate  parameters. 


impactor  experiment. 


Our  results  indicate  Lh  at,  dunng  the  passage  of  the  initial 
compressive  wave,  planes  2, 4,  and  5  (all  at  45°  to  the 
loading  axis,  as  shown  in  Figure  1)  received  sufficiert 
shearing  strain  to  eliminate  their  cohesive  strength.  The 
renewed  loading  that  occurred  during  the  passage  of  the 
shear  wave  brought  the  shear  stress  on  the  fifth  plane  (45° 
between  X  and  Z)  to  yield  again,  but  did  not  cause  slip  on 
any  other  plane;  hence,  this  plane  controlled  the  subsequent 
yielding  and  flow.  With  the  coefficient  of  friction  of 
H  *=■  0.1 1  (a  surprisingly  low  value),  the  model  represented 
well  the  measured  motions.  The  computed  normal  and  shear 
stress  and  the  yield  strength  on  the  plane  are  shown  in 
Figure  3. 
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FIGURE  3 

Normal  and  shear  stress,  yield  strength,  and  plastic 
deformation  on  plane  5  in  slanted  impactor  test  in 
aluminum  nitride. 

4 .  PENETRATION  CALCULATIONS 

Penetration  calculations  were  performed  to  study  the 
ceramic  behavior  near  the  head  of  the  penctrator.  In  the  test 
configuration  shown  in  Figure  4,  the  tungsten  penetrator  has 
a  Icngth’todiamcter  ratio  of  10  and  impacts  at  1 .6  km/s. 

One  sample  of  the  computed  results  is  shown  in  Figure  5. 

At  80  jis  the  penctrator  has  been  heavily  deformed  and 
nearly  stopped.  Surrounding  the  penetrator  is  a  rubblized 
zone,  partially  fractured  zone,  and  an  intact  zone,  The 
contours  labeled  '  fracturing  material"  indicate  the  transition 


from  intact  to  fully  fragmented.  Here  it  is  clear  that  the 
model  is  not  following  indmdual  large  cracks  (which  may 
occur  in  actual  penetrations),  but  represents  only  a  general 
rubbltzation  of  the  material. 
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FIGURE  4 

Geometry  of  penetration  calculations. 

The  locations  and  sizes  of  these  zones  correspond  well 
with  observations  of  ceramic  targets  following  penetration 
tests.  Behind  the  penetrator  (top  of  the  figure),  the 
granulated  ceramic  closes  up  the  hole  made  by  the  penetrator 
(as  observed).  The  flume  of  highly  expanded  ceramic 
powder  that  rose  behind  the  penctrator  is  not  shown.  The 
most  important  microscopic  material  property  governing 
depth  of  penetration  was  found  to  be  the  friction  between 
granules  in  the  rubble  ahead  of  the  penetrator. 
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FIGURE? 

Computational  result  at  80  ps  for  penetration  into  AIN, 
showing  intact,  fracturing,  and  rubble  regions. 

5.  CONCLUSIONS 

The  micromechanical  model  describing  ihe  processes  of 
intergranular  sliding,  dilatancy,  and  compaction  appears  to 
be  appropriate  and  useful  for  describing  the  highly 
nonlinear,  anisotropic,  and  history-dependent  deformation  of 
fragmented  beds.  For  the  case  of  long  rods  penetrating 
ceramic  armor,  the  model  predicts  that  the  intergranular 
friedon  governs  the  penetration  depth.  This  result  suggests 
that  experiments  such  as  the  pressure-shear  impact  test  of 
Shockey  and  KJopp1 1  should  be  used  to  examine  candidate 
ceramics  for  armor  applications. 
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Analysis  of  tilt  in  the  high-strain-rate  pressure-shear  plate  impact  experiment 
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An  analysis  is  presented  of  the  effect  of  tilt  in  the  high-strain-rate  pressure-shear  plate  impact 
experiment.  Normal  and  shear  tractions  and  the  velocity  difference  across  the  specimen  are 
obtained  from  frec-surface  particle  velocity,  tilt  angle,  closure  direction,  and  impact  velocity 
data.  The  effect  of  tilt  in  high-strain-rate  pressure-shear  plate  impact  is  shown  to  be  small 
unless  the  velocity  difference  across  the  specimen  is  itself  small. 


The  effect  of  tilt  in  experiments  involving  either  normal 
impact  or  pi  ess u re-shear  impact  of  two  parallel  plates  has 
been  analyzed.1,2  No  results  have  been  presented  however, 
on  the  effect  of  tilt  in  the  high-strain-rate  pressure-shear 
plate  impact  experiment,  shown  in  Fig.  l/  in  which  a  thin- 
specimen  plate  is  compressed  and  sheared  between  hard, 
clastic  flyer  and  target  plates,  and  the  stress  and  strain  rate  in 
the  specimen  arc  inferred  from  particle  velocity  measure¬ 
ments  at  the  back  of  the  target.  Because  of  the  tilt,  the  waves 
generated  are  inclined  with  respect  to  the  plate  surfaces.  In 
this  communication,  an  analysis  is  presented  of  the  effect  of 
the  inclination  of  these  waves  on  the  calculation  of  stress, 
strain,  and  strain  rate  in  the  specimen. 

Figure  2(a)  shows  the  pattern  of  waves  generated  dur¬ 
ing  tilted  impact.  The  specimen  is  idealized  as  a  slidi  ng  inter¬ 
face  of  zero  thickness.  The  flyer  and  target  arc  taken  to  be 
linear  clastic.  The  x,  direction  is  normal  to  the  impact  plane, 
and  the  x?  direction  is  in  the  direction  of  closure.  The  angle 
from  thex^  direction  to  the  shearing  direction  is  denoted  by 
fl  (Fig  2(b)).  With  an  assumed  supersonic  closure,  a  nor¬ 
mal  {p)  wave  is  generated  in  the  target  at  angle  t/',  and  two 
shear  (sv,sh)  waves  as  generated  at  angle  VV  For  small  an¬ 
gles,  t//(  —  ac,/u0,  where  a  is  the  tilt  angle,  ct  is  the  wave 
speed,  and  u()  is  the  component  of  impact  velocity  norma]  to 
the  plates.  For  a  typical  pressure-shear  experiment, 
a  =  10“ 3  rad,  w0  =  200  m/s,  and  c,  —  6000  m/s,  and  sc 
^,^30X  10~?  rads  1.7*. 

Jn  the  interpreiation  of  the  experiments,  the  measured 
normal  and  transverse  particle  velocities  in  region  9  arc  re¬ 
lated  to  the  stresses  and  particle  velocities  in  regions  3  and  4. 
Let  p''be  the  unit  polarization  vector  of  the  wave  traveling 
in  to  region  ;  from  region/  Let  [  u],J  designate  the  jump  in  u 
from  regions  \  toy:  u‘  —  u\  The  usual  subscripts  des¬ 

ignate  vector  and  tensor  components.  The  particle  velocity  is 
v  and  the  stress  is  u.  The  jump  in  particle  velocity  across  a 


wave  front  can  be  expressed  ns  follows: 

[vr'HM'  (1) 

where  A  >  J  are  amplitudes  to  be  determined.  For  isotropic 
materials,  the  polarization  vectors  are  normal  to  the  wave 
front  for p  waves  propagating  at  speed  c,  and  parallel  to  the 
wave  front  for  sv  and  sh  waves  propagating  at  speed  c:.  Fora 
wave  propagating  in  the  direction  n  in  the  xrx2  plane  at  a 
counterclockwise  angle#  from  the  positivex,  axis,  the  three 
polarizations  are 


'cos  6' 

r  —  sin  #' 

o' 

sin  6 

-H*  _  , 

1  p  — 

cos  d 

0 

.  o  j 

L  o  J 

k 

for  the  p  w'ave,  the  sv  wave,  and  the  sh  wave,  respectively. 

The  respective  notarization  vectors  for  the  stresses  are 
given  by 


Laser  Interferometer 


FIG.  1 .  Schematic  of  lhe  high-stram-raic  pressure-sheur  plate  impact  exper 
imenl 
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FIG.  2  Schematic  of  tilted  impact  of  inclined  plates.  (a)  Wave  pan em  gen- 
eraied.  lb)  Orientation  of  closure  direction  relative  to  shearing  direction. 
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where  corresponds  to  Icrn.t722,(r}j,t72itcrl3,crl2]r, 

p  is  the  Oyer  and  target  density,  and  77  =  c2/cv 

The  stress  and  particle  velocity  in  the  various  regions  in 
Fig.  2(a)  can  be  expressed  as  the  sum  of  jumps  across  succes¬ 
sive  wave  fronts.  Thus 


v3  —  v1  +  Apl2jtpL1  +  A  Iu2  Jp,1,2'J 


(4) 


V4  =  V6  +  A  '’‘■y6  5  +  A  “Ayj5'‘  +  A  >h 5  V* 5'4. 

The  surface  tractions  on  the  two  sides  of  the  specimen  are 
given  by  similar  relations,  with  <7, , ,  <7, 2,  and  a,  3  replacing  u,, 
u,,and  v3  ands„r6,  and s}  replacing p,,p2, and p},  respective¬ 
ly 

Tractions  are  expected  to  be  continuous  across  the 
boundary  between  regions  3  and  4  once  sufficient  wave  re¬ 
flections  have  occurred  to  establish  a  homogeneous  stress 
state  in  the  specimen,  The  flyer  and  target  are  initially  stress 
free  and  the  target  is  at  rest.  Thus 

A  pl2sf!2  +  A  w2,V/’2'1  +  A ,h 2'V,*2'3  -  A  5 

-  A  ‘^'V5'4  -  A  M5'4jf 5'4  =  0 

A  +  A  4  A  ,h  lV"']  -  A 

~  A  w>Asfy*~  A  ,h  s  *siH  3  4  =  0  ( 5 ) 

A  r'-2#'  2  +  A  "“if43  +A,h  13s‘,h2  >  -  A 

-  A  "'5  *s?iA  -A‘”  y*sf  y*  =  0. 

The  normal  velocity  and  the  transverse  velocity  perpendicu¬ 
lar  to  the  shearing  direction  are  expected  to  be  continuous 
across  the  specimen,  Thus 

Apy2pf-2  +  A  "^pf*  +  A 

-  A  p65pf?  -  A  aiApf*  -A*  y*p\l,i* 

-  v*  =  0,  (6) 


_ I 

A  '“■2(P?1 2  cos  n  -  pf-2  sin  11)  +  A  u,2  \p?  i  cos  Q  -  p\v2  i  sin  11 )  +  A  * yi(p\h2i  cos  11  -  p\h  2'3  sin  11 ) 

-  A  cos  a  -  pf 5  sin  11) -A  uiHpfA  cos  11  -  pf 4  sin  11)  ~A  >hy\pf 5  4  cos  fl  -  p\h  5*  sin  Cl ) 

—  u*  cos  ft  +  sin  11  =  0.  (7) 


The  velocity  in  the  shearing  direction  is  expected  to  have  a  jump  across  the  specimen  of  yh,  where  y  is  the  shearing  rate  to  be 
determined  and  h  is  the  specimen  thickness.  Thus 

A  pl2(jf2' 2  sin  11  +  pf-2  cos  11)  +  A  “2\ pf 3  sin  11  +  pf1  cos  11  )+A  ,l,21(p? 2i  sin  11  +  pf  2  3  cos  11 ) 

-  A  <*■ 5  {pf 5  sin  11  +  pf' 1  cos  11 )  -  A  ^A(pf*  sinfl  +  pf 4  cos  11 )  -  A  *  y\p'2h  y*  sin  11  +  pf 54  cos  11 ) 

—  v\  sin  11  —  u)  cos  11  =  yh.  (8) 

Equations  (5)-(8)  relating  tractions  and  particle  velocities  on  the  two  sides  of  the  specimen  provide  six  equations  for  relating 
the  six  unknown  amplitudes  Apl  2,A  'v2i,A  ,h 2,\  A  pb',t  A  "*■*,  and  A  'hiA,  and  the  unknown  velocity  jump  yh. 

A  seventh  equation  is  obtained  by  relating  the  wave  amplitudes  to  the  free-surface  velocity  measured  in  region  9; 
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T9  =  v1  +  A  '' VU  +  a  "Wp"1-3  +  A  +  A  +  A  ,v8-V8-7  +  A  "’8-7p"'8'7  +  A  *‘My  uo 

_j_  ^  WIO.»pMil0.9  ^  lA  la9p!'1  ,0>9.  (9) 

Equation  (9 )  contains  six  additional  amplitudes,  A  p2,8,  A  ,u8,7,  A  ih*'\  A  plU0,  A  ™10*9  and  A  lh  l0'\  which  can  be  eliminated  by 
considering  diffraction  at  the  target  free  surface.4  Diffraction  of  the  1,2  p  wave  from  the  target  free  surface  provides  three 
relationsfor  the  amplitudes  A  p2*fA  SV*-\A  ih*'\  and  A  p1,2.  Similarly,  diffraction  of  the  2,3  ju  and  sh  waves  from  the  rear  surface 
provides  three  relations  for  the  amplitudes  A  pllJ0,  A  °*9,  A  l/Tl0,9,  A  n,1'\  and  A  Ih2'\  For  a^sina,  application  of  these 
relations  to  Eq.  (9)  gives 

AP\.*^ppui  sjn  n  H-/?5,,2cos  H  —  P?2'8  sin  ft  —  p?'*  cos  H  —  4^,  (p^8,7  sin  H  +  p?*J  cos  ft)  ] 

+  Sinn  -h^5;2*3  cos  n  +477^(/7flU0  Sinn  +/7plut0  cos  H)  —  /j™10-9  sin  n  —  /?5U,°9  cos  H  ] 

+  A thl* [p?lm*  sin  n  4- p*}2* cos  n  + ,0’9  sin  n  +piMMcos  n  J  =  v\  sin  n  4-  v]  cosH,  (10) 


where  v]  sin  fi  +  tij  cos  il  is  the  measured  transverse  veloc¬ 
ity.  FromEqs.  (5)-(8)  and  ( 10),  yh  and  the  unknown  am¬ 
plitudes  can  be  obtained  from  the  measured  velocity  compo¬ 
nents  in  region  9  and  the  projectile  velocity  of  region  6.  The 
tractions  on  the  specimen  can  then  be  evaluated  from  rela¬ 
tions  analogous  to  those  of  Eqs.  (4). 

As  an  illustration  of  the  above  analysis,  consider  an  ex¬ 
ample  in  which  an  AIN  specimen  ( thickness  h  —  0,577  mm) 
is  impacted  between  B4C  flyer  and  target  plates  (p  =  2.51  g/ 
cm3,  cx  =  13.8  mm//is,  c2  =  8.5  mm//ts)  at  a  skew  angle  of 
0  »  20*.  The  impact  velocity  is  VQ  =  672  m/s,  and  the  tilt  is 
a  =  4  mrad.  Closure  occurs  nominally  in  the  shearing  direc¬ 
tion,  with  n  =97°.  The  average  target  frec-surface  trans¬ 
verse  velocity  is  167  m/s.  From  these  data,  the  wave  ampli¬ 
tudes  and  yh  can  be  obtained  from  Eqs.  (5)-(8)  and  (10), 
resulting  :n  a  shearing  rate  of  1. 1 3  X  105  s~‘.  From  the  am¬ 
plitudes,  one  calculates  a  specimen  normal  stress  of 
—  10.938  GPa  and  a  specimen  shear  stress  of  1.776  GPa.  If 
no  tilt  correction  is  made,5  then  one  obtains  a  normal  stress 
of  —  10.945  GPa,  a  shear  stress  of  1.781  GPa,  and  a  shearing 
rate  of  1 .09  X  10s  s~ 1 .  The  maximum  error  of  approximately 
3.5%  is  in  the  shearing  rate. 

The  magnitude  of  the  error  depends  on  n.  If  H  is  0*  so 
that  the  shearing  is  normal  to  the  closure  direction,  the  mo* 
tion  used  to  derive  the  shear  stress  and  shearing  rate  is  car¬ 
ried  bys/?  waves  uncoupled  from  the  p  waves.  Thus  a  tilted 
normal  wave  has  little  effect  on  the  shear  stress  and  shearing 
rate.  When  FI  is  near  ±  90°,  most  of  the  motion  used  to 
derive  the  shear  stress  and  shearing  rate  is  carried  by  su 
waves  ’  ’oled  to  the  p  waves  by  diffraction.  In  this  case,  the 
error  L  .ill  small  but  is  larger  than  in  the  FT  =  0°  case.  If  one 
sets  H  =  O’  in  the  example,  then  one  obtains  a  shear  stress  of 
1.779  GPa  from  which  the  uncorrected  value  of  1.781  GPa 
differs  by  only  0.1%,  and  the  corrected  and  uncorrected 
shearing  rates  are  identical. 

Theerrors  are  small,  even  in  large-tilt  case  considered  in 
the  example,  because  the  excess  x2  component  of  the  free- 
surface  velocity  is  th  e  same  as  the  excess  xz  component  of  the 
flyer  velocity  that  is  approximately  the  x ,  component  of  flyer 


I - 

velocity  uf  multiplied  by  the  tilt  angle.  The  contribution  of 
av\  to  the  free-surface  transverse  velocity  is  usually  small.  In 
cases  in  which  the  shearing  rate  is  low,  however,  contribu¬ 
tions  of  ow*  to  the  free-surface  transverse  velocity 
=  ul  sin  H  +  v\  cos  H  can  lead  to  large  relative  errors  in 
the  shearing  rate.  This  effect  can  be  made  evident  by  consid¬ 
ering  the  fractional  derivative  of  the  equation  for  the  shear¬ 
ing  rate  which,  for  a  —  0,  becomes 

dy  _  —  dv f< 

r  "o  -  ' 

where  u0  —  sin  ft  +  v*  cos  ft.  If  thei>0  —  vk  is  small,  dy/ 
y  could  be  large. 

The  tilt  error  in  high-strain-rate  pressure-shear  plate 
impact  has  been  shown  to  be  small  except  when  the  shearing 
rate  is  low.  The  above  example  shows  that  the  excess  target 
free-surface  transverse  motion  induced  by  the  tilt  is  very 
small.  Because  the  transverse  momentum  resulting  from  fly¬ 
er  tilt  is  very  small,  the  small  effect  of  the  tilt  on  the  trans¬ 
verse  motion  of  the  target  free  surface  is  not  unexpected. 
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5Onc  way  is  to  simply  let  it>,  —  0  in  Eq.  (5)-(8)  and  ( 10). 
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